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HAUPTAUFSÄTZE 


Theorie des Resonanzschwingungsdämpfers. 
Von ©. Föppl in Braunschweig, Wöhler-Institut. 


n Abb. 1 sitzt auf einem Balken ce, ein Motor mit der umlaufenden Masse m,. Wenn der 

Motor mit der bestimmten Eigenschwingungszahl », der Anordnung e,, m, umläuft, und 
der Rotor etwas exzentrisch gelagert ist, macht e, besonders große Ausschläge (Resonanz- 
schwingungen). Man kann diese Ausschläge wesentlich vermindern, wenn man einen Dämpfer 
c, m, von etwa der gleichen Eigenschwingungszahl », = n, anbringt, der infolge der großen 
Nachgiebigkeit der Feder ce, zu großen Ausschlägen £&, der Masse m, aufgeschaukelt wird und 
der etwa durch Reibung R=k Se Energie vernichtet. Resonanzschwingungsdämpfer dieser 
Art haben heute schon eine große Bedeutung für die Abdämpfung der Drehschwingungen 
von Diesel-, Auto- oder Flugzeugmotoren-Kurbelwellen. Man wird ferner solche Dämpfer mit 
großem Vorteil verwenden können, um Biegungsschwingungen etwa von Flugzeugen zu 
dämpfen. Im Nachfolgenden wird auf diese Einzelheiten der Praxis nicht eingegangen, 
sondern nur ein möglichst einfaches Beispiel behandelt. 




















C, 

| tg 7 u | pP 

die A , S 
mM = 
’ I ZW; 

= Cy Cz 

[RA25121) [RA2S1Z2] — &— -—-£&, — 
Abb. 1. Abb. 2. 


Abb. 2 zeigt die schematische Anordnung, die wir der Betrachtung zugrunde legen. Die 


Hauptanordnung e,, m, hat die Eigenschwingungszahl », 


Vı/ce 
| \ ', wenn ce, zugleich auch 


j m, 

die Federungskonstante angibt. Der Dämpfer e, m, ist auf e, m, abgestimmt‘). Wenn m, 

sroß ist gegen m,, ist angenähert — = = . An m, greift die erregende Kraft P= P,cos(iwt - a) 
? 2 Mm, : 


- 


1) 0. Föpp!l: Ing.-Arch. 1931, S. 347. 
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an, die genau im Tempo der Eigenschwingungszahl der Gesamtanordnung wirkt. P ist al: 
Resonanz erregende Kraft um W" gegen £, phasenversetzt. Wir untersuchen den Beharrungs 
zustand. Es ist &, = &,sin(»f—a), wenn mit &, der Größtausschlag der Masse m, be 
zeiehnet ist. An m, wirkt ferner etwa durch Scheuern an der mit m, fest verbundene: 

d(£&,—£,) 


Stange s die Reibungskraft R=k Fr 
( 
&, 5, £&, der beiden Massen gegeneinander ist. Statt der von der Zeit abhängigen Weg: 


5, &,sin(ot a) und &,— £,, sin (ot — 5) verwenden wir im nachfolgenden die Werte : 
und &,. Es ist also &,=£&,-+2.. 


‚ die verhältnisgleich der Relativgeschwindigkeii 


In Abb. 3 ist ein Vektordiagramm wiedergegeben, in dem die Vektoren £,, und &,, ein 
getragen sind. Die Indizes Null geben an, daß die Größtwerte gemeint sind. Das Vektor 

diagramm dreht sich mit der Winkelgeschwindigkeit &. Di: 

augenblicklichen Werte £, und £&, erhält man durch Projektion 

der Vektoren auf die Horizontale. Die Zeit t=0 falle zusammen 
Aw&,, mit &,—=0, so daß &,=£,,: sin wi ist. 







Infolge der Reibung wirken die Massen m, und m, auf- 
einander mit einer Kraft, die gleich ist kwE£,,cosot und 
deren Vektor auf dem Vektor £,, senkrecht steht (Abb. 3). 
Wir zerlegen die Vektoren in Abb. 3 in eine Bewegung 4, 
die mit &,, gleichphasig ist, und in eine Bewegung DB, die 
um 90° gegen A phasenversetzt ist. Den Größtausschlag £, 
der Masse m, zerlegen wir entsprechend in die beiden Kom- 
ponenten &,, und &,,, von denen &,,, In Richtung der Be- 
wegung A und &,,, In Richtung B liegt. 

Wir betrachten den Beharrungszustand, der bei einer bestimmten Reibung an der 
Schwingungsanordnung nach Abb. 2 auftritt. Im Beharrungszustand sind die dynamischen 
Grundgleiehungen für die Bewegungen A und B je für sich erfüllt. Wir können also je zwei 
Gleichungen für die Massen m, und m, nach den beiden Richtungen A und B aufstellen, bei 
denen die Kräfte gleich den auftretenden Beschleunigungen gesetzt werden. Wir beachten, 


[RA25123 


Abb. 3. 


a d? E, . Pa Y (\,-p . 
daß der Größtwert der Beschleunigung Ta gleich &,,: ©” und der Größtwert der Reibungs- 
8 1R a 
D dE, . > . U. . 
kraft % j eleich ko £,, Ist und erhalten für die Masse m;,: 
( 

m,(&ro + S10) © — 0, &,, =0 (Bewegung A) . . . . 2.2...) 
m, —KkE,,® 0 (Bewegung B) . . . 2.0...) 


und für die Masse m: 


. 9 = “ Sı2 P 
m, + ec Er P.—-=0 (Bewegung A) . : 2.2.2.2... ß8) 
10 
= ni) = Sı10 
DE FRE nl FE Me 2 IT’ Be u Tg 0 (Bewegung BD) . . .. 2.0.4). 
>10 


Gleichung 1 bezieht sich auf die Masse m,, deren Größtweg in Richtung A gleich &,,+ &,, 
ist. Die Federkraft, herrührend von der Feder e,, ist verhältnisgleich mit &,. Die Bewegung I 
der Masse m, (Gl. (2)) ist gleich &,,,:c0osot#t. In dieser Richtung wirkt die von der Reibung 
herrührende Kraft, die gleich £,,:k-w-cos mE ist. 

An der Masse m, greift sowohl die Federkraft e,. die verhältnisgleich &,,, bzw. &,,, Ist. 
als auch die Federkraft e, an. Der Ausdruck e,£,, tritt allerdings nur bei der Bewegung 4 
auf (Gl. (3)), während die von der Reibung herrührende Kraft nur die Bewegung B beeinflußt 
(Gl. (4)). Endlich ıst in den Gl. (3) und (4) noch die äußere Kraft P zu berücksichtigen, die 
gegen &, um 9%" phasenversetzt ist. Die Kraftkomponente für die Bewegungen A ist deshalb 
verhältnisgleich der Wegkomponente £,,, in Riehtung B und umgekehrt. Wir hätten noch 
eine fünfte Gleichung aufstellen können, die angibt, daß die Je Schwingung im Dämpfer dureh 
Reibung umgesetzte Energie ak: © £,,’ gleich der zugeführten Energie 7 P,£,, ist. Durch diese 
Gleichung wird P’ in eine lineare Abhängigkeit von k gebracht, die wir als gegeben annehmen 
wollen. 

Mit Hilfe der vier Gleichungen können wir für einen gegebenen Fall sämtliche Größen 
bestimmen. Nehmen wir z. B. an, es sei £,, gegeben, so können wir die Größen &,,0» &ss: 
k und «» berechnen. Eine dieser vier Größen, nämlich ®, können wir mit sehr guter An- 


) 
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Py 
u 
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‚äherung, für den Fall daß m, klein ist gegen m, und nur dieser Fall hat große praktische 
Bedeutung — , berechnen. Es ist nämlich: 
15) 0) C, - 
a ae he te Far 5 
Mm; 


; ist dabei ein Wert, der mit guter Annäherung gleich 1 gesetzt werden kann. Der genaue 
Wert von a hängt von der Abstimmung des Dämpfers ab. Es ist insbesondere zu beachten, 
laß sich die Gesamtschwingung in eine solche vom ersten Grad, für die a etwas kleiner ist 
als 1, und in eine solche vom zweiten Grad mit a größer als I aufteilt. 

Den Wert aus Gl. (5) setzen wir in Gl. (1) ein und erhalten 


C, | a” 


110 m m,” ro a? 


ebenso verbinden wir die Gl. (2) und (5) und erhalten 
l: k e Sro 


120 sro m, aye,m, 


cz 
- 


Der resultierende Größtausschlag £&,, ist also gleich 


k * u Br 


1 0) \ 5) = m: E D) 
=10 >110 + £&,» ro a? 0°C, m, 





Die im Dämpfer auf eine Schwingung umgesetzte Arbeit ist gleich 7: Größtkraft k- &,, 
mal Größtweg £7,,: 


C, 
ak «| s 
C „, 
A nk&E,.w=nk£&,.u 2 — — ER an ı)), 
ro ’ ro Mm, ( 1? 10 . ( 
\ a?’ a’ ım,C, 


Durch Einsetzen des Wertes £,, aus Gl. (S) haben wir eine Gleichung zwischen a, k und &,, 
erhalten. Wir setzen &,, konstant und denken uns k und infolgedessen auch die erregende 
äußere Kraft P verändert. 

Wenn %k gleich Null ist, wird im Beharrungszustand keine Arbeit im Dämpfer umgesetzt. 
Die äußere Kraft P ist also Null. Wenn k unendlich groß ist, bewegen sich die beiden 
Massen stets gleich. Der Relativweg £, der beiden Massen m, und m, gegeneinander ist 
Null. Es wird ebenfalls keine Arbeit im Dämpfer umgesetzt und P’ ist Null. Wir suchen 
den Wert von k, zu dem der Größtwert A bei konstantem Wert &,, gehört. Durch Differen- 
tiation von Gl]. (9) erhalten wir 


da N ' a?\® R,’ N 
s - (10). 
dk a” a?c,m, | 
Daraus findet man den günstigsten Wert %k, für die Reibung: 
N 1—a’ —: 4 
Be I RE N ET NE N el ©: 
P une & Say, | 


KEın Vergleich von Gl. (11) mit den Gl. (6) und (7) zeigt, daß die günstigste Reibungsziffer %k, 
dann erhalten wird, wenn &,, gleich &s0., d. h. wenn der Vektor &, um 45° gegen den 
Vektor £,, verschoben ist. Dieses Ergebnis stimmt mit dem früher erhaltenen Ergebnis einer 
geometrischen Betrachtung überein’). 

Den Wert von a können wir aus Gl. (1) und (5) oder Gl. (6) entnehmen. Es ist: 


S WET ; ® 
2 ro Sro*Sımw 
tt RE (12). 
rot Sı10 Sro* Sımo 3 x + | 
” , r RR i ’ E % | | 
Für den Fall der günstigsten Reibung ist &,0 = &so v2 &.. Mit x haben wir das mittlere 
| öl 1 | | 

Aufpendelungsverhältnis : - der ungedämpften Schwingung bezeichnet. Das negative Vor- 


>110 
zeichen gilt für die Schwingung vom zweiten Grad, auf die wir die Betrachtung erweitern. 
Die Arbeit A,, die bei der günstigsten Dämpfung im Dämpfer umgesetzt wird, erhalten 
wir aus den Gl. (5), (9). (11) und (12): 
A 


0 E20 T 4 
0 T ( | A")C, ro — r 4 | C, Sro . . . . . . . . . . . . . ( 13). 


°), O0. Föpp]l, Ingenieur-Arehiv 1931, S. 350. 
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Wenn die Masse m, des Dämpfers klein ıst gegen m,, können wir 1 gegen = ver- 


re 
nachlässigen. Aus Gl. (11) wird dann: 
l; I \ c,m, it: 
7 V(x + Dr | ec, Im, 2 ar a Eee 0 


ebenso wird aus Gl. (13) angenähert 
ee 


Die Annäheruneswerte Ila und 13a sind schon in ähnlieher Form in einer früheren Ver- 
öffentlichung mitgeteilt worden’), während die genauen Werte in Gl. (11) und (13) bisher 
nicht bekannt waren. 

Wir haben uns nun mit den Schlußfolgerungen zu befassen, die auf Grund der Gl. (11) 
und (13) für die praktische Berechnung des Dämpfers gezogen werden können. Wir müssen 
berücksichtigen, daß die eine Schwingung, die bei der Anordnung e, m, auftreten kann, durch 
das Aufsetzen des Dämpfers in zwei Schwingungen, nämlich die vom ersten und die vom 
zweiten Grade, zerfällt. Es hilft uns dabei nichts, wenn wir den Ausschlag der einen 
Schwinzungsform (etwa den für die Schwingung vom ersten Grad) entsprechend niedrig 
halten, wenn gleichzeitig die Schwingung vom zweiten Grad große Schwingungsausschläge 
der Masse m, ausführt, also wenig gedämpft ıst. Es ıst uns die Aufgabe gestellt, beide 
Schwingungen gleich stark zu dämpfen oder die Dämpfungsarbeit A, in beiden Fällen für 
eleiche Werte £,, gleich groß zu machen. Die Annäherung, die wir in Gl. (11a) und (13a) 
oereben haben, sind deshalb auch die für die Praxis wichtigen Werte, weil wir ja die 
eünstigste Dämpfung %, nicht für die Schwingung vom ersten Grad anders wählen können, 
als bei der Schwingung vom zweiten Grad. Wir müssen deshalb den Mittelwert einhalten, 
der gleich dem Annäherungswert nach Gl. (11a) ist. | 

Die je Schwingung im Dämpfer umgesetzte Formänderungsarbeit ist nach Gl. (13a) 
verhältnisgleich dem Ausschlag £,, der Masse m, und verhältnisgleich dem Relativweg £,, 
der Massen gegeneinander. Wir wollen bei der Schwingung vom ersten Grad die gleiche 
Dämpfung haben wie bei der Schwingung vom zweiten Grad. Für gleichen Ausschlag £,, 
müssen wir also gleiche Relativwege £&,, bei den beiden Schwingungen haben. Bei der 
Schwingung vom ersten Grad schwingen die Massen m, und m, stets in gleicher Richtung, 
bei der Schwingung vom 2. Grad gegeneinander. Im einen Fall subtrahieren, im anderen 
Fall addieren sich die Wege &, und £&, der beiden Massen zum Relativweg &,. Wenn wir 
mit s, das Aufpendelungsverhältnis 

>110 
mit s, den gleichen Wert für die Schwingung vom zweiten Grad bezeichnen, dann muß für 
aleiche Dämpfung s, — 1 gleich sein s,—+1. Wir müssen die Feder ce, oder die Masse m, 
des Dämpfers so berechnen, daß diese Bedingung erfüllt wird. Die Unterlagen für diese 
Berechnung sind im Ingenieur-Archiv 1931, Seite 348, angegeben. 


der ungedämpften Schwingung vom ersten Grad und 


Endlich ist noch darauf hinzuweisen, daß in den Fällen der Praxis m, klein ist gegen m.. 
Infolgedessen ist &,, groß gegen &,. Wir können angenähert £,, durch £&,, ersetzen. Auf 
alle Fälle ist der Mittelwert von £&,, für die Schwingungen vom ersten und zweiten Grad bei 
richtiger Bemessung des Dämpfers im obigen Sinne gleich £,,, da ja der Dämpfer so abge- 
stimmt sein muß, daß &,, für die Schwingung vom zweiten Grad ebenso groß ist wie &,, für 
die Schwingung vom ersten Grad. 

Mit den vorstehenden Ausführungen ist es möglich, einerseits den Dämpfer richtig zu 
dimensionieren, d. h. die Federstärke ec, zu einem gegebenen m, und zur gegebenen An- 
ordnung e, m, zu bestimmen und andererseits die günstigste Reibungsziffer k, zu ermitteln. 
In vielen Fällen der Praxis liegt allerdings nicht die einfache Anordnung der Abb. 2 vor, 
bei der die Feder e, an einem Ende festgehalten ist und am anderen Ende die Masse m, 
trägt. Bei den Drehsehwingungen bei Kurbelwellen ist es z. B. so, daß an beiden Enden 
der Welle Massen sitzen und daß ein Knotenpunkt in der Mitte der Welle vorhanden ist. 
Wenn die Masse m, des aufgesetzten Dämpfers entsprechend klein ist gegen die Schwung- 
massen auf der Kurbelwelle, dann kann man den Knotenpunkt als festliegend annehmen und 
erhält eine ähnliehe Anordnung wie in Abb. 2. Aus der praktischen Messung, wie weit die 
Kigenschwingungszahlen vom ersten und zweiten Grad voneinander liegen, kann man einen 
kückschluß auf die reduzierte Masse m, der Kurbelwelle erhalten. 251 


3,0, Föppl: Sehweizerische Bauzeitung 1931, Band 97. 
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Schwingungsdämpfung 
unter Ausnutzung der Werkstoffdämpfung. 
Von @ünther Bock in Berlin. 


rehschwingungen in Kurbelwellen können durch Anbringen von Schwingungsdämpfern 

bekämpft werden. Bekannt ist der Junkers-Dämpfer, der die Schwingungsenergie 
dureh Flüssigkeitsreibung aufzehrt, andere Dämpfer arbeiten mit Oberflächenreibung. O. Föppl') 
nutzt die innere Reibung in festen Baustoffen zum gleichen Zwecke aus. Versuche mit solehen 
Dämpfern werden a. a. O.?) ausführlich beschrieben. Hier wird die Theorie dieser Dämpfer 
entwiekelt und die Berechnung ihrer Abmessungen angegeben. Die Entwieklune führt auf 
erzwungene Schwingungen zweier gekoppelter Systeme, von denen das eine aus einer mit 
Masse behafteten Feder (eben dem Dämpfer) besteht?). Die Systeme sind durch Kraftkopplung 
und Reibungskopplung verbunden. 


Dem für den Schwingungsdämpfer ausgewählten Baustoff kann man verschiedene Formen 
reben: Stabform einseitig eingespannt mit freischwingendem Ende oder beiderseits eingespannt, 
Vollstab oder Hohlstab, Scheibenform oder Scheibenringform, innen oder außen eingespannt. 
Hinsichtlich der Baustoffanstrengung und der Wärmeabfuhr ist die ungünstigste Form die 
innen eingespannte Scheibe. Bei der außen eingespannten Scheibe liegen die Verhältnisse 
wesentlich günstiger. Auch die Stabform ist nicht ungünstig. Da sich rechnerisch deren Ver- 
halten am bequemsten verfolgen läßt, beschränken wir die Untersuchung auf stabförmige 
Dämpfer. Die gewonnenen Erkenntnisse gelten in ähnlicher Weise auch für die anderen 
Dämpferformen. 

















Abh. 1. Abb. 2. 


1. Die Schwingungsgleichung des Dämpfers. Der Dämpfer habe die Form eines zylın- 
drischen Stabes, der an einem Ende eingespannt ist (Abb. 1). Es soll zuerst untersucht werden, 
wie sich der Stab verhält, wenn der Einspannstelle Drehsehwingungen um die Stabachse von 
bestimmten Ausschlägen und bestimmter Schnelle aufgezwungen werden. Für die praktische 
Verwendung eignet sich besser der beiderseits eingespannte Stab. Er schwingt in derselben 
Weise wie ein einseitig eingespannter Stab von halber Länge und braucht daher nicht be- 
sonders untersucht zu werden. 


Den Beziehungen zwischen den Spannungen r und den Formänderungen y dürfen wir 
nicht das Hooke'sche Gesetz 7= Gy mit konstantem Gleitmaß @ zugrunde legen, denn dessen 
strenge Gültigkeit bedingt, daß bei zu- und abnehmender Verdrehung im Stabe die gleiche 
linergiemenge gebunden und wieder freigegeben wird. Wir wollen aber gerade diejenige Energie- 
menge — wir kennzeichnen sie kurz als „Dämpfung“ — untersuchen, die der Stab nicht wieder 
als Bewegungsenergie herausgibt, sondern in seinem Innern in Wärme umwandelt. Deshalb 
gilt auch für ein Kräftepaar, das am freien Ende des ruhenden Stabes von Halbmesser r, mit 
dem Moment M angreift und im Stabe die Randspannung r, hervorruft, nicht mehr die Be- 


2M ’ i u 
3. Die wahre Randspannung bleibt vielmehr hinter diesem Werte etwas 
€ 0 


zurück. Wir wollen jedoch diesen Ausdruck als „rechnerische Randspannung r’, 


z\ehung Be — 


‚“ beibehalten: 
1,0. Föppl: Schwingungsdämpfer für Kurbelwellen. Ing.-Archiv 1950, 8. 222. 
2) G. Boek: Versuche mit Schwingungsdämpfern. Der Aufsatz erscheint demnächst. 


3) Über Dämpfer mit getrennter Masse und Feder vgl. L. Gümbel: Verdrehungsschwingungen und ihre 


Dämpfung. ZVDI Bd. 66, S. 252 bis 256 und 281 bis 282 (1922). J.Ormondroyd and J. P den Hartog, The theory 
of the dynamie vibration absorber. Transactions of the American soe. of mech. eng., appl. mech. Bd. 50, 8.9 bis 15 
27). — J. P. den Hartog, foreed vibrations with eombined Coulomb and viscous frietion. Trans. Amer, soc. 


ech. eng. APM Ba. 53 (1931), S. 107 bis 115. 
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‚_.2M er : ss h ER 

r', ‚„ oder, auch für Hohlstäbe gültig, mit dem polaren Trägheitsmoment J des 
nr, 


Stabquerschnittes 


Vr 


’ e 0 


= 7 


Aus statischen Verdrehungsversuchen läßt sich für eine bestimmte Querschnittsform r’, in 
Abhängigkeit von der Randformänderung y, bequem bestimmen. Die y,r',-Kurve ergibt dann 
in einfacher Weise die wahre Randspannung r,') und die ein- und ausgeleiteten Energiemengen. 


Wir ändern jetzt die Größe des am Stabende angreifenden Momentes M und führen 

lanesam einen vollen Lastwechsel von + M zu M und wieder zu +M aus. Dabei ändert 
. 

I 
Messen wir stufenweise M und g und zeichnen die y,r',-Kurve, so erhalten wir einen ge- 
schlossenen Linienzug, die „Hysteresisschleife“ (Abb. 2). Die Flächen unter den Ästen der 
Hysteresisschleife entsprechen den von M abgegebenen und aufgenommenen Arbeitsbeträgen A, 
bezoren auf die Volumeinheit: 


0 


sich der Ausschlagwinkel 9 des Stabendes und damit auch die Randformänderung y,=4 


A “=; I(J , I E, 
" V V \ Mr j da fl \ N, ü ° N, " . \ P ’ e To 
"; . 
I I r Y 
Der Formfaktor o fr> >» berücksichtigt die Hohlstabformen. Für den vollen Kreis- 


querschnitt ist o == '/s. Die von der Hysteresisschleife eingeschlossene Fläche f, entspricht 
der für eine volle Schwingung aufgewandten und in Wärme umgewandelten Arbeit, der 


Dämpfung 9. Es ist Je ofh. 


Die Lastwechsel sollen nun mit bestimmter Schnelle » aufeinander folgen, wobei sich 
das Moment M sinusförmig ändern soll. Der Stab gerät dabei in Schwingungen. Von den 
abklingenden Eigenschwingungen sehen wir ab und nehmen für die erzwungenen Schwingungen 
zur Vereinfachung der Rechnung sinusförmigen Verlauf an. Wir ersetzen deshalb die Hv- 
steresisschleife durch eine Ellipse’) von gleichem Inhalt, die durch die Spitze der Hysteresis- 
schleife verläuft und hier den gleichen y,-Größtwert erreicht. Die Ellipse ist in Abb. 2 
gestrichelt eingetragen. y,,, T,, Sollen das zur Schleifenspitze gehörige Wertepaar bezeichnen, 
t’,,„ den Größtwert der Ellipse. Die Beziehung zwischen den Augenblickswerten y, und T’, 


zur Zeit ? können wir darstellen durch das Gleiehungspaar 


‚sinoft 7, r.sin(ot+a). 


nr nr 
»v Pa 


r r ® T ® ’ 
Zur Zeit wi=— ist ,=Yar; Fr ta tab. 


Wir führen 7’,, in die zweite Gleichung ein: 


’ 


zu ro sin (ot-+a) = Tr, (sınot+tgacosmfß). 


" cos a 


a muß so bestimmt werden, daß die Ellipse der Hysteresisfläche inhaltsgleich wird. Der 
Kllipseninhalt ist f, 7 7,,7.,f9@. Wenn wir die Dämpfung d = ofn in Beziehung setzen 


‚v1 
| | 
zu der im verdrehten Stabe aufgespeicherten Formänderungsarbeit « = 5 07,17, ),so erhalten 
r M au j » ’ 0 fin . 5 A 
wir die „verhältnismäßige Dämpfung y*: y en: ‚und damit die Hysteresis- 
( ’ 
-0Y T 


so) Yo o1 


[4 


| 5 
Häche fr 5 Y Yo To. y kann leicht dureh Versuche bestimmt werden. Aus fu = fer yyıT., fga@ 


ı, Vel. Föppl, Beeker, v. Heydekampf: Die Dauerprüfung der Werkstoffe. 8.12. Berlin 1929, Springer, 
oder Knacekstedt. Die Werkstofflämpfung bei Drehsehwingungen nach dem Dauerprüfverfahren und dem Aus 
schwingeverfahren. Berlin 1930, NEM-Verlag. 

5), Vgl. aueh z. B. Busemann, Diss. Braunschweig 1925. 

6, Dieser Ausdrnek entsprieht dem Flächeninhalt des unter den Verbindungsgeraden OA (Abb. 2) gelegenen 
Dreiecks. Der genaue Wert der Formänderungsarbeit a, der aus Ausschwingversuchen bestimmt werden kann, ist 


etwas größer als der hier benötigte Ausdruck. Folgerichtig müßte deshalb hier a’ und Y” geschrieben werden. 


\l« 


Di 


d 1 


lıe 


die 
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! 

/ . . \ . Ir . . . 1 .g 
letfga=, 5 Damit lautet die Gleichung für die Randspannung, wenn wir noch die Größe 
'— " einführen, die bei Gültigkeit des Hookeschen Gesetzes mit dem Gleitmodul @ 


N — 


Yoı 
lentisch wäre, r’ G'y [ein ot—+ 4 cosot)oder, da y y.‚sinwf dya y os of 
\entisech ale, T, Yoıl! ) 97, ) ‚aay, Yoı! 07) Zr Yoı cost, 
u» dy 
’ ‚, / je» 
T (i B = 2 . . . : P . . ; ö : ; ; ; i : ; | . 
0 797 rm dt (1) 


Wir betrachten jetzt, am Dämpfer Abb. 1 in beliebigem Abstande = von der Einspannung 
ine dünne Scheibe von der Stärke dx. An der Scheibe greift auf der einen Seite das 


\Ioment M, an, auf der andern das Moment 
oM 
+ M, rdx M.+ er dr. 
Ks ıst 
J J@G' vr OYyır 
an ' Y / 0X 
M,. — y T 0x Y ’UX + Ir m Ri N 


4 


; bezeichne den augenblicklichen Ausschlag an der Stelle x. Da g hier nicht von der Eın- 
spannung bis zum freien Ende gleichmäßig zunimmt wie bei der statischen Verdrehung, ist 
2 ua 0g 

jetzt für die Randformänderung zu setzen y, = r,, ‚. Damit wird 
i 0x 


Ri 


1 074g 
M.—J@ / N PSCR SR: 
pp 1 [2 I N hi «ol . ( 


d.ıM I@ e: y 0’ y | 
d.r ” Nur? It 0arot 


Die Bedingung, daß die beiden Momente M, und M,; ua, der Trägheitswirkung der Scheibe 
0°g - R 

dm \f: dlas Gleichgewicht halten müssen: 
AWE - 


0’ 
M.+M,+ax dm, j. 0, 


liefert, wenn u die bezogene Masse des Dämpfers ist, mit 
dm=uJdx und M.+4x=NM:t 3}. dı 
: a, 


die partielle Differentialgleichung der Dämpferschwingung: 


NM Ne 
Ö. 007g 
Ir Jdx = 0, 
Nr de uJda Ar 
IN? Pr \3 N2 
„I Pr, v 0g 049 
J( Bas mi — Ix F = Ir D, 
: O8? Irwoarot ai A OÖ pe | 
‚dh BO MER... OR 3.5. OR (3) 
Ox’  2rwmoxr’ot or ne 2 I  3ar 3 


Streng genommen müßte in Gl. (3) berücksichtigt werden, daß sich y und @’ in gewissen 
Grenzen mit y,, also mit g,., ändert. Bei den bisher ausgeführten Versuchen wurden in 
erster Linie Dämpfer aus Gummi verwandt. Bei Gummi ist die Änderung verhältnismäßig 
gering, daher dürfen für y und @’ konstante mittlere Werte eingesetzt werden. Wir kommen 
hierauf später zurück. Für andere Baustoffe muß diese Frage von Fall zu Fall entschieden 
werden. 


ur 
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Die Lösungen einer Differentialgleichung von der Form der Gl. (3) sind bekannt’). Wir 
hatten eingangs für den Dämpfer vorgeschrieben, die Einspannstelle solle eine bestimmte 
Schwingung ausführen. Wir nennen die Schwingungsweite der Einspannstelle y, und setzen 
für den augenbliekliehen Ausschlag 9x —- 0 = g,Ssin®#. Damit haben wir 2 Grenzbedingungen 
für die Lösung festgelegt. 2 weitere folgen daraus, daß das freie Ende stets spannungslos 
/ 


L 


\ 
i; no, 0% Eriz 
ist: für #1 muß in jedem Augenblick N 0 gelten. Mit diesen Grenzbedingungen ergibt 


sich für die erzwuneene Schwingung die Lösung 








( - i 
y 70 (asın af + peosmf) 
) 
a pP q x(Asınop.« + Becosop.r) + e v4(Asnop. ! UCOSsmpx) 
p ev 4x (Bsinop.# Acosopa)— e u sınop.«—+ Acoso px) 
n sin 2opl l MR | u ne (4). 
( een en 
I m 2 n 
B=e-?v«al] cos2mpI 
| mw —+ | u 
P .) ‚ 
) [ZB je 
"=et?"ral+ cos2wpl ’ 
! 
/ 
ei a m l- 
I) e' zoqlLe -oql + 2cos?2wpl Pie- 





Die Größen p und 4 (Dimension sece/cem) und die dimensionslose Größe w wurden zur Ab- 
kürzung eingeführt. Man kann sich leicht überzeugen, daß Gl. (4) sowohl die Differential- 
eleichung als auch die Grenzbedingungen erfüllt. 


. > u! ’ . . ui » . . . \ 
Setzen wir voraus, daß 2 klein bleibt gegen 1, so dürfen wir die vereinfachten Be- 


ziehungen einführen 


r 
/ 


’ ’ 


/ Mu ‘4 v er 
Ve an. 


u Y 
v | (1' 1 T pP 


T tı7 


Es ergibt sich z. B. mit y' = 0,6, dem höchsten bei den untersuchten Gummisorten gemessenen 
Dämpfungswert, 


y' u 


u — ] 1 +0, W919 1.004549 p 0,9966 e- EL 


4= 0,9943 


Y 


. . \ a . y 
Ebenso werden im foleenden Faktoren von der Form |1 + c(, = 


fortzelassen, sofern sie nicht in Ausdrücken auftreten, deren Glieder verschiedene Vorzeichen 
annehmen können. 





. wobei e< 1, stillschweigend 


2. Die Schwingungsform. Wenn wir 0 setzen, erhalten wir die ungedämpfte 
Schwingung. Es wird p =0 und 


7. 9,wosopa+ttgoplsnopa)snot . 2 2 2 222... (B). 
Vs .T . Pr » . “ > n . 
Für opl  , wird go: der Dämpfer schwingt in Resonanz zur Befestigungsstelle. 


Diese Eigenschnelle des Dämpfers bezeichnen wir mit oa, 


7t I (1 ü 
02) a ee he 5 
ap NV (9 
> > 
Die Eigenschwingungen höheren Grades &op/ 7, 57T... können hier außer acht bleiben. 


») 


’, Vel. z. B. Holzer: Die Bereehnung von Drehsehwingungen und ihre Anwendung im Maschinenbau, S. 188. 
Berlin 1921, Springer. 
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Vir können jetzt die beliebige Schwingungsschnelle » durch ihr Verhältnis zur Eigenschnelle 
es Dämpfers kennzeichnen, durch das „Verstimmungsverhältnis 2“: 


z7ı Zn &E 
() p I zZ „DA p / \ (dt) yp N e ] 
“) ZVDdA 
4 V ZU x zyv X 
DAGE DPE , 
a Fr 37.53 


Diese Ausdrücke führen wir in die Hauptgleichung 4 ein, gleichzeitig vereinfachen wir die 
Gleichung noch, indem wir die Exponentialfunktionen in Reihen verwandeln, deren Glieder 
‚öherer als 2. Ordnung wir vernachlässigen dürfen. Die Hauptgleiehung nimmt jetzt die 
(sestalt an 

















( P 2 
yo= j; (asınat-+Pcoswt) 
| P yaı\|. ra „.. Zn%E ( a (? y\” „an 
a +17 sin -2sinzr + cos + 608? 
Falsı 21 KT Ve i 
(“) 
’ zy . ZIC x „on ERE >» 
p 9 |sın "YA | N cos’ —5 |+ cos 51 > sn zn 
zıı zy\" 
na? ni 
D == 4c0s 5 714 
Wir betrachten zuerst die Schwingungsform im Resonanzfall. Für z = 1 wird cos „ 0, 
| zZ y' Av ’ 
ferner, da 387] 1, 
” l N ITX . y) ’ IT& 
7 COS ? sin 
16 1 21 2 997 
und 
(! — x NE , 8 ., NE 

ne COS sin @ — sin cos ot). 
Fre 7 21 vy 2] 


Das zweite Glied überwiegt das erste weit und bestimmt den Charakter der Schwingungsform: 
die Sinuskurve mit größter Steilheit an der Einspannstelle. Abb. 3 zeigt die einzelnen 
Phasen der Schwingung. Die Kurve der Größtausschläge ist gestrichelt eingetragen. Das 
IT 

= nach. 


freie Ende eilt der Befestigungsstelle um 


Dehnen wir unsere Betrachtung auf andere Schwingungsschnellen aus, z2<_ 1 und z >1, 
so erhalten wir Schwingungsbilder wie Abb. 4. Dort sind für den Zeitpunkt, in dem der 
Ausschlag des freien Stabendes seinen Größtwert erreicht, die Schwingungsformen bei ver- 
schiedenen Erregungsschnellen zusammengestellt. Wir sehen, der Sinuscharakter bleibt 
erhalten, doch verschiebt sich der Anfangspunkt der Sinuslinie. Für die ungedämpfte 
Schwingung y—=0 würden wir reine Sinusformen erhalten, deren Knotenpunkte auf dem 
/ylinder liegen oder außerhalb desselben, je nachdem z > 1 oder 2< 1. 


3. Resonanzkurve und Phasenverschiebung. Wir fassen jetzt lediglich das freie Ende des 
Dämpfers ins Auge. Mit #1 wird 
/ i zı\ z7ı 2 z u’ a Z7ı Po 2 
er (i + o | 3 )cos ) Pi n Ssın 5 / D ne at Pi URL: ; er 


Wir finden daraus die Schwingungsweite des freien Endes 


Po z 55 Yo j „an, [sw\ 
pı = D } A] +Pı „en, (ey 2 | 16 [cos? 5 +| 2 
Te 
* %)) 
VI en zp\? a ee 
\ cos’, + >) 
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' | en 
Im Resonanzfall 2° I wird cos „0, 
27 8 
E aa 
q 0 y 
2. B. mit y=05 97° 16g,. Das bedeutet: durch die Schwingung der Befestigungsstelle 


wird das freie Ende zu einer Resonanzschwingung von 16fachem Ausschlage angeregt. Diese 
eroße Formänderung ist es, die dem schwingenden Stab die Möglichkeit gibt, selbst bei ver- 


verhältnismäßig kleinem Ausschlage der Befestigungsstelle, einen nenneswerten Arbeitsbetrag 
aufzunehmen und in Wärme umzuwandeln und damit als Dämpfer zu wirken. Sie ist auch 
der Grund, weshalb als Dämpferbaustoff besonders Gummi in Frage kommt. 

















16 
P 
Yo 
70 
0° u u 
Wella 0) 
AN bh. w 














Abb. 3 und 4. Abb. 6. 


Nach Gl. (9a) hängt im Resonanzfall das Ausschlagverhältnis nur von y ab. Die ver- 
| - / 
hältnismäßige Dämpfung y ist demnach leicht durch einen Schwingungsversuch zu bestimmen, 
wenn g, und 97 gemessen werden. Vgl. Fußnote 2. 


Abb. 5 zeigt das Ausschlagverhältnis 97:9, °) in seiner Abhängigkeit von der Schwingungs- 
schnelle der Befestigung. Der Kurve liegt der Wert y =0, zugrunde, der etwa elastischem 
Gummi entspricht. Die scharf ausgeprägte Resonanzkurve legt die Frage nahe, ob nicht die 
beabsichtigte Wirkung des Dämpfers, die ja mit dem Ausschlagverhältnis in enger Beziehung 
steht, bei der Drehzahländerung zu rasch abklingt. Wir sehen, daß zur Beurteilung des 
Dämpfers der ganze Drehzahlbereich in Betracht gezogen werden muß. 


Wie Abb. 3 zeigt, eilt die Einspannstelle in ihrer Schwingung dem freien Ende voraus. 
Wir finden die Phasenverschiebung &, indem wir für z=0 und I! dg/dt gleich Null setzen: 


7 dj S zı 
ol, tot, tg”, 
u pP1 “ a 
I zı7ı JI zı 2 IT 
ef), of arte | -cig 5 .) E—- arte | S tg r 


Abb. 6 zeigt. wie e sich periodisch sprunghaft ändert. Innerhalb einer Periode ähnelt 
das Verhalten dem einer schwingenden Einzelmasse. 


s, In dem eingangs erwähnten Aufsatz von O.Föppl (Fußnote 2) wird das Verhältnis 7:90 als „Aufschaukelungs 


erhältnis“* bezeichnet. Wir wollen jedoeh diesen Ausdruck hier vermeiden und ihn für einen später auftauchenden 
Begriff aufsparen, der sieh auf die Gesamtanlage bezieht. 
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4. Dämpferkraft und vernichtete Arbeit. Für die Berechnung der Dämpferwirkung müssen 
vir die Dämpferkraft kennen, d. i. das Moment, mit dem der Dämpfer an der Befestigungs- 
telle auf deren Schwingbewegung antwortet. Wir erhalten, zunächst für eine beliebige 


ı) ‘telle x, aus der Hauptgleichung (4) 
\ Y >B y? . 48 \B 
a Ah, Ik m. op | 0g Po (‘ OB . 
sınol- cos ol © cos 1 q 
le dr D\O 1 v rn 7 l 237 D l er 7 pr sın of 
St 
l- nd setzen diese Ausdrücke in Gl. (2) ein. Es wird, für jeden Wert x gültie, 
| 
h a yv Op Op v Oa\| 
h ] Ja .. ( / l bzo e / E / | 
M T D sinotl, e- ea z. cos ot, = 2 = j 





Kür x==0 fassen wir die Klammerausdrücke mit D und I! zusammen und führen die Feder- 
konstante des Dämpfers ce, ein: 


/ R a y A pP l > p My N?) JG 
u - " C, 
D\O.x 0X =D D\o# PLAN) EX N I 
M.-„=sinwt-p,cau—coswti-p Cd! .: .: 2. 2 2,2020... (I0). 


Die beiden Glieder können wir als Projektionen zweier um — phasenverschobener Moment- 
vektoren My’ =g,cau und Ma" Y,cav auffassen. My schwingt gleichphasig mit 9,, es 


leistet keine Arbeit. My’ eilt dem Ausschlag 9, um nach und nımmt aus der Anschluß- 


> 
stelle die gesamte Dämpfungsarbeit auf, die im Dämpfer aufgezehrt wird. 


« und » können wir deuten als Verhältnis der durch die Schwingung hervorgerufenen 
Momente My und Ma’ zu dem bei statischer Verdrehung des Dämpfers um den Winkel y, 
auftretenden Moment e4Y,. Um uw und » zu bestimmen, differentiieren wir a und 5 und 
finden für 2 —0 


da | ( r zıp\® dp zı 
zıasınza | (zr sınzrı 
u. 7 | ey l 
zasınza ®) zy za -—-sSinzn 
u = - [Bi .) “ no. ö . ( | | ). 
„en (Zy\ “2 „er, (zwy\ 
4 cos > -( COS > -( 4 


Bei jeder Schwingung vernichtet der Dämpfer die Arbeit 


A Toy, My’ Tg n a A rer oo . . (12), 
im Resonanzfall 
In? 2 
A = a Cd q nl . 
/ 


Dieselben Werte für die Arbeit erhalten wir, wenn wir von der Dämpfungsarbeit einer 
Scheibe d.r ausgehen und über die Länge ! integrieren. 'Tun wir dies in 2 Stufen: #—=0...1/2 
und #=1/2...1, so finden wir, daß von der gesamten Arbeit die der Anschlußstelle benach- 
barte Hälfte des Stabes 82°/,, die äußere Hälfte nur 18°/, vernichtet. Dementsprechend 
erwärmt sich auch, wie die Versuche bestätigen, die angeschlossene Stabhälfte wesentlich 
mehr als die freie. 


5. Der verkürzte Dämpfer. Die soeben ausgesprochene Tatsache, daß die äußere Stab- 
hälfte nur 18°/, zur Arbeitsvernichtung beiträgt, legt den Gedanken nahe, diese Hälfte ganz 
oder teilweise wegzulassen und durch eine andere Masse so zu ersetzen, daß die Schwingungs- 
schnelle erhalten bleibt. Namentlich bei solchen Wellen, die verhältnismäßig langsam 
schwingen und demgemäß eine größere Dämpferlänge erfordern, kann eine Verringerung der 
Baulänge auf die Hälfte oder weniger erwünscht sein. In jedem Falle wird man eine be- 


»), Bei uw und » ist im Zähler ein Glied mit dem Faktor 2/16 vernachlässigt. Streng genommen wird u nieht bei 
=] gleich Null, sondern z. B. für 7 —=0,6 bei z= 1,23. 


BT 
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stimmte Masse am Stabende anbringen, schon um durch ihre Vergrößerung oder Verkleinerun: 
den Dämpfer bequem auf die Eigenschwingung der Welle abstimmen zu können. Wir frageı 
jetzt: welchen Einfluß hat eine solche Zusatzmasse auf den Ausschlag und die Wirkung des 
Dämpfers, und wie groß muß die Masse für eine bestimmte Verkürzung der Dämpferläng: 
gewählt werden ? 














Fr L ” 
G ma  |\Fa% | 
u 
RA24527 Io u 
Abb. 7. 


Wir dehnen unsere bisherige Untersuchung auf den Stab von der Länge I! aus, dessen 
einem Ende wie bisher die Schwingung mit konstantem Ausschlage y, aufgezwungen wird, 
und dessen anderes Ende jetzt die Masse m, trägt (Abb. 7). Zum Vergleich ziehen wir den 
gleich abgestimmten Stab mit freiem Ende und der Länge I, heran. Die Differentialgleichung 
Gl. (3) behält ihre Gültigkeit. Für «1 ist jetzt an Stelle der Spannungsfreiheit um 

w ( oe 
das Gleichgewicht zwischen dem Stabmoment M und der Trägheitskraft der Masse m. Be- 
dinzung: 


0” 


) Ve ( 
7) + Me 5 Be ne er 


v‚ Q 4 v 
J (1 ._ ) 
Q x zZ’ dr dt 


Damit ergeben sıch für die Hauptgleiehung (4), die im übrigen ihre Form behält, andere 
konstanten: 











| .. ® > N ) | ı t l ) ] ’ m. 
(1 ı)sınZzaoapt tr ZI COSZEZDO»P = 
/ > | Jyu@' 
l ); m I | . > .. | ) »9 iy 2 
b Di 5 sn2opl+ (w —t)cos2wpl+ e-2"al\w-+i’—, | 
u PL mv +] 
| I/w+1. | . |  , ty 2 
( 2 » Ssn2opl+ (m — it) cos!2!opl-+ e! warlw+r +. er dia). 
2 | In w+ 1, 
. m T | . ) ) .. TR } .) f I) 
I bi » sin2oqgl+2(w — i?) cos2mpT-+ (et?"al + e-?®al) (wi?) 
y .) 
ıy 
-(e+?vgl — 2ogl / >> 
+ (e 4 fi a) 2 3 
a: 1 





Die früheren Konstanten gehen hieraus hervor, wenn wir m, und damit i gleich Null 
setzen und ı wegheben. Wir betrachten zunächst wieder die ungedämpfte Schwingung und 
erhalten aus Gl. (4) und (14) 


2 ee i”) sin 2opil + »iCcos 2opl 


O=Q, .— — siınmpx + cosop.x|sın wt. 
d ! 4(cosopI ı sın pl) | r 1 
Die Resonanz 4 x tritt ein bei @© = i,—etgopl. Die zugehörige Schwingungsschnelle, die 


Eigensehnelle der ungedämpften Schwingung, nennen wir wieder &q. Wir vergleichen diesen 
Dämpfer mit dem zweiten, der keine Zusatzmasse trägt, aber vermöge seiner größeren Länge I, 
mit derselben Eigenschnelle 4 schwingt, und führen das „Kürzungsverhältnis 4* = I/I, ein. 
Für den ungekürzten Dämpfer fanden wir o&apl, z/2, für den im Verhältnis 4 gekürzten 


AT ATI 


gilt entsprechend map! wapıl, >. Damit wird i,=ctg „. Anderseits ergibt sich aus 
ir ij 7 @' _: 
(4l. (14). wenn wir in den Ausdruck für i für den Resonanzfall & = wa >] einsetzen 
u) au 
und /, durch 7 ausdrücken, 
m. ia ör E$ in m. 
i, ;* 
 JyuG zl u 2 Jul 2 ma 
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obei ma = uaJl das Massenträgheitsmoment des gekürzten Dämpfers ist. Die beiden 


‚jeiehungen für i, vereinigen wir: 
7 


/ 
cig — 


2 ma 2 


La m; 


\amit haben wir jetzt die Zusatzmassen m, bestimmt, die an dem im Verhältnis 4 gekürzten 
impfer angebracht werden muß, um ihn wieder auf &a abzustimmen: 


2 in n 
Mm. Mg, ctg 5 (15). 
u ’ ; 22) 
\ußerhalb der Resonanz wird mit 2 
j Da 
n. f in h 
n ir» — 2cte 6 
DA ’ I ( ) 


\bb. S zeigt das Anwachsen der Zusatzmasse mit der Kürzung. Beı geringer Kürzung ent- 
sprieht das Trägheitsmoment der Zusatzmasse angenähert dem des abgeschnittenen Dämpfer- 
pricht d g 8 


St iickes. 
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u 
My Pi 2X pen. © PR 
1. #T_— * 
je" > 
A 7 A 1 RA23SZB..0) A 7 
Abb. 8, 9 und 10. 
5 N in 2 
Bei sehr starker Kürzung darf für etg gesetzt werden: 
1 4 > 
A: Mm, = 75 Mi. Da). 
2 En? d ( L) 


Die Masse des Dämpfers tritt jetzt zurück gegen die Zusatzmasse. Wenn wir in Gl. (6) die 
4aJl und Gl. (15a) einführen, so finden wir 


% rJ W” n° 
l Mad 4 Er 


Größe ma 





@'.J 
Im; 


(G'J1 
Mg 


/ vr’ 
JT / T 7T 


M) / : .) 
21, u 21, 


und haben damit den Anschluß an das Gebiet der Einzelmassenschwingung gewonnen. 
Um die Auswirkungen des Dämpfers zu bestimmen, brauchen wir wieder die Momente 


u) = 














My g,cau und My" f,€carv und berechnen « und v auf dieselbe Weise wie früher: 
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Wir dürfen wieder »w durch 1 ersetzen und bei « im Zähler das dritte Glied vernachlässigen. 
Abb. 9 und 10 zeigen den Einfluß der Kürzung auf Ausschlagverhältnis und Dämpfer- 
kraft für den Fall ® = w,„, also z=1. Gl. (17) liefert hierfür 
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Das Ausschlagverhältnis berechnen wir für diesen Fall ähnlich wie unter Ziffer 3: 


AR 

r % sIN 

ed A e 
7 Po y tr +-sın/ırn 


Während der Ausschlag g; mit zunehmender Kürzung langsam abnimmt, wird die Dämpfe: 
kraft durch starke Kürzung wesentlich gesteigert. 
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ze RR dq u i 
Die größte Randformänderung y,=r, m der Anschlußstelle ist im Falle 2=1 
N 
( Y 
a 
| 
Y a . 0 N, l T . . v*. » . . . 
Für A=1 wird y, I wie sich schon aus der Sinusform der Schwingung mit dem 
1) 
Endausschlag 9; Po ergibt. 
/ 


6. Die Resonanzverschiebung des schwingenden Systems. Wir haben jetzt die Eigen- 
schaften des Dämpfers kennengelernt und wollen nun wissen, wie sie sich auf die schwingende 
Welle, an die der Dämpfer angeschlossen wird, auswirken. Wir untersuchen das einfachste 
schwingende System, eine Welle, die an einem Ende eingespannt ist und am andern eine 
Schwungmasse trägt, und vergleichen ihr Verhalten ohne und mit Dämpfer, Abb. Ila und b'"). K, 
Die übrigen Schwingungssysteme lassen sich in erster Annäherung auf ein solches Ein-Massen- 
System zurückführen, das in Eigenschnelle, Ausschlag und eingeleiteter Schwingungsenergie 
mit dem wirklichen System übereinstimmt. Man kann dann am Ein-Massen-System die 
Dämpferwirkung leicht überblicken und nötigenfalls eine genaue Durchrechnung des Mehr- 


hi 
Massen-Svstems nachholen. 
. y \ „Fu Jw . . . S R, 
Die Welle Abb. 1la, deren Federkraft ce ist, trägt die Schwungnfs#sse von der 
u 
Trägheit m. Auf die Masse soll ein mit der Schnelle & pulsierendes „erregendes“ Drehmoment 
einwirken, das Masse und Feder in Schwingung versetzt. Das rückdrehende Moment der d 
Welle ist - ey, wenn g den augenblicklichen Verdrehungswinkel des Wellenendes bezeichnet. d 
Wir nehmen an, daß die Masse der Welle gegen die Schwungmasse sehr klein ist und ver- v 
(; 
Dr h RR dag ’ 
nachlässiet werden darf. Auf die Masse wirke noch das dämpfende Moment — k TE Die 
( 
er : > ar r da 
Gleichgewiehtsbedingung der an m angreifenden Momente lautet damit M—cey- k lt 
d’gq .- h » 
mp» Wir betrachten eine erzwungene Schwingung und setzen 9=g,sinwt. Das 
( x . 


erregende Moment können wir in zwei Komponenten M, und M,.” zerlegen, von denen M,' 
7 
in gleicher Phase mit y, M,” um , dagegen verschoben schwingt. Der Augenblickswert des 


Momentes ist dann M = My sinot-+ M.”’cosot, der Größtwert M.=yM,.”—+ M,.””. Die 
KEigenschnelle der ungedämpften Schwingung — nennen wir @%. Die gedämpfte Schwingung 
hat bei o»,. nahezu ihren Größtausschlag. 

Jetzt setzen wir auf die Schwungmasse den Dämpfer auf, der auf die Eigenschnelle der 


Welle abgestimmt sein soll, Abb. I1b. Er übt auf m das Moment M,„, aus, dessen Größe wir 
aus Gl. (10) kennen. Die Gleichgewichtsbedingung der Momente lautet jetzt 


da dl? 4 
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Mit 4 y,sinof, dem oben entwickelten Ausdruck für M und mit Gl. (10) folgt 


MV, sin nf r M,"cos wtf eq „sin ef kog 0 cos mt + mw" q „sin ft 
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4 
vo 
ww, Foppl zeigt die Resonanzverschiebung dureh den Dämpfer an einer Welle mit gleichmäßig verteilter Masse: 
O0. Föppl: Grundzüge der technischen Schwingungslehre. 2. Aufl. 8. 110. Berlin 193), J. Springer. 
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\bb. 12 zeigt das Zusammenwirken der Momente im Vektordiagramm. 
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Abb. Ila und 11b. Abb. 13. 


Für die ungedämpfte Schwingung setzen wir k=-0 und „0 und damit auch » 0, 
Ks wird 
M, 


( z 
Te — mo: Cau 


Resonanz tritt ein („= x), wenn e - mw? -- cqw=-0 wird, also nicht mehr bei a, sondern 
hei zwei neuen Schwingungsschnellen &, und &, (Abb. 13). Die beiden Resonanzen ent- 
sprechen den Schwingungen 1. und 2. Grades, vgl. Abb. 14. Der Dämpfer selbst führt dabeı 
Schwingungen nach Abb. 4 aus. 


Für die gedämpfte Schwingung erhalten wir eine Kurve mit zwei Scheitelpunkten ın 
der Nähe von ®, und ®,. In Abb. 15 sind einige solcher Kurven zusammengestellt. Von 
dem Ausmaß der beiden Größtausschläge y, hängt die Brauchbarkeit des Dämpfers ab. Ehe 
wir die Ausschläge selbst berechnen, müssen wir ihre Lage bestimmen. Dazu formen wir 
Gl. (19) um: 
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Abb. 14. Abb. In. 


Da der Dämpfer auf die Welle abgestimmt sein soll, ist &4 = ©... Für m »*” schreiben 
wir 2’ mo’ı = 2’mw’.—=2z”c und sondern ce ab, 


V, 
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2,—- (19a). 
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Wir wollen die Betrachtung auf das Resonanzgebiet beschränken und führen die Resonanz- 


\ ; ‚ } . R i y' 
verschiebung # = 1 — z ein. x soll klein bleiben gegen 1, Jedoch immer noch groß gegen F )- 
’ ‚ u JI 
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Bei richtiger Bemessung des Dämpfers sind stets beide Bedingungen erfüllt. Für « und ı 
ergeben sich jetzt die einfachen Formeln 


y 
a x a {7 
u = — 0 = ” : —; Be 
ÄiAn+sinin _ y? a +sın/n y 
"+ * 7 ae? . 0} 
16 16 »? 
r. . er r P) ( . . M, " 
Wir bezeichnen vorübergehend den Klammerausdruck der Gl. (19a) mit A, also Do: -K. 
C 


Um die Scheitelpunkte der Resonanzkurve zu finden, setzen wir 


Op, | ° M, KLM OK N 

02 C Fi Tr 

Die Erregung M, ändert sich mit der Drehschnelle und bewirkt, daß die beiden Resonanz- 
ausschläge verschieden groß ausfallen. Wie Versuche ergeben, kann man diesem Einfluß 
durch eine geringe Änderung der Abstimmung begegnen, so daß sich beide Resonanzausschläge 
auf einen mittleren Wert ausgleichen. Wir machen der Einfachheit halber statt dessen die 
Annahme, daß sich der Größtwert des pulsierenden erregenden Momentes M, nicht ändert, 

oM, Ö OK 


® 4 u... hd » . r 
und setzen |, 0. Damit nun \ "—0 erfüllt wird, muß auch , gleich Null werden. 
(! 2 ( zZ . 


Wir führen GI. (20) in Gl. (19a) ein, differentiieren den Klammerausdruck K und erhalten aus 
IK 
( \ 
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0 die Resonanzverschiebung 
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Abb. 16. Abb. 17. 


Abb. 16 zeigt die Resonanzverschiebung in Abhängigkeit von der Dämpfermasse. Die 
Abbildung läßt erkennen: Sehr kleine Dämpfer rufen noch keine Resonanzzerlegung hervor. 


Erst bei einer bestimmten Dämpfergröße beginnt die Aufteilung. Mit weiter wachsender 


Dämpfergröße rücken die neuen Resonanzen mehr und mehr auseinander. Die Strichlinie 
deutet die Lage des kleinsten Ausschlages zwischen beiden Größtausschlägen an. 
Unsere Betrachtung gilt zunächst nur für einen bestimmten Größenbereich von x. Lassen 


yr \° 
wir die Bedingung =-] fallen, so ergibt sich für sehr kleine .“-Werte eine gewisse Un- 


symmetrie der beiden Kurvenäste der Abb. 16. Wır können auch hier die Symmetrie wieder 


herstellen, indem wir den Dämpfer um ein gewisses geringes Maß gegen das Wellensystem 
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erstimmen: ®&4 = @ (1 + 9). Die erforderliche Verstimmung d läßt sich berechnen, sie hängt 
w \° 
1) Vs. .. g* .. » / nd y» * 
on y, z und «* ab. Für größere Dämpfer, x. P ‚ wird d9=0. Wir brauchen uns im 
re: 
lzenden nur mit dem Gebiet der größeren Dämpfer zu beschäftigen. 

Auf den Dämpfer selbst wirkt sich die Resonanzverschiebung in der Weise aus, daß 
ie Endausschläge 9; ihren Größtwert jetzt nicht mehr bei &4 = @,„ erreichen, sondern in 
er Nähe der neuen kritischen Schwingungsschnellen. Dennoch behalten die für y, ent- 
‚ickelten Formeln ihre Gültigkeit. Sie drücken nämlich stets 9, durch 9, aus, geben also das 
\usschlagverhältnis g7,:9, an. Dieses erreicht seinen Höchstwert nach wie vor bei 4 
nabhängig von dem Massenverhältnis ma: m. 


7. Die Resonanzausschläge des schwingenden Systems unter dem Dämpfereinfluß. Die 
\uswahl des geeigneten Dämpfers. In Abb. 15 ist die y,-Kurve strichpunktiert eingetragen, 
auf der bei Änderung der Dämpfergröße die Scheitelpunkte der Resonanzkurve wandern. 
Die Abbildung läßt erkennen, daß mit wachsender Dämpfergröße die Resonanzausschläge 
orst rasch, dann langsamer abnehmen. Dieses Verhalten wollen wir näher untersuchen. Wir 
setzen Gl. (20) und (21) in Gl. (19a) ein und erhalten den Resonanzausschlag 
u Me | In Cd 3 ° / 7 


7 r ie 
r lyl 2 |y(la +? —1 al+ (2x -- 2)? | ce xy" iAn+sinin 
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ww 
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Wir können jetzt M,/[e als statischen Ausschlag 9, deuten, der sich unter der statischen Ein- 
wirkung von M, einstellen würde. Dann gibt der Klammerausdruck das Verhältnis des 
Resonanzausschlages zu dem statischen Ausschlage 9,:9, an. Wir können ihn deshalb als 
„Aufschaukelverhältnis »* bezeichnen, ' 
.) 
1s = ) : =. Graysn . . (23a). 
‘ vV2xlı (a+1)”? —1-al+ — 1? —1 

In Abb. 17 ist für verschiedene #-Werte das Aufschaukelverhältnis in Abhängigkeit von der 
Dämpfergröße nach Gl. (22), (23), (23a) dargestellt. Die Kurven verlaufen anfänglich sehr 
steil, später nur noch sehr schwach gegen die m,-Achse geneigt und nähern sich asymptotisch 


einem bestimmten Grenzwert. Wenn wir die Dämpfergröße so auswählen, daß « Br wird, 


treffen wir etwa den Beginn des flachen Kurvenverlaufes. Diese Dämpfergröße ist auf den 
3 Kurven der Abb. 15 gekennzeichnet. Wir wollen sie kurz die Normalgröße m.’ des Dämpfers 
nennen und bezeichnen die zugehörige Federkonstante mit ey. Es ist 
Cd y? an +sinia mı Zm?cı (2, 

c 2n’x An m I ce 
Bis zur Normalgröße m,’ nimmt mit der Vergrößerung des Dämpfers seine Wirkung stark zu. 
Kine weitere Vergrößerung über die Normalgröße hinaus hat kaum noch Einfluß auf die 
dämpfende Wirkung: die Resonanzausschläge nähern sich asymptotisch einer unteren Grenze. 

Diese Erkenntnis bestimmt nun den Gang der Dämpferbereehnung. Man ermittelt zuerst 
die Normalgröße m,’ und untersucht, ob dieser Dämpfer die Dämpfungsleistung ohne zu 
raschen Temperaturanstieg aufnehmen kann. Tut er das, so ist my’ die geeignete Dämpfer- 
größe. Mit ma’= u.J1 ist der Dämpferdurchmesser bestimmt. Die Dämpferlänge 1-41, war 
schon durch die Wahl von 4 festgelegt, da !, durch die Eigenschnelle der Welle gegeben ist. 
irwärmt sich der Dämpfer my’ zu rasch, so ist ma > ma nach der Wärmeaufnahmefähigkeit 
zu bestimmen, wobei die aufzunehmende Wärmemenge der Dämpfungsarbeit nach Gl. (12) 
entspricht. 


8. Versuche. Abb. (18) und (19) zeigen zwei Ergebnisse von Versuchen, die im Wöhler- 
Institut in Braunschweig mit Gummidämpfern durchgeführt wurden ''). Die Versuchsanordnung 
bestand aus einer Welle mit zwei verschieden großen Schwungmassen, die von einem Elektro- 
motor angetrieben wurde und über einen Kurbeltrieb eine an einem langen Hebelarm befestigte 
Schwungmasse in hinundhergehende Bewegung versetzte. Die dadurch wachgerufenen 
pulsierenden Trägheitskräfte dieser Schwungmasse regten die Welle zu Drehschwingungen 
an, die sich der gleichförmigen Drehbewegung überlagerten. An dem Wellenende, das die 
kleinere Masse trug, wurde ein auf die Welle abgestimmter Gummidämpfer angebracht, dessen 
Schwingungsausschläge an einer besonderen Ablesevorrichtung während des Betriebes ihrer 


11) Vgl. Fußnote 2. 
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Größe und Phase nach beobachtet werden konnten. Die Drehschwingungen der Welle mi 
und ohne Dämpfer wurden durch einen Geigerschen Torsiographen aufgezeichnet. In Abh 
(15) und (19) sind die Versuchsergebnisse den errechneten Resonanzkurven gegenübergestellt 
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Die Versuchsanordnung wich von den Voraussetzungen der Rechnung insofern ab, als 
bei einem schwingenden System von zwei Massen sich auf der Welle ein Schwingungsknoten 
ausbildet, der durch die Masse des aufgesetzten Dämpfers eine gewisse Verschiebung erfährt, 
je nach dem Schwingungsgrad in der einen oder andern Richtung. Ist die eine Masse 
wesentlich größer als die andere, so liegt der Knoten sehr nahe an der größeren Masse und 
wird durch den Dämpfer kaum beeinflußt. In diesem Falle kann für die Schwingungsrechnung 
die große Masse durch eine Einspannung im Knoten ersetzt werden. Das trifft z. B. für 
Kurbelwellen von Flugzeugmotoren zu. Bei diesen besitzt der Propeller eine etwa 100-fach 
erößere reduzierte Masse als die bewegten Teile eines einzelnen Zylinders. In der Versuchs- 
anordnung hatten beide Schwungmassen gleiche Größenordnung. Die Verschiebung des 
Knotenpunktes auf der Welle wirkt sich dann in der Weise aus, daß die beiden Resonanzen 
etwas näher zusammenrücken. In der Rechnung wurde die Knotenpunktverschiebung ver- 
nachlässigt. Abb. (15) und (19) zeigen deutlich den Einfluß dieser Vernachlässigung. 245 
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)ie Wirkung des Enddruckes auf die Biegungsschwingung 
eines Stabes mit innerer Dämpfung. 
Von K. Sezawa in Tokyo. 


1. Einleitung. Im vorigen Jahre hat Herr Muto') in dieser Zeitschrift eine Arbeit 
eröffentlicht, welche mit der inneren Dämpfung des schwingenden Stabes verknüpft ist. 
r hat eine Differentialgleichung, die von der von Holzer?) früher gegebenen Gleichung 
twas verschieden ist, unter der Annahme eingeführt, daß die innere Reibung des Stabes 
ich als lineare Funktion der Deformationsgeschwindigkeit verändert. Ich hatte aber bereits 
‚or einigen Jahren eine gleichartige Gleichung aufgestellt. Die Form der von mir erhaltenen 
‚leichung ist 
A 


0’ Y Vy,. y 
‘ A 2 . &e 1.2 ! 
BRENNEN: 4 ea 


- ist die Koordinate der Mittellinie des Stabes, y die dazu senkrechte Ausweichung ihres 
Punktes zur Zeit t, o die Dichte, E der Elastizitätsmodul, £ der Koeffizient der inneren 
/ähiekeit und k der Trägheitsradius des Querschnittes des Stabes. Es scheint mir, daß diese 
(Gleichung zuerst von mir ermittelt worden ist. Ich hatte jedoch die eingehende Erklärung 
dieser Gleichung unterlassen, denn man kann ihre Bedeutung ohne Schwierigkeit verstehen. 
Bald nachher hat Professor Suyehiro*) eine Ableitung für die Entstehung der Gleichung 
und ihre Anwendung gegeben. 

In der letzten Zeit habe ich ein Problem der gedämpften Biegungsschwingung unter 
Wirkung des Enddruckes untersucht. Dieses Problem ist deshalb wichtig, weil die Umstände 
der Schwingungserscheinungen bei der Wirkung des Enddruckes außerordentlich geändert 
werden, und überdies, weil die Schwingung des durch Enddruck beanspruchten Stabes mit 
dem Probleme der Knicksicherheit einer Säule eine wichtige Beziehung hat. 


2. Schwingungsgleichung. Das Problem, die Wirkung des Enddruckes auf die unge- 
dämpfte Biegungsschwingung des Stabes zu finden, hängt mit der gegenwärtigen Untersuchung 
zusammen, und die partikulären Fälle des Problems wurden von Cowley und Levy, 
Howland und mir untersucht’). Die Formulierung der Differentialgleichung der Biegungs- 
schwingung mit innerer Dämpfung rührt von mir‘) her. Wenn die Absicht dieses Versuches 
auch ist, wesentlich die Differentialgleichung auf die Lösung der Schwingungsprobleme an- 
zuwenden, werde ich doch zuerst den Prozeß, die Differentialgleichung theoretisch abzuleiten, 
darstellen. Wenn man die von der Elastizität und der Viskosität abhängenden Scheerungs- 

7) 
\j bezeichnet, dann ist die lineare 
( 
Verschiebung der Längeneinheit des Stabes durch die Gleichung 
0y OF o0F 
2: r 
of 0x Ox ot 


.g . nn am . . y ( 
kräfte ın der Fläche des Querschnittes x mit F und 


O0 A 2) 
vegeben; dabei ist A die Flächeneinheit des Querschnittes. Die Gleichung für die Winkel- 
änderung der Längeneinheit des Stabes lautet 


,„0y oO 0’y OF’ IY 


e* : „il ı + E, er \ 6 "+ +7 
Zul oro« Ba 0.x” a lainat! Of HA 


(3). 
. 0% 


wo T der Enddruck ist. Wenn man (3) in (2) einsetzt, dann ist die Gleichung des schwingen- 
den Stabes 


Be... Se. on To’y 9y 
0 + ER’ pre i ‚= ok’. r (4). 
“or O7 ofod«x are *+7 37037 
Wenn der Einfluß der Drehungsträgheit vernachlässigt wird, folgt die Gleichung 
"Yun 4,094 , T®y_o ; 
0 so tur Tech . a: er ee De 5 Ta. A 
"of 0 .r? otdx! Ada? 
1) K. Muto, Diese Zeitschr., Bd. 10 (1930), S. 346. 
®2), H. Holzer, Diese Zeitschr., Bd. 8 (1928), S. 272. 
3) K.Sezawa, Bull. Earthq. Res. Inst., Tokyo, Bd. 3 (1927), S. 50. - K. Sezawa, Rep. Aeron. Res. Inst., 
Tokyo, Bd. 4, Nr. 45 (1928), S. 91. 
#4) K. Suyehiro, Proe. Imp. Acad., Bd. 4 (1928), S. 263. — K. Suyehiro, Bull. Earthq. Res. Inst., Tokyo, 
Bd. 6 (1929), S. 68. 
5) W.L.Cowley u. H. Levy, Proc. Roy. Soe., London, Bd. 95 (1919), S. 440 bis 457. K. Sezawa, Jour. 
\eron. Res. Inst., Tokyo, Nr. 5 (1924), S. 39 bis 46 (in japanisch). R.C.J. Howland, Phil. Mag., [7!, Bd. 1 (1926), 
S. 674 bis 704. 
6) K. Sezawa, Rep. Aeron. Res. Inst., Tokyo, Bd. 4, Nr. 45 (1928), S. 107. 
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Die letzte Gleichung ist diejenige, welehe ich vor drei Jahren abgeleitet habe. Wenn T in (5 
null ist, stimmt die resultierende Gleichung mit dem von mir früher gefundenen Resultate 
überein. 


Um die Gl. (4) zu integrieren, nimmt man einen Zeitfaktor von dem Typus ei»! an. So 
erhält man den Ausdruck 





r. ” 
d’y BE; d?y op” 0 ’ 
4 7 ur ei - u )) 
d.r? F k? 1 ip£k? dx? E Jh? u ) p Eh? J 
Schreibt man in (6 
<= \ 
I 222 
oD’ MW \ M) I2 
4 0] au _ OP P i _ a 
a zur‘ nn (5) ! Ysın/iz . . (9. 
Eh? -+ipek® 5 E+ipe® TV 9 ö 
dann ergibt sieh 
I )* 
an | 
ee — te TUT ee AR 
JE 1° | 
[u 
Die Lösungen dieser Gleichung sind u 
# \ | | a + pP? Ad 4; \ | | oe ni pP? ta 
> r E I > N er u (11): 
A3 l "| I 
wobei I die Länge des Stabes ist. Das allgemeine Integral von (4) lautet W 
a —+-h’—a 6 a+p+a 
Y | Acos }) | e—+ B6o!} VVetnH X u. 
r e 
(12). 
ch a Aa ds a? "ta 
+ ( a ER 2+D&in! ud x\eint 
I I 
Durch Einsetzen dieser Formel in die Randbedingungen bestimmt sich die Beziehung 
zwischen dem Enddrucke und der Schwingungsperiode des Stabes. In den folgenden Ab- 
schnitten werden einige Beispiele gegeben, in denen nur die Hauptschwingung abgeleitet wird 
weeen der Wiehtiekeit für das Problem der Knieksicherheit der Säule. 
3. Die an beiden Enden gestützte Säule. Wenn die beiden Enden des Stabes gestützt 
sind, erhält man die Randbedingungen: 
\? 
O7, 
e=H+  y=d, En ee ee A 
2 J 0 .r? m 
der Wert von . wird von dem Mittelpunkte des Stabes an gerechnet. 
Die Kinführunge von (12) in (13), indem B, C und D in (12) null gesetzt sind, liefert 
fe) 
einen Ausdruck für die niedrigste Schwingung, nämlich 
ru a en an a a nt 
und, wenn man die Werte von a und p°? in obiger Formel aus (7) und (8) einsetzt, geht dies 
über ın u 
(1 . ‚.ntER pe Ek? n? 4 0 15) 
o»p° r 2” — I!» - )). 
O0} !: / j* /’ I? A ( 
Hieraus folgt 
ine n? (ER? \ rn j 
DE ee Be eg - y u: KU, 
2o(F + nk’) o’+a?k’)\ 7 1 J t o’(i?+ nk’) 1 
so daß die Gleichung der Ausweichung des Stabes, wenn die Einwirkung der Drehungs- 
trärheit vernachlässigt wird, lautet H 
Ä | ER Yen; "ii 7 ah, BR: in 
y==.|exp. = . cos 2 ). 
| 20" )/sinff oP\ PP ıj to? 1° | 





| 12, Heft5 « u; a h an 
*h ee 1939 Sezawa. Die Wirkung des Enddruckes auf die Bierungesschwingung 2717 





Es stellt sich heraus, daß der Dämpfungsfaktor nieht von demjenigen des Stabes ohne 














te nddrucek verschieden ist, daß aber die Schwingungsperiode in diesem Falle außerordentlich 
eändert wird. Wenn 
So T 7? El? nı® Zu k* 
1 z 401° 
ndet eine aperiodische Schwingung statt; und, wenn 
6). T n®’ElR? 
: er BE Be Gage, an 2 De) De 
A I? 
vird die Säule labiıl. 
) 4. Die an beiden Enden eingespannte Säule. Wenn die beiden Enden des Stabes ein- 
respannt sind, lassen sich die Randbedingungen 
l 0 ? 
a - yzyd, I u, re N a: 
2 0x 
0). 
schreiben. Durch Einsetzung des Ausdruckes von y aus (12) in diese Bedingungen, indem € 
und D in (12) null gesetzt werden, ergibt sich die Gleichung: 
[7,5 = a’ —+ ° a = — em ah’ +n 
1); | va’+rp’—a et — +] Ve’+p+targ | z | sr 
wo die Werte von a und ß° schon in (7) und (5) angegeben sind. 
Es ist sehr schwierig, die Gl. (21) nach p zu lösen, und daher habe ich die folgende 
mittelbare Methode gewählt. Die aus der Gl. (21) sich ergebende Beziehung zwischen « und /°? 
9) wird in Abb. 1 dargestellt. 
cccc A Dame Bang ME ' 
Pr | | | 
Ing MT | | 
| 
Ab- 800! 
rd 
700 
600, 
ıtzt | 
5o0\ 
400 
13), | 
300 | 
200 | | 
* | | I 
fert 200 | 
| | | 
0 | | | l rc 
20 WM AOC—H{3GO 4 BE TR E20 
14) [RAoa21] & 
| Abb. 1. 
lies 
Hieraus hat man näherungsweise die Formel: 
15). | 
Pi eh, a en ne 
| 19,4 
Durch Einsetzen von (7) und (8) in diese Formel ergibt sich die Gleichung: 
16), ERED RR: 
LE Be . 00ER? ‚.,[394EK’ 1 de 
op? (P 412,74) ip 12,74 | ” > Sn 
1g8- 
Hieraus folgt 
500818 12.74 3941ER  T 500° 227 r 
rm ) [77 nd “ 1} + 2} und ” | . . 677 und . .) = |, 
1). 1 20(? + 12,74 k?) 1? o(l?+ 12,74 4°) 1° | to?(1? + 12,74 K?)? 
Y r 
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so daß die Gleichung der Ausweichung des Stabes, wenn die Einwirkung der Drehungsträgheii 
vernachlässigt ist, lautet 


500 k? ’ cos] 1» 14 ku Ek? T | 00° €? 1° ] 





WEETYL i _ } a 

Se sin] ol I? A to?”l 
| | (25) 

/ 2 02 ER at Ata 
a ” +Pp " x BGof) 4 1 x 
g I I 
wo a und 5° bestimmt werden durch 
\ ep 

2 — . da ‚2. ‚BP ı 10) 
p 500 [1 197 a); sam AH © (26) 


indem der Einfluß von ipE& auf die Deformationskurve des Stabes vernachlässigt wird. Wenn 


T 3934ER? 300° 9" (27) 
rue f -s), 
A pP 12.74 x 4 ol" 
wird die Schwingung aperiodisch; und, wenn 
T 394 ER? 
HE... = (28) 


l 1° 


wird die Säule labil. 


5. Die an einem Ende gestützte und am andern eingespannte Säule. Die Randbedingungen 
in diesem Falle sind 


7, 
"—0) yo 3,0 
(29), 
7] 
d : == U) i 0 
a 1; Y ) dr 





indem man das gestützte Ende als den Anfangspunkt von x wählt. 
Wenn man die Werte von A und B im Ausdruck (12) null setzt und die resultierende 


Gleichung in die Bedingungen (29) einführt, erhält man 


D 
N > 


VYva+Pp—actgl ver —a yva+p+ a Cte ] vatpfta=0. . . (8). 


Die aus dieser Formel sieh ergebende Beziehung zwischen a und 5° wird in Abb. 2 dargestellt. 











Hieraus hat man die Näherungsformel: 


| 


2? __OINR 
’» — 238 [1 1041 [2 


ee 


Durch Kinsetzen von (7) und (5) in diese Formel ergibt sich die Gleichung: 


DIS Eh? I)2 ER? T 


op(P + 1LSR®) ip”, 11,79” j: 7 EEE >>} 


Iın 


wi 


H« 
W 
M« 
Da 
kri 


St3 


Me 


\h 


Zar 





1 
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‚ erhält man den Ausdruck : 


... 
ww 
.. 
ww 
—_ 


I23S Eh? Pe} 11,5 20,2 Ek? T 2387 h* 
pP 20(° 11,3 k?) 7? | ol + 11,574?) I? A to’(+11,8A? PT 


ie Gleichung der Ausweichung des Stabes, wenn die Einwirkung der Drehungsträgheit ver- 
achlässigt ist, lautet 





238° k? , cos] /11,8 ( DOD2EK? T\ 2382” 
y=e 21 IA . j Di} D) os 
J sin} 0 dä I’ A to®l 
(34), 
s ” a, a? ” 
<ICsin| Are 44 DSin! Far © .| 
| | 
‚o a und 5° bestimmt werden durch 
1 U; 
” 238 [1 . a): 20 = A I 5 0: MA 
P 101 = A EM gi 
Wenn 
T 22ER? 238° &°K* ar 
= . (3b) 
A 1? ILS X 401" 
\indet eine aperiodische Schwingung statt; und, wenn 
T 22ER? 2 
- 4], 


A j 
wird die Säule labiıl. 


6. Zahlenbeispiel. Die Werte von & für einige Materialien wurden bereits von Professor 
Honda und Dr. S. Konno’) durch Experiment gefunden. Mit Benutzung dieser beobachteten 
Werte von E habe ich die Fälle berechnet, in welchen das Material des Stabes aus Stahl, 
Messing oder Aluminium besteht und die Länge und die Dicke desselben 50 em und I em sind. 
Damit können wir uns einen Begriff von der Größe der kritischen Schwingung und auch der 
kritischen Knicksicherheit des Stabes machen, 

(i) Gestützt— gestützt. (il) Eingespannt—eingespannt. (iii) Gestützt eingespannt. 








i | ) (II) Belastung f. krit. (I) — (ID 
E & (I) Belastung f. krit. Dr: 
Metall 0 er I 4 x z Knicksicherheit (1) 
(C.G.S.) (0.6G.8.)  Scehwingung (ke/em?) nr n 
(ke/em?) un 
BEER area (i) 670 35 i) 670 0.52 
7.8 2x 1012 5% 108 (ii) 2670 — 72.4 (öi) 2670 2,71 
(iii) 1370 — 17.5 (iii) 1370 1.28 
Messing .... (Ü) 336 0,29 (Ü) 8336 0,09 
8,5 1 x 1012| 1,5 X 108 (ii) 1340 6.00 (ii) 1340 0.45 
(iii) 690 1.45 (iii) 690 0.21 
\luminium . (Ö) 235 0.26 i) 2535 0.11 
27 |0.7X 10120,8% 10% (ä) 940 — 535 (ii) 940 0.57 
(iii) 472 1,29 (iii) 472 0,27 


Es stellt sich heraus, daß die Schwingung aperiodisch wird, wenn der Enddruck des 
Stabes irgendeinen Wert unterhalb der Knicklast annimmt, und auch daß, je mehr die Enden 
les Stabes eingespannt sind, desto kleiner die der aperiodischen Schwingung entsprechende 
Belastung wird. Dieses Ergebnis hat eine Beziehung zu dem von Professor Suyehiro°) 
egebenen Satz über die kritische Länge der Schwingungswelle eines Stabes. Wenn die Haupt- 
schwingung auch allein im vorstehenden Aufsatz abgehandelt ist, so lassen sich die höheren 
ıarmonischen Sehwingungen des Stabes doch leicht in der gleichen Weise untersuchen. Es 
nuß jedoch bemerkt werden, daß die Schwingungen eines Stabes von äußerst kurzer Wellen- 
änge durch Einführung einer Methode *) diskutiert werden, welche von der Anwendung der 
zewöhnlichen Differentialgleichnng des schwingenden Stabes wesentlich verschieden ist. 203 

’") K. Honda und S. Konno, Phil. Mag., Bd. 42 (1921), S. 115 bis 123. 


s,K. Suyehiro,|.e. 
9») K. Sezawa, Verh. 3. int. Kongr. f. teehn. Mech., Teil 3 (1930), S. 167 bis 172 
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Berechnung der Oberflächentemperaturen geometrisch ein- 
facher Körper bei Abkühlung oder Erwärmung für den 
Bereich der Isolierstoffe. 

Von @. Pösehl in München. 

Mitteilune aus dem Laboratorium für technische Physik der Technischen Hochschule München !). 


Mathematische Umformune der Formeln; Reehnungsgang für Abkühlung einer Platte; Rechnungsergeb 


nisse und Anwendungsmögrlichkeiten der berechneten Kurven.) 


nter der Rationalisierung der Arbeit versteht man das Bestreben, mit einem Mindest 
Ü aufwand von Arbeit, Zeit und Geld einen möglichst günstigen Wirkungsgrad zu erzielen. 
Von ihr wird in der Praxis vielfach Gebrauch gemacht, und zwar oft so weitgehend, daß das 
Ergebnis der rationalisierten Arbeit den Eindruck der übertriebenen Eintönigkeit erweckt. 

Bei wissenschaftlichen Arbeiten findet die Rationalisierung noch verhältnismäßig selten 
Anwendung, obzleieh die Anfänge dazu schon weit zurückliegen, denn z. B. die Logarithmen- 
tafeln, die trigonometrischen Tabellen und andere Funktionentafeln entsprangen dem Gedanken. 
denjenigen, welche diese Funktionen bei ihren Berechnungen nötig haben, die Arbeit dadureh 
zu erleiehtern, daß man solche Funktionen einmal ausrechnete und sie dann im Druck der 
Allgemeinheit zugänglich machte. 

Bei der fortschreitenden Verwertung theoretischer Überlegungen in technisch physika- 
lischen Arbeiten mehren sieh die Funktionen, deren man bei den Berechnungen bedarf. In 
diesem Zusammenhang hat bereits H. Gröber?’) Funktionentafeln und Diagramme veröffent- 
licht, mit denen man imstande ist, die Temperaturen sich abkühlender Körper in einfachster 
Weise zu bereehnen. 

Gröber behandelte die zeitliche Temperaturänderung, welche Kugel, Zylinder und Platte 
bei Abkühlung oder Erwärmung erfahren, und zwar auf der Oberfläche und in der Mittte. 
Außerdem bestimmte er noch den Wärmeaustausch der drei Körper mit einer Umgebung von 
konstanter Temperatur. Dieses Problem war zur Auswertung von Funktionentafeln im Sinne 
der oben erwähnten Rationalisierung deshalb besonders geeignet, weil die Temperaturberechnung 
nur mittels unendlicher Reihen möglich ist und einen großen Zeitaufwand erfordert. 

Im besonderen ergab sich dabei, daß die Konvergenz dieser Reihen bedingt ist durch 
eine Kenngröße, die der Wärmeleitzahl und der Zeit proportional ist. Die Konvergenz wird 
um so schlechter, je kleiner der Wert der Kenngröße ist, also je kleiner die Wärmeleitzahl 
und je kürzer die Zeit. 

(‚röbers Rechnungen beziehen sich auf verhältnismäßig große Werte dieser Kenngröße, 
bis hinab zu einem Wert von 0,01. Infolgedessen besaßen seine unendlichen Reihen eine 
verhältnismäßig gute Konvergenz. Die nachstehend berechneten Tabellen ergänzen die 
(röbersehen Reehnungen im Gebiet der kleinen Werte der Kenngröße und sind daher 
anwendbar auf Materialien von kleiner Wärmeleitzahl, also für die Wärmeschutzstoffe, und 
auf den Beginn des zeitlich veränderlichen Wärmeaustausches. Diese Ausdehnung bedingte 
eine äußerst umfangreiche Rechenarbeit, weil die unendlichen Reihen wegen ihrer schlechten 
Konvergenz teilweise bis zu 30 und mehr Gliedern ausgewertet werden mußten. 


A. Mathematische Umformung der Formeln. 


Es empfiehlt sich wegen der schlechten Konvergenz der unendlichen Reihen die von 
(röber anzerebenen Formeln so umzuformen, daß die in ihnen auftretenden trigonometrischen 
und Besselfunktionen nach Möglichkeit eliminiert werden. Dies ist übrigens auch bei guter 
Konvergenz der Reihen von Vorteil und gibt den Formeln gleichzeitig ein wesentlich an- 
schaulieheres Bild. 


I. Platte: Für die Abkühlung bzw. Erwärmung einer Platte von der Dieke X unter 
Annahme einer einachsigen Wärmebewegung (Wärmeströmung nur in Richtung einer einzigen 
Koordinate) gibt Gröber für die Temperatur an der Oberfläche ©, folgende unendliche 
Summe an: 


4 . et 
Y sın ayR IR? r., f 
9,=0. I 2 — e Be 5 a OR Te 
1 9% + sın 07; Co8 d;; | 


') Herrn Geh. Reg.-Rat Prof. Dr. Dr. Inge. e. h. Oscar Knoblauch, meinem hochverehrten Lehrer, bin ich 
für sein unermtüdliehes Interesse und seine freundliehe Förderung zu herzliehstem Dank verpflichtet. 
-), HM. tröber: Die Abkühlung und Erwärmung geometrisch einfacher Körper, Z.d. V.d.]1. 69. (1925), S. 70», 


\ 


ni 
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iierbei bedeutet a die Temperaturleitzahl: «a — ), = Wärmeleitzahl, ce = spez. Wärme, 
spez. Gewicht, t= Zeit. 


Diese Formel (1) gilt, wenn die Platte mit der überall gleichen Übertemperatur (bzw. 
U'ntertemperatur bei Erwärmung) ©. plötzlich in ein Medium kommt von konstanter Temperatur, 
‚lie wir als Null wählen können, da wir nur Temperaturunterschiede betrachten. 


Die in der Formel auftretenden ö, sind Wurzeln der transzendenten Gleichung: 


DERBER FEED . alien nn A 


A r .. hd BET} .. ET} . 5 . 
wenn h= . das Verhältnis von Wärmeübergangeszahl « zu Wärmeleitzahl 4 bedeutet. Die 
/ 
anzerebene Gl. (1) gilt für einseitige Abkühlung der Platte. Bei beiderseitiger Abkühlung 


D iin. a 
setzt man „, statt X in die Formel, wobei dann D die Dicke der Platte ist. 
Man benützt nun Gl. (2), um in Gl. (1) die trigonometrischen Glieder zu ersetzen. Aus 
(Gl. (2) folgt: 
Ö,; Sin dj, 


COS 0% = h X 
und eingesetzt in Gl. (1) erhält man: 
AyR: al 
I so. d.2 d r ’ ‘ X2 
4, 9, No» O7; sın Ök ia F 7 v2 9, 3 a j Ä 
ki eg sın" 07;\ ‚| 
Nö, -+ Ör X\——. +; 
h \[ ON h X a2 sin’d;, AhX, 
| 
Da nun —=1-+ cig’ 
sin?’ 
dr 
und aus Gl. (2) etg? Od; 7 Y 
so ergibt sich 
| d2 at dı2 at 
v Je ze .y eh > 
oJ, = 6 - 4. > h X an = - ge 
I; | ] Y | Rn" 0%" | 0% u h A (1 4 h AN ) 
en HXNETUNX 
„al 
, O2 T, 


. y E . ’ 
Ersetzen wir den Ausdruck 2h 2 durch die Abkürzung P,, so ist es 


k-109v +hX(l+hX) wi 
also möglich, die Oberflächentemperatur ©, zur beliebigen Zeit f aus der Anfangstemperatur ©. 
mit Hilfe der Funktion P, zu berechnen. 


al 

" O1? . 

a we _ 
D—=2hX) : I N 
0 l En Ö% h % (A h X) ( ) 


Der Vergleich der Gl. (3) mit (1) zeigt deutlich die Vereinfachung und Vergrößerung der 
Anschaulichkeit, welche durch die Beseitigung der trigonometrischen Funktionen erreicht 
worden ist. 





Die große Zahl der Veränderlichen (Wärmeleitzahl, Wärmeübergangszahl, spez. Wärme, 
spez. Gewicht, Dicke X, Zeit), von denen ®, abhängt, können zu nur zweien zusammen- 


» » .. ® r ..g at r . 1 . . 
vefaßt werden, nämlich zu den Kenngrößen v: und AX, denn 6, ist nach Gl. (2) eine reine 


. 


Funktion von hX. 


In ähnlicher Weise lassen sich die Temperatur der Mitte der Platte 9, und die bis 
zur Zeit £ von der Oberfläche F der Platte an die Umgebung übertragene Wärmemenge 
durch die Funktionen P,, und Y ausdrücken, so daß 


In = I. Din 
und 
g9=FXcyG,P, 


w 
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wobei ®,, 
eebracht wurden: 


4 04? 4 
. Y > “L“ l 
Dm=2hX, ABER | J 

| nr 102 +hA(l-+hX) cos 9; (+), 

af 

dj . e 

A > . X2 

P-2(hX) L Ian (5). 


._‘ 0,” 107° ur dh X)”] E= 0; h X 


2. Kugel: Durch einen analogen Rechnungsgang bekommen wir für die Oberflächen- 


temperatur der Kugel vom Radius AR aus Gröbers Formel 


‚_SinP; —r; COSYvE —’kK? 75 sin vj. 
Oo, , \ ) I, I, } P R / +. D, 
u . ° sın v,, COS vr, ”] 
durch Umformung 
r at 
Y € 2 ni 
DD =2hR z a en 
= 2° hRÜ-hR) (0) 
En, ar. 0 99 SEEN BB. Wurzel der Gleichung sind: 
rel Alt . >: 5 Er rn ae 


Auf ähnliche Weise bestimmt sich, ausgehend von Gröbers Formeln für die Temperatur in 
der Mitte: 
iR — 9, «D 


m 
„at 
vi 5.” R? 
2) e v]; 
ö.=2hR) — BER. 
u. re —hR(i—hR) sin v, 
und für den Wärmeaustausch: 
| 
0=— R’ncy9,. 
.) 
12 af 
4 | e MR 
W—6(h R): - r ’ 
NN = yo —hR(—hR)) 


3. Zylinder: Bei der Abkühlung des Zylinders vom Radius A können wir für die Tempe- 
ratur an der Oberfläche, ausgehend von der Formel: 


L at 
v ri “ HU? 5, 
,=9.1,=9. So | I (ur) r RB? J (up). 


IN l 4 I (ur) J (ur) 


worin ./, und ./, Besselfunktionen erster Art 0. und 1. Ordnung sind, durch Zuhilfenahme 


der Gleichune: 


MEN 2-4 50 ee ee rn (5), 
deren Wurzeln die u, Hss..... 4... Sind, schließlich auf die Form kommen: 
2 at 
2 u ?7 
, e 
BD =2hR, Ra Te ee Ra: 


1 ' U“ B- (h RR)” 


Durch Elimination der Besselfunktionen mittels der Gl. (S) ist es also möglich, der 
Funktion ?, eine vereinfachte und ähnliche Form zu geben, wie den Funktionen ®, nach 


Gl. (3) und (6). In entsprechender Weise gelangt man durch Elimination der Bessel- 
funktionen zu folgenden Gleichungen für die Temperatur in der Mitte: 


In . Du . 
„at 
ve ”w | 
Don 2 hı BZ 


1 | 1," 1 (h RR)” J (1%) 


und Y durch Umformung der von Gröber angegebenen Formeln auf die Form 


- 


il 


e 


n- 
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ınd für die ausgetauschte Wärmemenge: 
0=R’nLcy9.!V. 


al 
, uk? 
Zn R: 
re 


vo Ar? (up? + (hR)P) | 
B. Rechnungsgang für Abkühlung einer Platte. 


Aus den von Gröber veröffentlichten Kurven und Tabellen läßt sich leicht erkennen, 
daß die von ihm errechneten Werte für die Temperatur in der „Mitte“ und für den „Wärme- 
austausch“ auch für kleine Werte der Wärmeleitzahl / und damit für kleine Werte der 

r fi . . . rn\ . . Min 
Temperaturleitzahl «a: je ausreichen, da sich die Temperatur in der Mitte und der Wärme- 


ME 
" ne a : ’ ed ad 
austausch erst nach längerer Zeit merklich ändert, so daß der Wert der Kenngröße ve bzw. Rn: 
RW u 


für praktische Fälle in dem von Gröber berechneten Bereich liegt. Für die Oberflächen- 
temperatur aber ist die Änderung auch bei sehr kleinen Zeiten schon so groß, daß die 
u a ’ 2 a af 
kleinsten von Gröber angenommenen Werte der Kenngrößen ve 0,02 bzw. RR: 0,01 stark 
4 L” 
unterschritten werden, wenn wir den anfänglichen Verlauf der Oberflächentemperatur unter- 
suchen wollen. 
Im folgenden sei der Rechnungsgang für die Bestimmung des Verlaufs der Ober- 
» .. . .. e . ä jr a r 
[lächentemperatur der Platte in Abhängigkeit von den Kenngrößen v: und AX kurz an- 
zedeutet, der auch die Schwierigkeiten und den großen Zeitaufwand für die Berechnung er- 
kennen läßt. 


Aus der Formel 


t 
dr? a 
k“xa 


Oh X S © $ 
uEnd e ERLERNTE ME ET ae Be ee ee 


erkennt man ein starkes Abklingen, also eine gute Konvergenz der Reihe, wenn der Einfluß 


der e-Funktion groß ist, d. h. für große Werte von ve' Da die Reihe aber gerade für kleine 


b at ' em : 
Werte von y; Zu berechnen war, mußte eine große Anzahl von Gliedern der Summe bestimmt 
werden, damit die Vernachlässigung der Restsumme keinen Einfluß auf die gewünschte Ge- 
nauigkeit mehr hat. Um nicht zu wenig, aber auch nicht zu viel zu rechnen, muß man die Rest- 


a r . 
summe für jedes Y: und AX abschätzen und aus der Abschätzung die Anzahl p der Summen- 
glieder bestimmen, die für eine geforderte Ge- 
nauigkeit der Summe notwendig ist. (Vgl. die D 
folgenden Abschnitte) Aus der Genauigkeit ru 


N # > | 
der Summe und der Anzahl der Summenglieder | 
folgt die Genauigkeit des Einzelgliedes, das die | 

, | 
| 
| 











Form hat: 
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somit die erst noch zu berechnende Größe d, enthält, deren Genauigkeit durch Differentiation 
dieses Einzelgliedes g,, nach ö, gewonnen wird. Mit der sich ergebenden Genauigkeit ist 
damit die notwendige Dezimalstellenzahl jeder Wurzel gegeben. Dann können durch ein 
Näherungsverfahren die p ersten Wurzeln der Gl. (2) ösind=-hXcosod bis auf die notwen- 
dige Genauigkeit bestimmt und in die Formel eingesetzt werden, um nun die Reihe numerisch 
zu berechnen. Diese Rechnung ist für verschiedene Werte von RX zu wiederholen. Man 
erhält dann eine Kurve für = PD, mit IX als Abzisse bei konstanten ” als Parameter 
C 4 


'vel. Abb. 2). 
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1. Abschätzung der Restsumme. Wenn man die einzelnen Glieder der unendlichen Summen ind 


als Ordinaten in Abhängigkeit von d, als Abszisse aufträgt und die Endpunkte durch eine 
Kurve verbindet, erhält man für den Bereich der Restsumme ungefähr folgendes Bild. 

Die Funktion sind h Xcosd hat die Eigenschaft, daß der Unterschied zweier auf- 
einander folgender Wurzeln sich um so mehr dem Wert z annähert, je größer k ist, so daß man vol 
im Bereich der Restsumme für deren Abschätzung die Abstände der aufeinander folgenden en 
Werte 9, alle gleich = setzen kann. Wie aus der Zeichnung ersichtlich, läßt sich dann die 
kestsumme 


af 
“ OK? 
X: 
Ih X > Er 
rer thAl+ BR 
dureh das Flächenintegral auf folgende Art abschätzen: wi 
BR z „at 
’ 02 = . 02 r. 
\ onX N E X? 2hX E N? 10 
7% = n x. 
u 1 i N O, .. Jh N\ ( | 1 In = ) 7ı 0" == hı N (1 u hı 2 ) 
hi) | ' oT% 
» r +) 
> ra 


Die Fläche zwischen Kurve und Abszisse wird durch Streifen von der Breite z ersetzt, wo- 

durch man das Restelied &,, also die Summe der Ordinaten, durch Division des Flächen- 

integrals durch 7 erhält. W 
Um zu einem integrierbaren Integral zu kommen, ersetze ich das Integral unter ver- Gi 

schiedenen Substitutionen, die leicht verifiziert werden können durch ein sicher größeres, 


z 4 * .. ‚Ir A ) Pr r . 
wie ım folgenden angedeutet, unter Benutzung der Abkürzungen a für und A für X (I+AN). 


X? 
M 
- R: T- \ 
\ 3hXl e-?dö D2hX[ e-e:+eP 3hXe? ([e-% 
u | 7 = ıLß T \ 2:21 2 ) ») 1 E ) \ 2 dz \ 
Rn =/ na „asyae-p Say] 2 ei 
Oy Ten P I? 3 
h 
daraus 
I. 
v h X e“P \ e "du 
m | 7 OÖ) f N u 
fi du? an» ( 
e . ” * ” ” 
Das Integral \ du ist ein bekanntes und tabellarisch schon ausgewertetes bestimmtes 
A . 
( 
A 


Inteeral ’). 
— 


Um für $, die für die Praxis hinreichende, auch von Gröber angenommene Genauig- | 
keit von O,01 zu erreichen, muß der Fehler von P, kleiner als 0,005 sein. Rechnet man nun 
für die Vernachlässigung der Restsumme 3, die Hälfte des Gesamtfehlers, dann kann die 
andere Hälfte des Fehlers der gerecehneten Summe selbst zugeteilt werden, die ja durch die 
Ungenauigkeit der Wurzeln ö, auch nicht vollkommen genau berechnet werden kann. Die 
vernachlässigte Restsumme muß demnach kleiner sein als 0,0025. Aus dieser Bedingung läßt 
sich die Anzahl » der zu berechnenden Glieder bestimmen. Es ergab sich dabei, daß z. B. 


° ° ‚ nz [2 £ . ‘ » Y ® 
bei einer Kenneröße 0,0003 bis zu 36 Summenelieder zu berechnen waren. 


X? 
2. Erforderliche Genauigkeit der Wurzeln ö,. Die Genauigkeit des Einzelgliedes ergibt sıch, 
0.0025 a u “ 
wenn die Anzahl der Glieder gleich p ist, zu dgyı gr Durch Differentiation des Einzel- 
s ) 


yo, nach Ö, läßt sieh die erforderliche Genauigkeit dö, und damit die gewünschte 
Dezimalstellenzahl erhalten: Aus 


eliedes q 


folet 
n ») ») ao 3 2 +.) > [44 Rh) „2 
Ph _ a ET HP) 2er 


do, (0,” + pP)’ 


') Jahnke-Emde: Funktionentafeln mit Formeln und Kurven 8. 21 (Integral-Logarithmus) Leipzig-Berlin 


1000, Teubner. 
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nd daraus 


] 5 j dq 0 (0%: — p)” 
ur Ih X09,e ed? (1 + ad PB) 


0.0025 
vobeı dag ok p 


‚nau berechnet werden. 


zu setzen ist. Im ungünstigsten Fall mußte $, bis auf 5 Dezimalstellen 


3. Berechnung der Wurzeln der Gl. (2) 


sınd=hXN Cosd. 


Auf graphischem Wege oder durch Schätzung erhält man den ungefähren Wert von dj. 
Wenn man Gl. (2) in der Form 


Ö 
h X 


ete od 


»raphisch darstellt, sieht man, daß man die Wurzeln als Abszissen der Schnittpunkte der Ge- 
Ö r R 

raden 9 = und der Kurven y==ctgö erhält. 

« u. « 

Die Verbesserung der entweder graphisch oder durch Schätzung ermittelten Werte der 
Wurzeln 6, erfolgt durch ein Näherungsverfahren. Es hat sich hier als praktisch erwiesen, 
(1l. (2) in logarıthmischer Form anzusetzen: 

log h X — log ö,, = log tg d).. 
Man bildet nun die Differenz auf der linken Seite der Gleichung, indem man den geschätzten 
Wert Ö, einsetzt und erhält aus einer trigonometrischen Tabelle für diese Differenz unter 
lor te d, ein genaueres Ö,;, das man nun wieder in die Differenz einsetzen kann, um daraus 
einen noch genaueren Wert für ö, zu bekommen. Dieses Verfahren kann beliebig oft wieder- 
holt werden. Meist genügt eine ein- bis zweimalige Wiederholung. 


af 


Für die Oberflächentemperatur von Kugel und Zylinder für kleine Werte von RE 
L 


war 


der Rechnungsgang ähnlich wie bei der Platte. 


C. Rechnungsergebnisse und Anwendungsmöglichkeiten der berechneten Kurven. 


Für den Verlauf der Oberflächentemperatur bei Abkühlung und Erwärmung wurden in 
Ereänzung der Gröberschen Kurven je fünf weitere Kurven berechnet, die im Bereich von 


a af 
x: bzw. np: Zahlentafeln 1 


bis 3 geben die Werte für die Platte, den Zylinder und die Kugel wieder. In den Abb. (2 bis 4) 
sind die errechneten Werte graphisch aufgetragen und liefern die gezeichneten Kurven. Die 
Abszisse der Darstellungen ist im logarıthmischen Maßstab gewählt worden, um den aus- 
veedehnten Bereich von "X bzw. hR auf ein Zeichnungsblatt zu bringen. 

jede Abbildung noch je eine benachbarte Gröbersche Kurve eingezeichnet. 


0.0003 bis 0,01 des charakteristischen Parameters liegen. Die 


(sestrichelt ıst ın 








» * .. ee) 
Z a h ] enta fe ] » Temperatur der freien Oberfläche der Platte D, = . 

= 
hX 04 0.5 1 4 10 2) a0 100 200 >00 1000 2000 
at zu u ” Pan 2 m 
ys — (), 0003 0,999 0,996 0,9850 0.927. 0,833 0,705. 0,468 0.287. 0.157. 0,066 0,053 .0.017 
0.0010 0,997 0.983 0.965 0,872. 0,726 0,555 ..0,309. 0.171. 0,088 0,036 0.015 0,009 
0,0025 0.995 0,975. 0,947 0,809 0,615 0.441 0,211. 0,111. 0,057. 0,023. 0,011.0.006 
0,0050 0.993 0,963. 0.926 0,747 0522 0,336 0.154 0,079. 0,040 0,016 0,008 0,004 
0.0100 0,989 0.948 0,897 0.670. 0.428 0.256 0.111. 0.056 0.0285. 0,011. 0.006 0,003 
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’ , er , 
Zahlentafel 2. Öberflächentemperatur des Zylinders ®, S* 
w 
hR 0] 0. 1 4 10 >) DV 100 200 >00 1000 2000 
at : h = ae o En . f 
n: 0.0003 0.999 0,001. 0,981 0,926 0,831 0.703 0,465. 0,286 0,156 0,064 0,032 | 0,016 
>: 
0.0005 0.008 0.088 0.975. 0905 0.789. 0.643 0.394 0.230 0.122: 0.050 | 0,025 0,012 
0,0010 0.996 0,983 0.965 0,868 0.720 0.551: 0305 0,168 0.087 0,054 0.017. 0,009 


> 0945 0.815. 0611 0.421: 0.205. 0.107 0.054 | 0.022. 0.011 0,006 


0,0025 0.094 0.97 
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Abb. 4. Oberflächentemperatur der Kugel bei Me =1!9: % Ie Palma h R) . " 
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. i “ , a7 
Zahlentafel2. Oberflächentemperatur der Kugel ® pr 

5 
h R 0.1 0.5 1 4 10 >) 50 100 200 1000 
t - nn i . 

er 0.0003 0.999 0.991 0.981 0.925. 0.830. 0,701. 0,463 0.284. 0.154 0,031 
?2 

0.0005 0.998 0.988 0.975 0.905 | 0,757 0.638 0392. 0227 0.131 0.026 

0.0010 0.996 0.982 0.0964 0.868 0.717 0.545 0300. 0,164. 0.084 0.015 

0.0025 0.994 0971. 0.944 0.805. 0.603 0.413 0.197. 0.103 0.052 0.010 

0.0050 0.991 0.59 0.920 0.734 0.502 0.315 0.137. 0.071: 0.035 0.007 


Das Anwendungsgebiet der Kurven beschränkt sich auf die Voraussetzungen, unter 
lenen die Differentialgleichung des Wärmeaustausches im nicht-stationären Zustand den 
rtlichen Randbedingungen und den Anfangsbedingungen angepaßt wurde. Sie gelten also 
streng nur für die erwähnten Körper, wenn sie aus dem Zustand der örtlich und zeitlich 
konstanten Temperatur heraus abgekühlt bzw. erwärmt werden. Wenn die physikalischen 
körperkonstanten bekannt und als von der "Temperatur unabhängig angenommen werden 
können, ist es möglich mit Hilfe der Tabellen oder Tafeln rein rechnerisch, somit ohne Ver- 
suche, das Problem der Abkühlung (Erwärmung) zu lösen. Außerdem ist Voraussetzung, daß 
der Körper mit der konstanten Temperatur plötzlich in das Medium mit der niedrigeren 
(höheren) Temperatur gebracht wird und daß das Medium seine Temperatur unabhängig vom 
Vorgang beibehält. 

Kine andere Anwendungsmöglichkeit ist dadurch gegeben, daß man aus den experimentell 
ermittelten Temperaturen und aus den Schaubildern umgekehrt auf die Körperkonstanten 
schließt. Dies läßt sich, wie Gröber erwähnt, besonders für die Wärmeübergangszahl « 
durehführen. Die Bestimmung der Wärmeleitzahl / auf ähnliche Weise würde wohl auf 
erhebliche Schwierigkeiten stoßen, da 4 in beiden charakteristischen Größen, sowohl in 
” . E wie auch in AX “X vorkommt. 

X? cyX°® / 


Der unverhältnismäßig große Arbeitsaufwand, den die exakte Berechnung von Wärme- 
problemen im nicht-stationären Zustand erfordert und die stark bedingten und dadurch stark 
beschränkten Anwendungsmöglichkeiten auf immer nur ein einziges Problem, lassen es als 
erwünscht erscheinen, Näherungsmethoden anzugeben, die es ermöglichen, einen ver- 
änderliehen Wärmevorgang, wenn auch nicht exakt, so doch mit einer Genauigkeit zu 
beschreiben, die für praktische Zwecke ausreichend ist. Nachstehende Arbeiten verfolgen 
dieses Ziel. Von A. Haltmeier?) wird eine Annäherungsformel errechnet, die für die Aus- 
kühlungswärme von ebenen und zylindrischen Wänden gilt, wenn diese aus dem stationären 
Zustand heraus abgekühlt werden. Durch andere Überlerungen und unter Verwendung einer 
Fehlerfunktion kommt ©. Krischer?) auf eine ähnliche Annäherungsformel für den Wärme- 
verlust bei Auskühlung von ebenen und zylindrischen Wänden. 

Die Zahlentafel der Fehlerfunktion ist von W. Esser‘) angegeben, die mit der Arbeit 
von OÖ. Krischer zusammen erschien ‘). 


Zusammenfassung. Vorliegende Arbeit gibt für den Verlauf der Oberflächentemperatur 
bei der Abkühlung von Platte, Kugel und Zylinder, in Anlehnung an eine Arbeit von Gröber, 
weitere Schaubilder und Tabellen für den Bereich der kleinen Wärmeleitzahlen bei zugleich 
kleinen Werten der Zeit. Nach einer mathematischen Vereinfachung der Formeln für die 
Oberflächentemperatur, die Temperatur ın der Mitte und den Wärmeaustausch, die den Vor- 
gang streng wiedergeben, wird der Rechnungsgang für die Oberflächentemperatur beschrieben. 
Diese kann nur durch eine unendliche Reihe ausgedrückt werden, die gerade für das behandelte 
(sebiet sehr schlecht konvergiert und daher einen großen Arbeitsaufwand erfordert, der ja 
immer mit der streng mathematischen Lösung nicht stationärer Wärmeprobleme verbunden ist. 

224 

1), A. Haltmeier: Auskühlung ebener und zylindrischer Wände aus dem Beharrungszustand. Dissertation, 
Darmstadt 1926. (R. Oldenbourg, München.) 

5) 0. Kriseher: Die Auskühlung gerader und zylindrisch gekrümmter Wände aus dem stationären Zustand 
heraus, sowie die Anheizung derselben bei Zuführung einer konstanten Heizleistung. Dissertation, Darımstadt 1928. 
(Berlin 1930, J. Springer.) 


6) W. Esser: Wirtsehaftlichkeitsbereehnung isolierter Rohrleitungen und ihre wärmeteehnisechen Grundlagen. 
Dissertation, Darmstadt 1929. (Berlin 190, J. Springer.) 


7) Von den genannten 3 Abhandlungen erhielt ich erst Kenntnis, nachdem die oben mitgeteilten Reehnungen 
bereits abgeschlossen waren. 
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Die Berechnung kavitationssicherer Tragflügelprofile. 
Von ©. Schmieden in Danzig-Langfuhr 
Einleitung. 


n den letzten Jahren haben Kavitationserscheinungen bei technischen Problemen, besonders 

bei rasch umlaufenden Schiffspropellern, erhöhte Bedeutung gewonnen. Die üblichen Pro- 
file besitzen im Gebiet großen Auftriebes bzw. günstiger Gleitzahl als charakteristische Eigen- 
schaft eine Unterdruckspitze an der Flügelnase, die bei größeren Geschwindigkeiten dort 
leicht den Druck unter den Dampfdruck sinken läßt, so daß Kavitation einsetzt, deren un- 
eünstire Wirkungen sowohl auf den Wirkungsgrad wie vor allem auf das Material des Pro- 
pellers hinreichend bekannt sind. Bei Luftschrauben sind ebenfalls zu kleine Drücke uner- 
wünscht. da dann die Geschwindigkeiten leicht die Schallgeschwindigkeit erreichen, was 
wiederum eine Abnahme des Wirkungsgrades bedingt. Beseitigt man durch Verminderung 
des Anstellwinkels die Unterdruckspitze, so sinkt im allgemeinen der Auftrieb stark und die 
Gleitzahl verschlechtert sich erheblich. Es liegt somit folgende Fragestellung nahe: wie muß 
das Profil gestaltet werden, damit es bei einem praktisch zulässigen Anstellwinkel kavitations- 
sicher ist. aber trotzdem einen hinreichenden Auftrieb liefert? Da sich eine theoretische 
Untersuchung auf ideale Flüssigkeiten beschränken muß, muß die Frage der Gleitzahl dabei 
offenbleiben. 

Der praktisch zulässige Bereich der maximalen Geschwindigkeit am Profil liegt bei 
dünnen Profilen etwa zwischen 1,2r, und 14», (v„: Betrag der Anströmgeschwindigkeit) bei 
diekeren Profilen zwischen 15 v, und 18r,. Um einen möglichst großen Auftrieb zu erhalten, 
wollen wir die Profilform so zu bestimmen versuchen, daß der der zulässigen maximalen Ge.- 
schwindigkeit entsprechende Druck auf einem endlich langen Bogen der Kontur, der die 
Profilnase in sich enthält, erhalten bleibt. Als obere Grenze dieser Bogenlänge wird man 
etwa 40°/, der Flügeltiefe zulassen dürfen, da sich sonst der notwendige Druckanstieg auf 
einem zu kurzen Stück der Kontur vollziehen müßte, was starke Ablösungen und damit einen 
eroßen Widerstand des Profils ergeben würde. Die Aufgabe kommt somit darauf hinaus, zu 
einer qualitativ gegebenen Druckverteilung das zugehörige Profil zu bestimmen’). 

Durch die Forderung konstanten Druckes auf einem Stück der Kontur läßt sich die 
Aufgabe auf ein Problem der konformen Abbildung zurückführen, das dem in der Theorie 
der unstetigen Strömungen auftretenden sehr ähnlich ist. Im ersten Teil der vorliegenden 
Arbeit werden solche Profile untersucht, die teilweise von Geraden begrenzt sind; im zweiten 
Teile werden überall gekrümmte Konturen allgemein behandelt. Der dritte Teil bringt die 
Anwendung der allgemeinen Formeln auf zwei Klassen spezieller Profile, wobei die Rechnung 
bis zur expliziten Bestimmung von Profilkontur und Druckverteilung durchgeführt wird. 

Es ist mir eine angenehme Pflicht, an dieser Stelle Herrn E. Pohlhausen und 
Herrn G. Flügel für das Interesse, das sie meiner Arbeit entgegengebracht haben, zu danken, 
ebenso der Notzemeinschaft der deutschen Wissenschaft, durch deren Unterstützung die um- 
fanereiehen numerischen Reehnungen zum großen Teil erst ermöglicht wurden, schließlich 
Fräulein H. Just, Herrn P. Podlieh und Herrn G. Frischmuth für ihre Hilfe bei den 
numerischen Rechnungen. 


I. Bezeichnungen. Wir stellen zunächst die Bezeichnungen für die wichtigsten auftreten- 
den Größen zusammen. 


weg-tiy komplexes Potential, g Geschwindigkeitspotential, 
y Stromfunktion. 
z=ie+iy komplexe Koordinate in der Strömungsebene. 
dm we“ b; 
U. — iv, komplexe Geschwindigkeit, ®x, *, Komponenten 
(2 : 
nach den Achsen. 
dz u. * a Ai 
big, In iQ: y Betrag, © Riehtung der wahren Geschwindigkeit; 
dm ‘f 5. 
4, Geschwindigkeit auf der freien Stromlinie. 
v_ (Geschwindigkeit der ungestörten Strömung. 


1) Danziger Habilitationsschrift. 

?, Herr Weinig hat in einem Prager Vortrag, dessen Anszug (ZAMM (9) 1929. 8.507) mir erst nach Abschluß 
meiner theoretisehen Untersuehungen bekannt wurde, ein Verfahren skizziert, aus einem vorgelegten Hodographen das 
zugehörige Profil zu ermitteln. Dies Verfahren ist prinzipiell auch zur Berechnung kavitationssicherer Profile ver 
wendbar: jedoch dürfte für diese Profile unser Verfahren schneller zum Ziele führen, da es von vornherein darauf zu 
reschnitten ist. 
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F Zirkulation, o Dichte; + 
A=P,+iP, Auftrieb, P,, P;, Komponenten nach den Achsen. 
t,u,ı komplexe Koordinaten. 


W=a+tib; u =a,+tiß,: r,—r,e” Bildpunkte des Unendlich-Fernen der z-Ebene in 
den entsprechenden Abbildungsebenen. 
7 


0. BB: 9; 9,0} Winkel im Bozenmatß. 
“rs n 
#4, ri Konstante. 
k Krümmung. 
Rn, Real- bzw. Imaginärteil eines komplexen Ausdruckes. 


Konjugiert komplexe Größen sind durch UÜberstreichen gekennzeichnet. 


A. Die Grundlösung. 


2. Die Abbildung der [-Ebene auf die 1- bzw. u-Ebene. Die in diesem Teile zu unter- 
uchenden Profile seien von folgender Gestalt (Abb. I): zwei geradlinige Stücke, die den 
Winkel 5 miteinander bilden, und von denen das untere den Anstellwinkel « besitzt, seien 
ohne Kniek durch einen Kurvenbogen verbunden, so daß die Neigung der Tangente überall 
mit Ausnahme der Hinterkante eine stetige Funktion der Bogenlänge ist. Um den gekrümmten 
Teil der Kontur festzulegen, schreiben wir auf ihm Konstanz der Geschwindigkeit vor. 

In sehr großer Entfernung vom Tragflügel wird sich die Strömung wie eine Parallel- 
strömung verhalten, die von der Strömung eines Wirbels im Nullpunkt überlagert wird. Wir 
können somit die asymptotische Entwicklung des komplexen Potentials für sehr große z so- 
fort angeben; es wird: 





2 | 
wemv_ 2 ei I a (ı). 
e 2 
Machen wir die Substitution 2° _ , so verhält sieh in der Umgebung des Nullpunktes 
-ı 
der zEbene 
v. il 
WEB N en Eh AE- Year ME 
2, Pe / 
&-fbene 
Ö N 
z-Ebene BEE SEEN EN in 
C - 
p „ 
P] U u a: Pe 17 
mr, ( F_ 17% 
B A Vo, u u ku ORDENS NEESEENEUGEEENDESRSEGEERBEEESERE D 48 
[RArsoZ4 ” 4 
Abh. 1. | Abb. 2. 


Es ist also zur Erzeugung unserer Strömung im Unendliehen eine Doppelquelle vom 
Moment », und ein Wirbel mit der Zirkulation — /’" anzubringen. Bildet man den Außen- 
raum des Profiles in der z-Ebene konform auf eine andere Ebene ab, so sind nach Obigem in 
dem Bildpunkt von z= x eine Doppelquelle und ein Wirbel anzubringen, deren Moment bzw. 
/irkulation zwar nieht mehr dieselben sind wie in der z-Ebene, sich aber immer aus ihren 
lortizeen Werten bestimmen lassen. Da wir diese Größen nachher doch auf einem anderen 
Wege berechnen werden, gehen wir hier auf diese Transformation nicht näher eın. 

Von der gesuchten Strömung können wir noch folgende Aussagen machen: im allge- 
nneinen werden zwei Staupunkte vorhanden sein, von denen der eine in die Hinterkante 
zelegt werden muß, um glattes Abströmen zu erreichen, der andere bei positiven Anstell- 
winkeln auf der unteren Geraden liegt. Auf dem gekrümmten Teile der Kontur ist die Ge- 
chwindigkeit konstant; ihr Betrag ist allerdings zunächst unbekannt und seı mit q, bezeichnet; 
ıuf den Geraden ist die Riehtung der Geschwindigkeit gegeben. Es liegen also, wie oben 
chon erwähnt wurde, dieselben Randbedingungen wie in der Theorie der unstetigen Strö- 
nungen vor; deshalb wollen wir fernerhin das Stück der Kontur, auf dem die Geschwindig- 
‘eit konstant ist, als „freie Stromlinie* bezeichnen’). 

3) Aus der umfangreichen Literatur über unstetige Strömungen sei als für unsere Theorie wiehtigste Arbeit 
enannt: M. Brillouin. Ann. d. Chim. et Phys. (8) 23. 1911. Weitere Literaturangaben in dem Bericht von 


.Jaffe, ZAMM (I) 1921. S. 398. 


14 
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Da sich alle Randbedingungen auf die Geschwindigkeit beziehen, führt man zweckmäßie 
als unbekannte Funktion ein: 
S dz u; 2 dw 2 ' i ? j 
ing, =in—-+i0 9 — % = Rıiehtung der Geschwindigkeit . () 
dm 7 dz 
In der S-Ebene bildet sich die Kontur und ihr Außengebiet auf den in Abb. 2 gezeichneten 
Bereich ab. Damit die Abbildung eindeutig ist, muß offenbar q, größer als jede andere Ge 
schwindigkeit, speziell also auch größer als v, sein, da schon für ,=v, der unendlich fern 
Punkt der z-Ebene auf den Nullpunkt der Z-Ebene abgebildet würde; das widerspricht der 
Kindeutiekeitsforderung, da dieser Punkt schon das Bild eines Punktes der freien Stromlinie 
ist. In dem z = entsprechenden Punkte, also in Z,=In . 0, haben wir jetzt eine 
Doppelquelle und einen Wirbel anzubringen und das komplexe Potential der von ihnen in 
dem Bereich der Ö-Ebene erzeugten Strömung zu berechnen; dabei sind Moment und Zirkulation 
so zu wählen, daß beide Staupunkte dieser Strömung ins Unendliche fallen, wie es in Abb. 2 
angedeutet Ist. 
Um dieses komplexe Potential leicht angeben zu können, bilden wir den Bereich der 
-Ebene mittels der Schwarzschen Formel auf die obere t-Halbebene ab, und zwar so, daß 
sich die folgenden Punkte entsprechen (Abb. 3): 


A>t!=1: B>lt=»: (>1:=-0: D>l=x?: E>t-7°. 


Bei B und © haben wir Winkel von „, bei A und D Winkel von Null und bei # einen 


Winkel von ?rz. Also wird: ie 
dz e\ (t-4)dt 
dw vtt— )(t—»’)' 


« 


c=ing, (3). 


Es empfiehlt sich, um die Wurzel fortzuschaffen, die t-Ebene auf die «-Ebene abzubilden 
durch die Transformation (Abb. 4): 


vVi=u. 
t-fbene 
n u-fbene 
| P) 4 oo r2P) 
I\ Ki 
\ \\ IN 
N u En 
— \ in 
Od \#? 2 u* “ 0 YA Sr ” 
BED EI B C DE AB 
RAıso 73 [RA 18074 
Abb. 3. Abh. 4. 


Hierdurch wird die obere t-Halbebene auf den ersten Quadranten der «-Ebene ab- 
vebildet:; aus (3) erhält man jetzt nach einfacher Partialbruchzerlegung: 
fl Dar, DE u + x] 


„In - „In 


| C Sa). 
u 1 —-"u—1 »l—x) u 2 re 


Füru sv wird = i(a-a); beim Umlaufen des Punktes « = I nimmt Z um ir ab, beim 
Umlaufen von «= x nochmals um ©ß. Hieraus erhält man die Gleichungen: 


| j? CRe—a B 
C ai; z . ee 5: 
| — x? x | 2° 7 
und daraus endgültig: 
E u+1 BB, u+zx,. 
C = In + In un 5. Beer. | u © 1 
T | ı u X 


3. Bestimmung von w in der t- und u-Ebene. In der f- bzw. u-Ebene ist jetzt das kom- 
plexe Potential für die folgende Strömung zu berechnen: an der 2=» entsprechenden Stelle 
sind eine Doppelquelle und ein Wirbel anzubringen, die reelle Achse muß Stromlinie werden 
und schließlich müssen die Staupunkte in t=z° und t-1 fallen. 

Wir suchen zunächst die zweite und dritte Bedingung zu erfüllen; spiegeln wir Doppel- 
quelle und Wirbel an der reellen Achse, so wird diese Stromlinie, das komplexe Potential 
hat also die Form: 


‚fe iy ei; 1.3 f,\ jr 
‚ ( m 
[2 \/ t f , ) 7 In t tJ . . . . . . . . . D (»)). 


Ri 





se 


m 


SP 


Sti 
lie 


se] 


Ko 





en 


ab- 


3a). 


Im 


>h) 


om- 
elle 


den 


pel- 
tial 


(>). 


g nd 12, Heft 5 y . Fr . . . rn ” > 
Oktober 1939 Schmieden, Die Berechnung kavitationssicherer Tragflügelprofile 29] 





\* Pr n ar ° . r „T 
Eine längere Rechnung, auf deren Wiedergabe wir verziehten wollen, zeigt, daß y x , 


ein muß; um die Formeln zu vereinfachen, wollen wir daher von vornherein y=- , wählen. 


Setzen wir f,-a-+tib und die Konstante C, die nur den Maßstab beeinflußt, gleich 


ins®), so wird: 


Ih ı N ! dA ih ” 
w= 2 m n ze re De N er » 
t-a® + 2n"t-artib I 
n diesem Ausdruck sind noch drei freie Konstanten enthalten, von denen wir zwei dureh 
ER . dw ’ u R 
lie Forderung festlegen, daß für tx? und t=] dt 0 sein soll: aus (5a) erhält man: 
dm Ib(it a !' h 
— +2 
dt It — a)” + b°]? 2a(t — a)” +0? 


u .. dw . 
Für reelle t wird T reell: es verschwindet, wenn: 


Y Y % v 
t- a” IN, —2tt—-a) + I =0; 5 N 
also für: 
| . 
ta) r\4#] 1—(1,b)). 
dw, ) 
Andererseits soll lt für {=x° und !=1 verschwinden, so daß 
+ x? | 1 — x? l 
= a+-,; , vi—(I, b)? 
> r > | 
sein muß. Aus diesen Gleichungen folgt: 
/i | l-+x° 2\2 
a r h? ( 5) (6). 
a7 12 /} a Pf ja 
| en | > 
ls bleibt noch die letzte Bedingung zu erfüllen: da dem Punkte 1, 2=x entspricht, so 


u dw u ; 
muß für t=H, Fri werden, also der Imaginärteil von © verschwinden. Der z=x ent- 


sprechende Punkt der «u-Ebene sei u,—=a,+iP,, wobei 


a? » |} a®+b-+al: Bm - I] a@+b-—al. 
Somit folgt aus (4): | 
9°) 9 OR. 
I()=aretg; _ ie pe N Le, 
0 Po JT X Ag Po 
Der Bogen liegt hierbei, wie man sich leicht überlegt, in dem Bereich von — z bis 0; damit 
(7) eine Lösung hat, muß a,’—+ß,°<_ 1 sein, da nur dann das erste Glied kleiner als ei wird. 
Für die bisher noch unbekannte Geschwindigkeit q, erhält man: 
Rıö,) In Je . In \% he, 5 In Ce L e Pu (S). 
Ur 2 (a,--1”’+P, an —-H”’rP, 


4. Schließungsbedingung und Auitrieb. Aus e& erhält man durch Integration nach w die 
Strömung in der z-Ebene; wenn die Lösung einen physikalischen Sinn haben soll, muß sich 
(die der Kontur in der Ö-Ebene entsprechende Profilform schließen, d. h. das Konturintegral 


sein. Um diese Bedingung zu erfüllen, müssen noch zwei weitere Gleichungen zwischen den 
Konstanten bestehen. 


4) Streng genommen müßte man Ü negativ wählen; wir ziehen es aber vor, mit einer positiven Zirkulation zu 
"eehnen und den so entstehenden Vorzeichenfehler erst in den Sehlußformeln zu korrigieren. 


1y* 
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Das Konturinteeral können wir schreiben: 
$ ’ 5 dz dw ' 
GdZ . (IM 
i J dw du 


oder mit Rücksicht auf (4) und (Ha): 


. P da | | | E ur —ü E: 
Ö dz ® 2 Ar s u - 7 seztudu., 
, I dw \u’— u) (u u.) ar (u tt, tu“ a) 


Dabei ist das Inteeral über den in Abb. 5 gezeichneten Integrationsweg zu erstrecken; das 
Inteeral über den Viertelkreis mit unendlich großem Radius kann man hinzufügen, da es 
Null ereibt. Z/ und «, sind so bestimmt, daß die geschweifte Klammer an den Polen von 


von der Ordnung Eins verschwindet; die einzige Singularität des Integranden innerhalb 





dm 
4 a | "vo ; 
der Kontur lieet also bei « = w,. Hier ist rue Die ersten Glieder der Entwicklung 
SF 
dz . 
von u nach Potenzen von u - u, lauten: 
dn 
dz HM, be 5, (d (dz\\ 
u Hu u,)| u | 
dn Ü_ 10, rldu\dw)fu ug 
W > ist al 
egen "_ Ist aber: 
dz . 
Id (/dz\| I Id dz\) I /dE 
| ) | In(y, ) | 
'd „ \d Im, u Hy Ür \d [A d Ih EZ Hy ı : d u „H HU, 


Hiermit erhält man für das Residuum an der Stelle «,: 


MR i dE 64"7 
Abb. >. Ne) (n,)= | z 
u 2u,v, \(dulu=w 27, 





Dieses Residtuum muß verschwinden, so daß wir die Schließungsbedingung in der folgenden 
Gestalt erhalten: 


I (/dö\ I 


v.\dujlu—u 7 


(9). 


Da die linke Seite im allgemeinen komplex sein wird, erhalten wir zwei weitere 
Gleichungen, so daß sämtliche Konstanten festliegen, wenn eine von ihnen, etwa a oder p, 
vorgeschrieben wird; darüber hinaus besteht natürlich auch noch die Möglichkeit, daß es 


überhaupt keine physikalisch sinnvollen Lösungen unseres Gleichungssystems gibt. 


Den Auftrieb gibt uns das Integral: 


url as r+in 
l; 


Führen wir wieder « als Integrationsvariable ein, so erhalten wir): 

\ dım\? ’ pl“ w\ dz dw ' Be dw ' 

44 ; AU u dt, 
| dz | K dz dv du J dz du 
j 
dw dz 

Dieses Interral unterscheidet sich von dem eben untersuchten nur dadurch, daß / statt In 
dz ( 


im Integranden steht; berücksichtigt man die obigen Formeln, so erhält man für das Residuum: 


Nur N =. il’v, 
eh te,) 5%. Ün Fr FOR In 
Hieraus sowie aus (9) folet für den Auftrieb: 
pP. P.=z2iUta)afr,; KH=efo,; A=eb , «+ 4 3 IM 


‚ Weren der Integrationsriehtung in der «-Ebene vgl. Anm. 4, S. 291. 
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5. Gleichungsdiskussion. Mit der Abkürzung schreiben sich die beiden (9) ent- 
7 N 


;prechenden Gleichungen: 


a,t-a)+Pp,b , za,’ -a)+Pp,b ya’ + b’ 
(1 a?” +b  n (z""—a”’—+b? „fi „2\2 
| Da ar 1 Be 
ab—B,(1l—a), x a,b— B,(x"— a) 


= pe - = : 0 
1—-a”+Db n (z’- a” —+b’ 





er Ansatz b=x liefert z—=1, also ein unbrauchbares Profil, da beide Staupunkte in der 
IIinterkante zusammenfallen; für den andern Grenzfall sehr kleiner b wird in erster Näherung 
ı—= 2”, a,—%, P,=0 und damit aus (11), (7) und (8): 


HM = 





a3 nm; A P=V; y=I%, 


a das Bild des Unendlichen auf den Rand des Bereiches rückt, so erstreckt sich die Kontur 
ins Unendliche: man erhält demnach einen Parallelstreifen, dessen beide Geraden dureh die 
[reie Stromlinie verbunden werden. 


Für kleine endliche b erhält man endliehe Profile, die aber wegen des zu großen 4, 
praktisch unbrauchbar sind; einschließlich der Glieder zweiter Ordnung in b wird: 


re: Ze a 78 Ss do o Ms; 32 , 
= org b?; ‘eg mi 0 3 b; In rn In 3 G b° In b+5 b°. 


6. Berechnung der Profilnase.. Da es für einen Vergleich mit den später berechneten 
Profilformen von Interesse ist, soll hier doch noch die Form der Profilnase für den Anstell- 
winkel « ==9" berechnet werden: für diesen Wert wird: 


Fu Tat 
! I ) | )0 | l 
In Im 1; x u . 
is 2 7.300 7 
Die Gleichung des Profilkopfes erhält man durch Auswerten 
(les Integrales 
iu 
\ dz dw | 
’ AU. RÄ1307 
dw du ei CBr 
v Abh. 6. 


Auf die Wiedergabe dieser Rechnung verziehten wir und geben gleich das Resultat an: ver- 
nachlässigt man Glieder der Größenordnung 5’, so erhält man unter Verfügung über den 
Maßstabfaktor C: 


| ) | 3 " Ir 


E R . — ent > 
Si 3lroprd n=arcigv+] ,, te ro (12). 


Dabei liegt der Koordinatenanfang in dem Endpunkt der oberen Geraden, die &Achse zeigt 
nach vorn, die 7-Achse nach unten. Der Parameter » läuft von O bis ©: die Profilnase er- 


.. u. . .p ® » .. T » 
hält man für » 2, die größte Breite des Profiles beträgt ,. In dem entsprechenden Maß- 


stab beträgt die Länge der Geraden etwa 80, die untere Gerade ist um 1,5 länger als die 
obere, die Profilnase liegt nur um 0,1 vor dem Endpunkt der unteren Geraden (Abb. 6). 


B. Profile mit überall gekrümmter Kontur. 


7. Allgemeiner Ansatz für Ü. Die im ersten Abschnitt behandelten Konturen sollen im 
‘olgenden als Grundlösung bezeichnet werden und alle sich darauf beziehenden Größen mit 
lem Index „0* versehen werden. Um auch gekrümmte Konturen zu erfassen, gehen wir 
‚weekmäßig durch die Transformation 


ui SER SE Ve Or a HB EEE Bo Ein ER er Fe 
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von dem ersten (Juadranten der «-Ebene auf die obere Hälfte des Einheitskreises der r-Eben« 
über (Abb. 7)°); die imaginäre Achse der «-Ebene geht in das Stück der reellen Achs: 
zwischen I und +1, die reelle Achse der w-Ebene in den Teil des Einheitskreises deı 
r-Ebene oberhalb der reellen Achse über. Es wird also die freie Stromlinie auf das Stück 
der reellen Achse, die übrigen Teile der Kontur auf die obere Hälfte des Einheitskreises ab- 
eebildet, und zwar so, daß w=1 in r=i und w=x in r=ei%o transformiert wird, wobei 


ı-fbene l ei Yo ( 
A gi tg /,' 
' 1 + e' Yo “5 
\\ £ 
> ist. Setzen wir (13) in (4) ein, so folgt nach kurzer 
A) 5  wischenrecehnung: 
u N . | N 
. ‚T+i .T—eo-!M. , 
RE .+ .. Inei#o i.- ri —a).,: GH 
in ‚ t—eim 
‘ 
| a ;- - Wir verallgemeinern unseren Ansatz und setzen: 
B Ränoz 0 L 
Abb. 7. (== u, tiZl Be : . . ee . . F . . (15). 


Dabei wollen wir Z so wählen, daß die Zuordnung der Ränder von t- und ©- bzw. z-Ebene 
nieht gestört wird. Für reelle r muß der Realteil von £ verschwinden, also Z rein reell sein. 
Wir können somit nach dem Schwarzschen Spiegelungsprinzip die Funktion Z in die untere 
Hälfte des Einheitskreises fortsetzen, indem wir für konjugiert komplexe r Z konjugiert kom- 
plexe Werte erteilen. Die einzigen Singularitäten, die © haben kann, nämlich die den Stau- 
punkten entsprechenden, sind schon in £, enthalten, so daß Z in und auf dem Einheitskreis 
holomorph sein muß; es kann also durch eine im ganzen Bereich einschließlich des Randes 
konverzente Potenzreihe dargestellt werden: 


/= n (, TA s , ; : ö R ; r ö R s ; R R . ö . . (15a). 


Setzen wir r=rel”, so erhalten wir aus (15) und (15a) bei Zerlegung in Real- und 
Imarinärteil: 





SL 
Tr | I+2rsnd tr, 1] 1 —- 2rcos(#+9,)+ r? Fr 
In ‚In z ’ z " a,r"sinnd 
7 2 | 2rsınd + r° In | 2rcos(d —o,)+r? 
| ' 
. (16) 
.) y 
.«) ( ) | 
(-) asın-+t A Aa +m, 9, + (m, D,)-+ )y a,r"cosn ) 
2 n n Pr 
0 
mt 
‚ I+rsınd . Il -rsın d ‚,  rsın®+sing, „  rsind—sıng, 
tg, tom, =: go, : io, = 
rcos u) ! r cos ı) a rCoSU --.C0OSY, -.. rCOSU - COSgY, 


8. Ungleichungen für die a,. Die Konstanten a, sind in gewissen Grenzen willkürlich, 
so daß man sie benutzen kann, um Profile verschiedenster Form zu erhalten. Jedoch müssen 
sie eine Reihe von Ungleiehungen erfüllen, damit die Lösungen einen physikalischen Sinn 
haben: außer der sehon genannten Konvergenzbedingung müssen die a, unter allen Um- 
ständen so beschaffen sein, daß an keiner Stelle der Ebene die Geschwindigkeit q, über- 
schritten wird, wobei noch zu beachten ist, daß q, innerhalb der in der Einleitung genannten 
Schranken liegt. 


( . “2 r . .. 
Der Ausdruck In " muß also immer größer als Null sein; für Punkte nahe der Wand 
1 
(ir >1) liefert diese Bedingung speziell die Ungleichung: 


| Il + sınd | | cos (d-+Y,) ’ 
In ——— + . n a,sinnd >00 . 2... T). 
p | sın)  2n In | cos (d — Y,) u (17) 


| 


6) Die noehmalige Abbildung der Ö-Ebene auf die r-Ebene ist zuerst von Levi-Civita in die Theorie der un 
stetigen Strömungen eingeführt worden: man erreieht dadureh eine recht erhebliche Vereinfachung der Formeln. Für 
unsere Theorie wäre diese Abbildung aueh ohne diesen Vorteil noeh mit Nutzen einzuführen, da man durch sie eine 
Abbildung des Profiles auf einen endliehen Bereieh erhält, was für die späteren numerischen Rechnungen von großem 
Wert ist. 
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ir reelle r wird natürlich In Io -(), wie es sein muß; für Punkte mit kleinem Imaginärteil 
4 


-7 +ir”; 7’<1) müssen wir q9<T q, verlangen, oder mit anderen Worten: 


dd 
3 <o. 
dr A = 


urch Ausrechnen ergibt sich: 


F- 
.) ‘) ‘1 v 
2 2 sınyq 
—;,t+ en Ferm fi N BE. ee 
I+T? n1—2rcoso,tT? 


it den Spezialfällen ’): 


R l 1,9 
N N Ay: zZ etg “ Ehe En jr D)na,„<1} - tg E Ir -—-1I) .„ . (188) 
u ad un a and 


Steht hier das Gleichheitszeichen, so besitzt das zugehörige Profil eine freie Stromlinie mit 
endlicher Krümmung an beiden Ansatzstellen, in allen anderen Fällen wird die Krümmung 
an diesen Punkten unendlich. Der Beweis hierfür wird im Abschnitt 12 erbracht. 

Nennt man den dem Unendlichfernen der z-Ebene zugeordneten Punkt der r-Ebene 
',r,„e!%, so ergibt sich für das Verhältnis der Anströmgeschwindigkeit v, zu der Ge- 
schwindigkeit g,: 


Yy 1 1l+2r,snd, + r,/? | l—- 2r,cos(d, +o,)-+r,’ 5 
In ‚In re + iu ) ? z \, a,r,"sinnd, . (NM). 
v 2 l—2r,snd,+r? 2n 1-—2r,csd,—o)tr? 


1 
Der Wert dieses Ausdruckes muß innerhalb der angegebenen Schranken liegen. 


9. Hauptgleichungen zur Konstantenbestimmung. Schließlich müssen die a, noch so be- 
stimmt werden, daß die (7) und (9) entsprechenden Gleichungen erfüllt werden, d. h.: (a ist 
dabei mit in a, aufgenommen, da es jetzt keine physikalische Bedeutung mehr besitzt). 


an 
.) 2 1 .) R . Y 
Mr = /) = I, *% ‚ 
Jh: arc tg € „+ —aretg ; ——— t+r+ \ na,nr"cosnd =0 . . (2). 
B "-a— in) MW u A A a . 
0 
I | en 4 ne | EN 
IT U, \ .d Ulu— Zr ‚d T/t=t d Uu u) j { f j > 1 ; i u 
Aus (13) folgt: 
lI+iu i—-u dr 24 
T — u; 
| iu it+u du (tt + i)° 
Ferner wird: 
F- F. 
dZ R ” i— u\® 
"na, \ 7 a,| 
| 
Somit: 
F 
‚(dZ dr 2 Wi di — ut] 
| r - M a,| f 
dr duju=w (+ Ur +, 
PR ’ ad l 1. A 22. 
“ührt man noch die Abkürzung / u tg , ein, so erhält man nach einigen Um- 
) a 
Iormungen an Stelle von (11): 
et RICH BB 
2ya 4 b° % + P, tg .) 2 n A, Po ctg ») | 
FT. 
/ | ’ i— u, \ 1 
+y4b?+(1— x’) ya?’ + b? x | > n a,l. - 
(a ) MN U, rn \d + U, 

\ | (22) 
| AA m, x u ss). 
Kar ig .) Bu) + n A, etg .) B.) 

f SF 
| a K— u, \® 1 
+ y4b’+(1- °P - vya’—+b? y _ 2 nnd F ") )=0 
Yu,li+ N "\-+u, ; 





*), Für die unstetigen Strömungen wurden die (17) und (18) entsprechenden Ungleichungen von Brillouin]. «e, 
beeleitet. 
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Daher ist 


ee Bu 71 | ) N . » 
ya —(a,(1—-a)+ß,b+2b) -i(ka,b— B,(1l—a)+2a) . . . (22a). 


“(ti (++ pr] R u 


u 0 


10. Berechnung der Profilkontur. Mit den Gleichungen des letzten Abschnittes ist Z voll- 
ständie bestimmt. Um die Gestalt des zugehörigen Profiles in der 2-Ebene zu erhalten. müssen 
wir das Integral 


| \ dm 


dr (23) 
(fo h dr 
auswerten: die freie Stromlinie erhält man, wenn man r von I bıs +1 laufen läßt, die 
übrigen Teile der Kontur für r=e'’ mt 0=-J=— nr. 


"ür «die Durehführung dieser Rechnung müssen wir w von « auf r umreehnen: wir 


erhalten aus (5) mit y° , und Unterdrückung einer unwesentlichen additiven Konstanten: 
ZHIEZAZER | ı , „.d+:) If 13 | \ 
‘ . ns 1 RS — tm e 
% di) kr, nd) a) Er, der 
(24). 





2 (t-—-r,) (1 —rr,) 
Durch Differentiation folgt hieraus: 
dw .((+7r)% | | | | (1-+7,)° | | | | 
dı Ben ns) le —r)” A—ır))) 41— rd) —r® (A— rer) 


(24a). 





’ T[1—-2r,+7,/ l— 2rr, +7, | 

Für reelles r und für re!” wird dieser Ausdruck reell; da er aber in dieser Gestalt 
wenig brauchbar ist, schreiben wir ihn auf eine explizit reelle Form um. Eine längere 
/Zwischenreehnung ergibt: 


2 na |(t—r,)(l 7T,) (r—t,)(l—rr,) 





dm 65h Tl Hr” + r)| 4 
wur e- u. un: u 
dr Ita® -b’)(l+- Tr" +2all r-+(l ı)*]? 
3»). 
dır Ib I, sin29 [cos dl +z’) Al ”)] 
- a Zu vSsv—[n 
dı) |* (l + COS U)" + (COS ı) (1 —+ (l a))”]’ 


Der Zähler der rechten Seite der zweiten Gleichung verschwindet, wie es sein muß, für 


"==g, und 9, : außerdem für #0 und #=7, entsprechend dem Verschwinden der 


ersten Gleiehune für 7 -I und 7 l. 


Auf die weitere Durehführung der Integration wird im dritten Teile dieser Arbeit noch 
näher eingezangen. 


11. Berechnung des Auftriebs und der Druckverteilung bei anderen Anstellwinkeln. Die 
bisherigen Entwieklungen ergaben unter anderem die Abbildungsfunktion des Außenraumes 
unseres Profiles auf die obere Hälfte des Einheitskreises der r-Ebene bzw. auf die ent- 
sprechenden Gebiete der u#- und #Ebene. Wollen wir Auftrieb und Druckverteilung bei 
anderen Anstellwinkeln bereehnen, so müssen wir die Richtung der Anströmgeschwindigkeit 
in der z-Ebene um den entsprechenden Winkel ändern. Gehen wir durch unsere Abbildungs- 
funktion zu einer anderen Ebene, z. B. der f-Ebene, über, so müssen wir in dieser Ebene die 
Achse des der Parallelströmung entsprechenden Dipols drehen; weiter muß natürlich die 


Strömune in der F-Ebene so beschaffen sein, daß die reelle Achse Stromlinie und der der 


Hinterkante zugehörige Punkt ?°- 2° Staupunkt bleibt. Diese zweite Forderung bedingt eine 
2 e > ‚ ® er ® .. . T 
Anderung der Zirkulation 7. Da wir für die Anströmriehtung parallel der .-Achse y= ,, 


eesetzt hatten, machen wir den Ansatz (s. (B)): 


MR In na, ya Aue De AR ET 5. a hr "Sale EA 








co} 


Mt 


\\ 


oll- 


SelI 


+)) 


lie 


wii 


en: 


ı). 


lt 


l 
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d erhalten daraus: 





‚beosy (t -a)siny 1 ta ib) 
— r N - 
” - (t— a” +? 2x (Nt—a)+ib 
(26a). 
dm “ sın y KK a” —+b?]| 2 (t a)|ıt a) sıny beosy] 9 I I 
dt E [((t— a)’ + B’P "Bnltt— a)” +? 


'lierin ist #, als Bild des Unendliehen als gegeben zu betrachten. 


Damit {= x?” Staupunkt wird, muß folgende Gleichung erfüllt sein: 


E“ 
(2? a)(I'’b + sın y) I (2? a)lb eos y+b®(l bh sın y) 0: E- > 5 
Daraus folet: 
an 2(z” — a)b cos y i (x? — a)? — b? 
1,9 —- rn‘ siny, >; u, 
(x dt)“ + 5? (7° a’ —r h: 
| ’ h u 
Setzen wir noch , cetg , , so wird: 
ee —ıu 2 
EB ir Ye Od rl Re ee Eee ae 26 A 


Für den Betrag und die Richtung der Anströmgeschwindigkeit ergibt sich allgemein, 
wenn man mit », das komplexe Potential bei der Anströmriehtung Null bezeichnet: 


dw dw dt dw dt dw, 
u: dt dze dt dw, dz' 


I. 
dv 4 hd .. . * . > 
a "—r,, so folgt mit Rücksicht auf (Aa) und (26a): 
dz 
dm nv ur u 
» | „It a)sın y h cos y] U .„ (COS Y ı sın y): u. € ir. 
ai F 


ist also auch die Richtung der Anströmgeschwindigkeit. 


Den Auftrieb gibt uns das Konturintegral: 


: d a = Py+iP;. 


(Genau wie im ersten Teil rechnen wir das Integral auf « als Integrationsveränderliche um 


und erhalten: 
| 
i | /dw\. /du\ dz | /dw\°? - 
ee Pla) ae) ne 5 
| ( 


H 
In 
Auf dem Rande der Kontur lieet keine Singularität, da ın den Polen von dz vn 
di 


höherer Ordnune verschwindet. Die einzige Singularität liegt also bei «= «,; an dieser 
Ä 0 


‘ . ( iu a y . 
Stelle besitzt nach (9) foleende Laurent- Entwicklung: 
du 


lı 
dz ! | p ! 
‚+ Pu —u,). 
du 2u,v, (u — u,) ’ 
Weiter ist: 
d Th fZ iy ei T nr | l | | \ 
.); I 
-ın 2 212 2 = 212 ) 2 p, 2 2 . 
du lu u, WW —ü) 7 HM, u TE 


Somit ergibt sich für den Integranden: 


„[dw | 2 
>Din,r, (m “|, = + — uw” P,luw— u). 


ae dw i og RR 
“ür die Größe (m | erhält man die folgende Entwicklung: 


du 
dw IE m + 
v—ı / Hu - uw) P,u u, 
| DT Tr 2u,u u, / | dei | 
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und daraus durch (Quadrieren das Residuum: 
1 LE mn Eu 2 ;. vr > iu 
Neil lu, — sv, E; a7. 


Also wird der Auftrieb: 


.) 


> i- 2 . on .f ” a = Un . ) ar ON 
P,+iP,=riRdlu)=v, !e'!= „ inlatp)e Be 
ö j 2 dw 
Zur Bereehnung der Druckverteilunge an der Kontur brauchen wir die Werte von j, am 
dz 


Profil; es ıst 
dw 
dw dw, dr 
dz dz dw, 
dı 
dm . 
Für 7 erhält man nach längerer Zwischenreehnung einen ganz ähnlich wie (25) gebauten 
Ausdruck: 


dv dm, 165° sın y(1 —r’) (l- MM +2all - PP +Hr (2a rl -+r)" 29 
CuoSY T 15) > 1) " f 0 | 4 N) 15) 1) . öe |. 
dr dı (’- a?’ +b [1—-r’+2al—P’+(l+T'(a’+b°)]’ 
Somit wird: 
d n dın, : 7 | r | „? R) (1 } r)®|\ .) 
z- dz \e0s7+ sın yetg > 1 ”+| > |» ir r er 3 
Für rei” erhält man hieraus: 
dw dw,f u | } 1— x?\?1-+cosV|] 
“ DONE 25 SI ar g+Tır | I? 2 | B 
dz "FU ah En; ur b > 2 cos» (1) 


a TER „T . am N. 
Die Klammer wird für 9 ‚ unendlich; da aber dort ' von derselben Größenordnung ver- 


( 
2 dz 

i . ar s ö o : ’ ” 
schwindet, bleibt F endlich. Die Lage des Jetzigen vorderen Staupunktes erhält'man durch 

( % 

Nullsetzen der geschweiften Klammer: 
b sin y 
vos ı, re 


u 
(1 COS u)(eos y + sin y etg > 


bsiny 


Mit „>0 geht 9, gegen 5, wie es sein muß. Da beine recht kleine Zahl von der Größen- 


ordnung 0,05 ist, so kann man wieder nach b entwickeln und erhält bei Vernachlässigung von 
Termen, proportional b*: 


cos u — | 2b; cg>—_» 


Damit wird (32): 


vos, 


m\ 

I 2bletey-+ectg 5) 

Setzt man yxb, so wird für y_ 0 cos #, x — 1, der Staupunkt rückt also für negative Anstell- 
u 


winkel rasch gegen 9,7: man rechnet leicht aus, daß dieser Fall für ctg y ctg 5 


eintritt. Für dieses y erhält man für die Geschwindigkeitsverteilung: 


COS ”—| — (OS 5) 


dm dm, : | (1 -+ 7 dw,y, 12 | z»??1-+ cos) 
- .- z i 
dz / dz ) l + 7° / dz \ y COS ı) 


Der Staupunkt liegt in der Ansatzstelle der freien Stromlinie, die Geschwindigkeit auf der 


freien Stromlinie ist anfangs sehr klein und bleibt immer kleiner als q,, so daß der Auftrieb 
stark verkleinert wird. Für y a4 wird er übrigens bereits Null. 


‘ 


Be ® 





IreW 
lech 


——— 


28), 


am 


ten 


A), 


0). 


I- 


NM 
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a h ji h u 
Setzen wir andererseits im Bereich der positiven Anstellwinkel speziell etgy = ctg 5, 


‚wird cos 9,» und: 


> 














aw cos dw, ») ip: + i a Schul, 4 cos an. ) | + | er ui 
= COS Y ve‘ / f E sv Ye > , “ 
dz ' dz 2 1+7 dz 2) cos ı) 
dm 


ür reelles r ıst (0; also erhöht sich die Geschwindigkeit an der Flügelnase in diesem 


dz 
"alle ungefähr auf das Doppelte und sinkt dann für >1 nahezu wieder auf q, ab. Mit 
vachsendem Anstellwinkel nimmt also die Geschwindigkeit an der Nase sehr stark zu, es 
‘tellt sich qualitativ die auch bei den üblichen Profilen bekannte Erscheinung des großen 
Unterdruckes an der Nase ein. 


12. Diskussion der Krümmung‘). Wir haben noch den Beweis für die Bemerkung des 
\bsehnitt 8 bezüglich der Krümmung der freien Stromlinie an den Ablösungsstellen nach- 
zutragen. Da wir durch die Untersuchung der Krümmung auf interessante Zusammenhänge 
„wischen der geometrischen Beschaffenheit der Kontur und den a, geführt werden, sollen die 
entsprechenden Formeln in voller Allgemeinheit abgeleitet werden. 


Für die freie Stromlinie ist In Io und y eleich Null, so daß 
q 85 


wo; (6=1®, 


dabei ist © die Richtung der Geschwindigkeit, also auch des Tangentenvektors der freien 
Stromlinie. Da die Geschwindigkeit dem Betrage nach gleich q, ist, wird, wenn s die Bogen- 
länge bedeutet: 


dw 0« ’ OR, 
ö a. do; P=G% Fr c. 


dz Os Os 


Für die Krümmung der freien Stromlinie erhält man daher: 


d I de do dw 


TR —r I rn ae ee a 
ds I dag lo dt ' dı vn, 
do GL. 
Da nun “ Faludagi 5 >= ist, ferner 
dc ( 2 8 sing, dzı 
) .- z 
dr I+®" nl 2rc0sy, + 1? dr)’ 
so geht (33) über in: 
n | 
| 2 2 sınq, ’ 
k=—q, > z .. ie nat" \ - dw 20.2... (34). 
L +7? n 1—2rcosp,+T” << 
) dı 
” dw a 
Da für = +1 er verschwindet, so ist im allgemeinen an beiden Ablösungsstellen die 


u Bas dT .. 
Krümmung der freien Stromlinie unendlich, es sei denn, daß auch Ar für = + I von der- 


selben Ordnung verschwindet. Damit das der Fall ist, müssen die beiden Gleichungen erfüllt sein: 


IZ | \ IZ I“ 
=+1 > —1- ete 1°: 7 I > ==] tg 
dr wu dr „3 


Diese beiden Gleichungen, die jetzt auch die Bedeutung des Gleichheitszeichen in (18a) er- 
kennen lassen, werden durch den folgenden Ansatz für Z befriedigt: 


j | | | un . nn . 
Z=tr(l+ „, cosec q o a. er etgy,+ n but" r n +98 e a ee 
0 


*) Die entsprechende Untersuchung für unstetige Strömungen wurde von Brillouin durehgeführt (I. e.). 
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Der Zusammenhang zwischen den alten und den neuen Konstanten ist also der folgende: 


| | 
A I+ COSEC -B.: 7 ete « -h 
n ”. : 79,0 ro s 


.)- 
(A) 


« 
Dn 


| n — 2 
(dl. h 5 DB: (d,, h h > 
, ; +) n r 





Während im allgemeinen an den Ablösungsstellen nur die Tangente stetig bleibt, bleibt füı 
diesen speziellen Ansatz auch noch die Krümmung stetig, was in einem glatteren Verlauf de: 
zugehörigen Profils zum Ausdruck kommen wird. 
Ebenso wie für die freie Stromlinie bereehnen wir jetzt die Krümmung für die übrigen 
Teile der Kontur: wieder ist: 
Re) d@Q) ds 
erse ds dd» ‘d ı) 


’ 


| | ( 
da jetzt 7 ei” zu setzen ist. Da y 0 ist, so wird mit der Abkürzung In ER 


‘1 
7 ds e” dg 
do e -.; ( " ( 
Os dV dd 
und somit 
d (-) dq 
k hey we 26 
do‘ da "dd vun. 
Unter Berücksiehtizung von (16) erhält man aus (36): 
1 .i . | 
cos U) | cos(V -« N =Nnaysinn‘ y _ 
I 0 | e | \, na,snnd dw ... (N. 
| — sın ”\l cost +9 A 
l dd 


Die Krümmung der Kontur wird immer unendlich für 9 = g,, unter Umständen aber auch 
noch für J=0 und J=7. Sie bleibt hier immer endlich, wenn statt der a, nach (35a) die 
b„ eingeführt werden. In diesem Falle ergeben dann (34) und (37) denselben Wert der 
Krümmung an den Ablösungsstellen, so daß die Stetigkeit der Krümmung eine Folge der 
wesentlich schwächeren Forderung nach Endlicehkeit der Krümmung an den Ablösungs- 
stellen ist. 

Aus dem Verschwinden des Krümmungsradius folgt übrigens keineswegs, daß Profile 
mit dieser Singularıtät unbrauchbar sind: im Verlauf der Kontur konımt sie vielmehr kaum 
zum Ausdruck. Nur die Druckverteilungskurve weist an diesen Punkten einen Knick auf, da 
orad p dort unendlich werden muß: ist nämlich vd die vektorielle Geschwindigkeit, so besteht 
zwischen ihr, dem Druckgradienten und der Krümmung der folgende Zusammenhang: 


| 
h: gi vgradp]. 


C. Berechnung spezieller Profile. 


13. Allgemeine Bemerkungen. In diesem Teile werden die bisher entwickelten Formeln 
auf die Berechnung spezieller Profile angewandt; dabei beschränken wir uns auf die Unter- 
suchung der beiden einfachsten Klassen von Profilen, nämlich den Fall a,, a,=-0, a„=0 und 
den Fall b,, b,==0, b,==0. Im ersten Falle erhalten wir Profile mit unendlicher Krümmung 
an den Ablösungsstellen, im zweiten Falle solche mit endlicher Krümmung an diesen Punkten. 
Es genügt, sich auf diese beiden Fälle zu beschränken, da diese Profile schon von praktisch 
brauehbarer Form sind und die Hinzunahme weiterer Glieder der Reihe entweder nicht viel 


bringt bei kleinen a, —, oder aber bei größeren a,„ — die Profilformen verschlechtert 
bzw. vollkommen unbrauchbar macht. Sind nämlich die folgenden Koeffizienten ihrer Größe 
do 


nach vereleiehbar mit dem ersten, so wechselt auf der Unterseite des Profils mehrmals 


Av 
sein Vorzeichen, und zwar treten diese Wendepunkte, wie durchgereehnete Beispiele zeigen, 


verade in dem Bereich 9 x , auf, also in der Umgebung des vorderen Staupunktes; solche 


Beulen“ sind natürlich unter allen Umständen zu vermeiden. Schon kleine Werte von 
do 


19° ohne daß Nullstellen innerhalb des Intervalles auftreten, bedingen eine ungünstige Aus- 
(ta 


ae 
is er 
y „s 4 
73 
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'Idung des Profilkopfes. An den Stellen J=n 
(-) dw 
dageeen nichts, da dort auch 19 verschwindet. 
I ’ \ (tı 
. das man mehr mitnimmt,. die erforderliche Rechen- 
heit so stark. daß eine evtl. zu erwartende kleine Verbesserung der mit wenigen a, be- 


I 


Außerdem wächst mit jedem «, 
schneten Profile den dazu nötigen Arbeitsaufwand nicht mehr lohnt. 

14. Der Fall a, -=-0; zwei Spezialfälle.. Wir untersuchen zunächst den einfachsten Fall, 
imlich nur a,==0 und a,=+0; um Anhaltspunkte über die Größenordnung der verschiedenen 
'ofilkonstanten zu gewinnen, diskutieren wir die beiden Spezialfälle b=x und b 0, für 
22) leicht in geschlossener Form angeben lassen. 


‘ie sieh die Lösungen der Gleichungen 


Für b=x ergibt sich: 
) r rc | r 24 ‘ | ex 
30: m, > 7,0080, Ber er r,sınd), I+ydc + 
. 44/9... 
| | m L—x N | | | I 1+2y2x-+x 
v2» 1+y2#2+x) (u, + 0) V=-xz(1+y2x-+x)’ 


\ ° 
Ylu,te,+ ©) 
‚ein, so erhält man nach Kürzen gemeinsamer Faktoren: 


Setzt man alles ın (22 
2y2x l+x | | 2 I+2V2r +x 
+ d,, 
1—y2r+x 1-yV2r+tr\ nn I+yV2r+x 
| i | Ad, 
I—-y2x+x n I+y22-+x 
Kliminiert man hieraus a,, so folgt: 
li+x2=y2x(äön—]1). 
(uadrieren liefert für # die quadratische Gleichung: 
„” —Sın!(n +10, 
I sein muß, erhält man daraus: 
1”? —1, 


Da 2 
z—=4n(n I) ylön’in 


Kinsetzen dieses Wertes in die obigen Gleichungen gibt für «a, 


| 


BE, 
Aus Gl. (19) und (20) folet: 
Yo | / 1\ | | 
in pr „(+ mi n In’ 
N 2% l—x 
(, 1 )aretg! e 5 
N, I-z I+y2x-+x% 
Unter dieser 


Für praktisch in Frage kommende Profile muß n=5 oder größer sein. 
- eine kleine Zahl, so daß man überall Reihenentwicklungen veı 


Voraussetzung ist also - 
wenden kann. Vernachlässigt man = und höhere Potenzen, so erhält man: 
| Ni | | I, 41 un 43 
) 0) ) ie ' 2. da = 65} = % 
Sn’ n Sn’ Sm" . In?” 24m 
3 | | we 3 | 
> ( = zum; 
I tin? 12n’) * 


In r =; 

"_ In’ 12n 
Damit sind alle Konstanten bestimmt; zu Jedem » gehört ein einziges Profil, das eine 
ganz bestimmte Orientierung gegen die Anströmrichtung besitzt, die durch die zweite Gl. (16) 


Wir führen die folgenden Winkelbezeichnungen ein: 


festgelegt ist. 


Vo 
ee 
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„Anstellwinkel* «= Riehtung der Tangente an der Unterseite der Hinterkante. 


Yy 
J r m N‘ 

N, (— 1)" a,, | ») Riehtung der Tangente am vorderen Staupunkt. 
A _. 

u 

N I)" «a, (7) Rıiehtung der Tangente an der Ansatzstelle der freien Strom 
DZ 4 D 

v linie an der Unterseite. 

N, (dy 4, (0) Del., aber an der Oberseite. 
ne Z 


Die Winkel sollen dabei so gezählt werden, wie es in der Abb. S angegeben ist. 


9,60-& 








Abb. 8. 


In unserem Fall erhalten wir für n >: 


| 8 “ | | 
at la 2? Un 
kaz | 
NE OT 


a liegt also in einem zulässigen Bereich, Jedoch werden die beiden anderen Forderungen, die 
wir für ein Kavitationssicheres Profil gestellt hatten, nicht erfüllt: die Geschwindigkeit auf 
der freien Stromlinie ist nur um 3 vH. größer als die Anströmgeschwindigkeit, der Auftrieb 
wird also viel zu klein, ferner erfolgt, da ©, (a) >, ist, die Ablösung der freien Stromlinie 


erst hinter der Nase. 


Für «en zweiten Spezialfall erhalten wir: 


b=0: a=x°: De “ U een, n 9Q 

und 

\ l — x”\° 

"P— r zw ..) 0.5 1,COSgy, 

4x 
TR l+x 2x l+x 21-2) +°) 
In In a: ‚in Ir 

V, | l+x° | H (1 x”), 

Da 7, e’vo ist, fallen die Bildpunkte der Hinterkante und des unendlich fernen Punktes der 


»-Ebene zusammen, so daß wir ein sich ins Unendliche erstreckendes Profil (Halbkörper) er- 
halten. Für die oben eingeführten Winkel ergibt sich: 


a . een, a,= a, (l+cosg,): 


7 
a»: | z d d,COSg 


IWW) a a a,(l -cosg,). 


Sämtliche Konstante sind als Funktion von x gegeben; für ein festes x erhält man also 
wieder einen ganz bestimmten Halbkörper. Wir wollen die Abhängigkeit von q, und a, von 
x näher untersuchen: für z=0 wird a, q, 1, ein Fall, der schon in dem oben behandelten 


Tan 


a 


... 


Kenne 


... 
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» für b 0 sich ergibt. Läßt man z zunehmen, so nimmt a, ab, bleibt aber zunächst 

ositiv. so daß die Unterseite konvex wird; q, nimmt dabei monoton zu. Für den ausge- 
| : ’ et 

siehneten Fall = wird a,=0; = 3v,, Unter- und Oberseite arten in Gerade aus (vgl. 
\bschn. 5). Wird x noch größer, so wird a, negativ, die Unterseite wird konkav, 9, >®3r,: Da, 
vje wir nachher zeigen werden, bei festem x mit wachsendem 5b q, zunimmt, so können wir 
‚ııs den beiden bisher behandelten Spezialfällen folgende Schlüsse ziehen, die natürlich streng 
ıır für die Klasse a,=+0 gelten. Angenähert werden sie aber nach den Bemerkungen des 
\bsehn. 13 auch im allgemeinen Fall gelten. 


Profile mit konkaver Unterseite besitzen ein q,, das sicher größer als die dreifache 
\nströmgeschwindigkeit ist, sind also nicht kavitationssicher. Profile mit a, = 1, also xx 0 
ind zwar kavitationssicher, liefern aber einen zukleinen Auftrieb, da q, > v, ist. Brauchbare 
Profile werden sich also nur ergeben, wenn x Me und b)< x ist, Da außerdem die Un- 
oleichungen (17) und (18) immer für a, 1 erfüllt sind, erscheint folgendes Verfahren als 
;weckmäßig: 


Man nimmt ein x» [ und ein b< x fest an und löst die Gl. (22), die in - und «a, 


linear sind, nach diesen Größen auf. Aus (19) erhält man q,, aus den übrigen Gleichungen 


M a ff ge (> i JT \ . . 
«, und die kennzeichnenden Winkelgrößen a; 9, (m); 4, | z ): J,(0). Das Ergebnis dieser 


Rechnungen enthalten die folgenden Tabellen, in die zum Vergleich auch die entsprechenden 
Zahlen für die Grenzfälle b=x und b=-0 aufgenommen worden sind. Für die beiden 


® 24 « .f JT . . de . . . 
wichtigsten Größen 4, und gr ist die Abhängigkeit von x und b noch in den Abb. 9 und 10 
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h 0,00 0,02 0,05 0.10 0,20 10,00 0,02 0,05 0,10 0,25 [0,00 0,03 0.05 0,08 0,30 


si” 
N IL 60” 19° 16° 4856’ 125°] 0° 5°0’ 14° 55° 37° 20° 130°] 0° 6°20° 11° 50’ 21036’ > 135° 
a ()" .00 32’ .90 36’ -6° 56’ .17° 0" 01’ -0057' 3057’ n.-170 ()9 v0 38’ 033’ .0° 11’ u .17° 


7,10) 1 0,055 0.0681 0,1093 0,1888 0,355 | 0,073 0,0809 0.1047 0,1656 0,407 | 0,090 0.0926 0,0993 0,1195 0,453 
JT 

2 | > 0,649 0.6977 0.7209 0,7074 0,645 10,567 0,6102 0,6435 0,6556 0,593 10.478 0.5344 0,5604 0.5843 0,549 

6er) | 1,353 1,4635 1,5511 1.6036 1,645 | 1,209 1,3013 1,3917 1,4768 1,593 | 1,052 1,1614 1.2201 1.2881 1,549 


1.144 1,234 1.355 1,557 ..19 1,232 1.335 1.437 1,628 2,1 1.354 1.481 1.5553 1,660 23 
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Aus diesen Tabellen und Kurven ersieht man, daß der Kantenwinkel hei festem 
n 


mit wachsendem 5b monoton zunimmt, während er bei festem 5b mit wachsendem 2 monotoı 
abnimmt. Die Geschwindigkeit q, ist eine monoton wachsende Funktion von b und 


Da also mit wachsendem #2 4, zunimmt, z ddareren abnımmt, kann man in gewissen 


Grenzen für ein Profil den Kantenwinkel und gleichzeitig die maxımale Geschwindiekei 
am Profil vorschreiben, ein Resultat, das für die praktische Anwendungsmöglichkeit unsere: 


Theorie von großer Bedeutung Ist. 
Schließlich sei noeh bemerkt, daß bei sämtlichen brauchbaren Profilen 9, (m) x, ist, 


so daß das Druckminimum in unmittelbarer Nähe der Profilnase auftritt. 


15. Numerische Berechnung von Profilkontur, Druck verteilung und Auftrieb. Für die Durch 
führung numerischer Rechnungen ist es zweekmäßig, für Zare tg 2=g, einen solehen Wert 
zu wählen, daß das Bild der Hinterkante in der --Ebene in einen brauchbaren Teilpunkt des 
(uadranten fällt; für die im folgenden berechneten Profile ist 9, = 18°. Die erforderlichen 
Integrationen wurden numerisch ausgeführt, da eine Berechnung mittels Reihenentwieklung 
sehr umständlich wäre, in der Umgebung der Hinterkante sogar praktisch unmöglich ist. 
Dieser Teil der Kontur bietet der Reehnung überhaupt die meisten Schwierigkeiten, da sich 
dort der Inteerand nieht nur sehr rasch ändert, sondern auch zahlenmäßiz sehr groß wird. 
Wie man nämlich aus (25) und den Tabellen des vorigen Abschnitts erkennt, verschwindet 

dm 
zwar der Zähler von 1 für 9 q.,, aber der Nenner wird ebenfalls sehr klein und zwar 
sehon in unmittelbarer Nachbarschaft der Nullstelle des Zählers um zwei bis drei Zehner- 
potenzen kleiner als der Zähler. Daraus erhellt schon, daß eine Reihenentwieklung des Inte- 
eranden an dieser Stelle kaum durehführbar ist. Bei der numerischen Integration kann man 
dageren diese Schwierigkeit durch eine hinreichend feine Teilung des Integrationsbereiches 
überwinden. 


Der Inteerand wurde daher bei sämtlichen Profilen an folgenden Punkten berechnet: 


auf der reellen Achse von =  |Ibısr 0,1 ın Schritten von O,l, von r Vlbısz H11 
. . j : rn 5 T . ’ . 
in Sehritten von 0,05; auf dem Einheitskreis von J nz bıs U » +9" In Sehritten von 9, 
„T 0 ’ : : . aan f TE ’ : 
von „) = +9 bis’ 0’ 30’ in Sehritten von 4° 30’ und von II’ 30° bıs t’ 30’ ın 


Schritten von 2" 15’. Der Integrand wurde also teilweise noch um zwei Teilpunkte über das 
: | - dz dz 

Integrationsintervall hinaus berechnet, um an Stellen großer Anderung von bzw. 19 noch 
{ (T dı 


die höheren Differenzen bis zur sechsten einschließlich bestimmen zu können. In der Um- 
gebung von 7 I(9 a) ist das nieht nötig, da der Integrand dort so klein wird und sich 
so wenig ändert, daß die Genauigkeit der Simpsonschen Regel ausreicht. Die notwendigen 
Rechnungen wurden mit einer Stellenzahl begonnen, die der Genauigkeit der mir zur Ver- 
fügung stehenden Tafeln entsprach”), d. bh. im allgemeinen sechsstellig, teilweise siebenstellig. 
Diese Anfangsgzenauiekeit ist notwendig, wenn man im Endresultat für jeden Funktionswert 
mindestens vier Ziffern erhalten will, da mehrmals die Differenzen nahezu gleicher Zahlen in 
die Reehnung eingehen. Zur Berechnung der Integrale wurde die Integrationsformel I ver- 
wandt'®), bei der man den Integralwert nur in jedem zweiten Teilpunkt erhält; war der Ab- 
stand von zwei so bereehneten Punkten der Kontur zu groß, so wurde noch ein Zwischen: 
punkt mit der Integrationsformel II") eingeschaltet. Die sechsten Differenzen dienten im 
wesentlichen zur Fehlerschätzune: in den Tabellen sind die Koortdinaten nur mit so viel Stellen 
angegeben, daß die sechsten Differenzen die letzte Dezimale maximal um eine Einheit ver- 
fälschen können. Eine sehr wirksame Fehlersehätzung gibt noch die Schließungsbedingung: 
da die Integrationen von der Nase aus sowohl entlang der Ober- wie der Unterseite des 
Protils durchgeführt wurden, kann man aus der Abweichnng der auf beiden Wegen berechneten 
Koordinaten der Hinterkante eine Abschätzung des Gesamtfehlers erhalten; dieser Gesamt- 
fehler beträgt bei allen berechneten Profilen ungefähr O,1 "so. 


»„, Havashi, Sieben- und mehrstellige Tafeln der Kreis- und Hyperbelfunktionen. Peters, Sechsstellige 
Tafeln der trigonometrischen Funktionen. 
10, Runze-Königr, Numerisches Rechnen. NS. 239. 


2) 8, a. OB, 3, 34. 








) s ‚f , ° . . . 3 . nm.; . kon x .; \ ‘ . 
\ Bank IR Ei e Sehmieden. Die Bereehnung kavitationssicherer Tragtlügelprofi« 30H 
| Oktober 1997 ERS, . u. Ne EBERLE ARE _— - 






























































1 uns 1 
| | 
“7 | 
JE 
„4 / 
J | / 
/ 
| R 
f > 
, l 
" 
nr} 
“ n 
nn = un 
1} 
| = = 
| = | ; | = 
| | e N N | 
Pr | > 2 
— N - 
| | er I 2; | — 
. | _ I Sm | u 
1 } o& . 
un I LI || : © a 
| ap“ I ’ > .. 
WW a . I = 
ze nn 
. l 2 m 
N | = Ss. r- 
Io | per n ER - 
5 ‘ u 
f Ü 
ft | Sn | | 
| S 4 
> 
| x 
Du 
of | | 
| 
| os 
u | || 
| NS Il 
m \ 
N = 
. Ä /I\\\e -|\ 
“4 NSS \ : 
SE — la _ a: 
I = 
SQ z 
4. 
LU. 
‚| 
() 
j k 
i 
- T 
2 
N 
LS | 
1} 
h 
I 
l- | 
h 
| 
N | | 
. = | Zi 
l'- = | - 
. | | 2 
r - 1) a 
_ _ | — 
. = Epei- 
t F- | go 
5 we 
I’- g 0% be 
n. gen S 
nn 
)- N R . R 
l z ” 
ll | | 
| 
N | 
1 | | 
E | 
N * | | 
SS 
2 ! 








ep 
[3 
1b 











5 ” er 5 
E | ri f r 
12 E7 j > 
be) h . 1 - 
| I 5 | -1 <c | 
Zen: > - An w sinsnnibeue 


r( 


u 








Ztsehr. f. angew, 


306 Math. und Mech 


Schmieden, Die Berechnung kavitationssicherer Tragtlügelprofile 





Die in den Tabellen angegebenen Koordinaten sind auf einen für die numerische 
kechnung bequemen Maßstab bezogen es wurde ın (5) € Ihr eesetzt (vel. (25)) —. der 
Ye 


1 

Koordinatenanfange liegt in der vorderen Ablösunesstelle der freien Stromlinie, die x-Achse., 

die mit der Anströmriehtung übereinstimmt, zeigt nach rechts, die y-Achse nach oben; in den 

Zeiehnuneen wurden alle Profile auf die gleiche Länge reduziert. In den Abb. 11 bis 14 sind 

7 s 2 - . . a 
aufgetragen, und zwar die zur Saugseite gehörigen ober- 


halb und die zur Druekseite eehörieen unterhalb der Profilkontur. 


außerdem noch die Werte von | 


) 


- 


I s x Y ° > 
Mn : auf diese Gerade als Nullinie 


Die strieh-punktierte Gerade entspricht dem Werte | 
‘fo / 


. ‚ . { 
bezogen, stellt dann die Kurve der Werte von | L 
\Yfo 
Druckverteilung im üblichen Sinne dar, und zwar bedeutet ein Verlauf der Kurve oberhalb 
der strich-punktierten Geraden Unterdruck, unterhalb dagegen Überdruck. 


Das vierte Profil (#025; b=0,05) ıst das zuerst berechnete, das hauptsächlich zur 
Aufklärune der Verzerrungsverhältnisse der konformen Abbildung dienen sollte: bei der Wahl 
von x ließ sich das Bild der Hinterkante nicht in einen Teilpunkt des Quadranten legen, so 
daß dieser Punkt der Kontur nur graphisch bestimmt werden konnte. Eine Abschätzung des 
Gesamtfehlers ist daher bei diesem Profil nicht möglich; da aber die Rechengenauiekeit nicht 
viel geringer wie bei den übrigen Profilen war, ist die Annahme wohl berechtigt, daß dieser 
Fehler von derselben Größenordnung ist. Bei diesem Profil wurde die Druckverteilung auch 
noch für eine Anzahl anderer Anstellwinkel berechnet. 


bıs auf einen Maßstabfaktor, direkt die 











Yo 18°) 
Pr 2040" u 1.1910, c 0524| u 1,2399, € 0,536 u 6" 1313 v9. ( 0.553 
Saueseite Druckseite Saugseite Druckseite Saugseite Druckseite 
r Y 7 Y JE Y 2 Y JE Y 1 Y 
V.O010 0B.OOSOL — 0.00OOL 0.0032 0,0013 0.0079 0.0002 0.0032 1 0,0015 0.0079 — 0.0003 0.0032 
0.0241 0.0616 +-0.0056 0.0351 0.0261 0.0609 +0.0042 0.0349 0.0271 VOBOH -++0.0035 0.0346 
0.1809 0.2334 0.0860 0.1485 0,1967 0.2283 V.O806 0.1478 I 0.2008 0.2257 0,078t 0.1457 
- >» - >» + >» 
0.7947 0,5905 01.3070 0.3010 0.8102 0.5683 0.255. 03014 0,8189 0.5525 0,2900 0,2943 
3.104 1.587 0.6133 0.4437 | 3.723 1.451 0.5944 0.4429 1.595 1,308 0.5825 0,4276 
S.451 >,585 1.3246 0.6770 | 8.057 2,210 1.2554 0,6723 | 7,227 1.791 1.2464 0.6345 
IS.411 3.925 3.258 1.092 16,170 2.949 3.126 1.065 12.685 1.991 2928 0.955 
34,909 5,029 9.986 1.884 127,002 3.099 9.077 1.747 118.020 1.668 7626 1.379 
58.016 1.867 19.842 | 2516 138.284 2.226 16.672 2,195 122.002 0.926 12.391 1.533 
. ” - » En} >» 
76.72 1.3 13.49 3.22 15.31 0,947 30,90 54 23,98 0.228 18.96 1.91 
106.8 0.50 94.33 3.26 54.24 0.786 | 50.60 2,43 6.20 0 —0,583 | 29,18 1.30 
134.5 2.45 134,9 2.45 60,73 2.14 60,73 2.14 27,11 1.19 21.11 1.15 

















Betrachtet man die bereehneten Koordinaten, so fällt auf, daß sich beim Übergang von 





einem zum anderen Profil nur die Koordinaten der Punkte in der Umgebung der Hinterkante 
stark ändern, während die Änderungen am Profilkopf nieht sehr bedeutend sind. Es scheint 
demnach, daß eine starke Veränderung des Hinterteiles des Profils die Strömung am Profil- 
kopf nieht wesentlich beeinflußt. Trifft diese Annahme zu, so könnte man auch dem dicken 
Profil, das in seiner jetzigen Gestalt natürlich unbrauchbar ist, durch geringe Veränderungen 
des hinteren Teiles der Saugseite eine strömungstechnisch verwendbare Form geben. Die 
Entscheidung über die Riehtigekeit unserer Vermutung muß dem Experiment überlassen werden. 


Den Betrag des Auftriebes erhalten wir aus (28); bei der Anwendung dieser Formel ist 
zu beachten, daß bei ihrer Ableitung die Konstante € in (5) gleich Eins gesetzt worden war. 


Setzt man den Wert für Z" aus (6) ein und nennt die aus den Tabellen zu entnehmende 
';, Der Doppelpfeil bedeutet bei der Saugseite das Ende der freien Stromlinie, bei der Druckseite die Lage des 
vorderen Staupunktes. Die Tabellen sind stark gekürzt. 


Im 


be 
be 


wi(* 


che 


ZU 





“2 


AN 


eh 


R 
12, Heft 5 . a j . ’ _ vor .2 s aan 
Ba mr 1939 Schmieden, Die Berechnung kavitationssicherer Tragflügelprofile 30%; 











'rofillänge I,, so erhält man für den Betrag des Auftriebs (y = Riehtung der Geschwindigkeit 
ven die #-Achse): 

u 16 b° | u Eu 
A=o:2nlv, sin (u+y): 1 2rolv; sin(u+y)-e. . (88), 


2\ J 0 


| "+ ) 


obei 1 die wirkliche Profiltiefe bedeutet. Der erste Faktor ıst gerade der Auftrieb einer 
henen Platte von der Tiefe !, die den Anstellwinkel «+ y besitzt'’), der zweite Faktor e 
t von der Größenordnung Eins und ist ebenso wie der Wert von « in den Tabellen 


nzegeben. 


16. Der Fall b, =-0; ein Spezialfall. Die bisher berechneten Profile leiden noch an zwei 
helständen: erstens wirkt sich die Singularität der Krümmung besonders bei diekeren 

’rofilen ungünstig aus, zweitens ist die freie Stromlinie so lang, daß sieh der Druckanstieg 
ıf einem verhältnismäßie kurzen Stück der Kontur vollziehen muß. Man muß also ver- 
chen, die Länge der freien Stromlinie zu verkleinern; es ist zu erwarten, daß Profile mit 
ndlieher Krümmung an den Ablösungsstellen diese Kigenschaft besitzen werden. 

/uerst untersuchen wir wieder den Grenzfall b,=0, der sieh ins Unendliche erstreckende 
Konturen liefert. Dabei beschränken wir uns auf den «den letzten Abschnitten entsprechenden 
2. a | 
"all DB, ==0. Drückt man die Werte von «a,, a,, a, nach (35a) dureh b,, „ aus und trägt das 

ELLE 
I») 


Kreebnis in (22) ein, so erhält man aus diesen Gleichungen: 


n - a=1-+b e,=0; = 6 


* 
ww 


ı /1 +’ (1 -+%°) | 
! 
I | | | 4° 2x 
| I+sıng, 
ın To ! In Pr “ * _a,sing, —a,sin3gQ,. 
vn 2 | sınYy, 


Der Zusammenhang zwischen 5b, und z ist ein ganz anderer wie der zwischen «a, und x 
ım vorigen Abschnitt; nahm dort a, mit wachsendem # monoton ab, so wird jetzt b), an 
beiden Grenzen des Intervalles für z positiv unendlich; dazwischen muß b, ein Minimum 
besitzen, das ungefähr bei 9, 40° liegt (vgl. Abb. 15). Der Wert dieses Minimums ent- 
scheidet über die Möglichkeit von technisch brauchbaren Profilen dieser Klasse; wird nämlich 

b,.;n größer als 0,5, so wird 3a, größer als «a,, das Profil bekommt also in der Umgebung 
des vorderen Staupunktes eine Beule. Nun wird zwar bin ungefähr 0,3, der Bereich der 
‚ulässigen Werte von x bzw. g, wird jedoch stark eingeschränkt. Alle Werte von g,, die 





b, >0,5 liefern, scheiden von vornherein aus; der übrigbleibende Bereich unterhalb der 
vestrichelten Geraden in Abb. 15 wird noch weiter verkleinert, weil schon sehr kleine 
dr 
20° 
% A 
7,0 2 
u, ——— ’ 
0° 75° 30° 45° 60° 75° =%, 2 15° 30° 4% o 50° 7% r. 
Ranoi e Rämozn 
Abb. 15. Abb. 16. 


1?) Vel. Fuchs-Hopf. Aerodynamik. 8. 72. Gl. (92). 
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Werte von 2 möglichst zu vermeiden sind und die Geschwindigkeit 4, innerhalb der b« 
(ta 


kannten Sehranken bleiben muß (Abb. 16). Wenn wir noch das sogleich abzuleitende Resulta! 
vorweenehmen. daß ’n mit wachsendem l hei festem y.4 zunimmt, so bleibt für Yo nur de 
Bereich von > 30" bis = 50° übrig. Für die numerische Rechnung kommen innerhalb diese: 
Intervalles nur eine Anzahl diskreter g,-Werte in Frage, die einem brauchbaren Teilpunk 
des Quadranten entsprechen. 


Für von Null verschiedenes b bzw. « wurden die Gleichungen (22) in den beide: 
Fällen 4, 93S°15’ und g, 49" gelöst. Das Ergebnis zeigen die beiden folgenden Tabelle: 
und «lie Abbildungen 17 und 18. Besonders aus diesen Abbildungen sieht man, daß di: 
Kurven für den Kantenwinkel wie für die Gesehwindiekeit auf der freien Stromlinie nahezu 
oeradlinie verlaufen; nur für sehr kleine «. weichen sie von einer Geraden ab. 





Go 
v4 j 

15 P 

y;45° 
38915 =. 79°716' 

ia % 30°, y,-36°% 
14 20° 

12 70° 

10 | | j ER j | gr 

g0 n- y0 g° g° we 4 0 2° 14 9 5° g® 10° 

Aırso 71 ; MAraozuy) 

\bb. 17. Abb. 18. 


Für die zu y, 98°15’ gehörige Profilserie wurden die Koordinaten der Kontur, die 
Geschwindiekeitsverteilunge und der Auftrieb numerisch bestimmt. Die Unterteilung des 
Integrationsbereiches war im wesentlichen dieselbe wie bei den im vorigen Abschnitt be- 
reehneten Profilen: nur in der Umgebung der Hinterkante wurde sie nochmals verfeinert. 
Das Ergebnis der Integrrationen zeigen die Tabellen und die Abb. 19 bis 22. Besonders auf- 
fallend ist die geringe Empfindlichkeit des Profilkopfes gegenüber großen Änderungen der 
Umgebung «der Hinterkante Die Länge der freien Stromlinie beträgt maximal 20°/, der 
Gesamtlänge, so daß diese Profile im wesentlichen allen Forderungen genügen, die man an 
ein kavıtationssicheres Profil stellen kann. 








y 38° 15’ 45° 

f 20 40° j* 6° 90 240 4° 6" U 

a 724’ 11" 40’ 180 38° 30048’ 5° 20’ 8’ 10’ 12° 53° 16° 11' 

N 

a" 1" 1° 40’ 3° 5° 40’ v1’ 00 13’ 0016’ 0 30° 
ZN 0,3624 0.1006 (0.4667 0.874 0,4167 0,4360 0,4684 0,4935 
+, 0,5081 0,5210 0.5400 0,5682 0.4204 0,4317 0,4485 0.4595 
I (a 1.2741 1.2786 1.2540 1.2898 1.1982 1.2080 1,2227 1,2326 

—e 1.386 1.145 1.526 1.630 1.5838 1.590 1.670 1.714 
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=" 15 
u 2040’ 1.386 v0, c=0,516 u 4° 9% 1,445 vv, c=0,531 
Saurseite Druckseite Saugseite Druckseite 
rl Y Mn 4 LE U z Y 
0.0017 V.OOTS 0,OOHU 0,0124 0.0017 VOOTS V,OOHO 0,0124 
DRLIEIE DOGS 0.0679 0,0679 0.0031 VOGHIN 0.0676 0.0687 
“ >» + ® 

(LUSSA 0.2234 0.3608 0.2381 0.0012 0.2224 0,3704 (1,2440 

0,3625 0,4928 0.7124 0,3007 0,3701 0.4871 0,7108 0,4011 

1.1561 0.0796 1.5581 0.6916 1.1704 0.0541 1.5439 0,7077 

2.6691 1.5871 2,4982 0.9621 2.6708 1.5094 1.3158 1.4087 

1.1756 2.0105 1.4402 1.4029 1.1297 1.8685 8.691 2,140 

- > 4 » 

6.572 2.451 6.227 1.736 6.385 2.193 13,219 2,700 
14.79 3.114 3,176 2.197 13,54 2,445 21.429 3.431 
32,00 3.224 24.29 3,847 26.43 1.860 36.75 4,293 
52.81 2.54 16.06 36 10.63 0,599 61.95 4.81 

113,5 1.04 110,1 6,98 63.35 2,24 19.8 1.68 
176.1 1.39 176,1 1.935 79.3 4,68 
Yo 838° 15 
u 6° 0, 1.5260, € 0,5 u 9° u 1.693907, « 0,574 
Saurseite Druckseite Saugseite Druckseite 
l 4 r 4 r 4 l Y 

0017 VOOTD 0.0059 0.0124 0,0016 VUOOTS V.OOGO 0.0125 

VUO0O3D VOHGIS V,UHSB 0.0699 0,0049 0.0619 0,0695 0.0716 
P) » - > 

0.0954 0,2210 0,3707 0.2524 0.1018 0,2192 0.3698 0,2640 
0.3816 0.4785 0.7055 0,4150 0.3901 0.4654 0,6917 0.4320 

1.1887 0.0136 1.5097 0.7251 1.2084 0,8478 1.1340 0,7361 

>. 6460 1.3541 1,0468 1.3860 2.5607 1.1853 3.2882 1.2986 
3.9975 1.6458 1.699 1.986 3.6923 1.3066 6.018 1.695 

- > - > 

2,962 1.803 16.105 2.182 2.157 1.277 10.55 2.017 

7.803 I.SUN 31.84 135 6.3753 1.139 15.36 1.079 
11.404 1.595 36.41 3.01 s.412 0.756 16.39 1.01 
18.00 HU 11.678 0.132 
32.28 2.040 15.52 1.524 
36,1 3,01 16.39 - 1.90 

Zusammenfassung. Es wird ein Verfahren zur Berechnung kavitationssicherer Profile 


entwiekelt und auf die explizite Bestimmung einer Reihe von praktisch brauchbaren Profilen 


aneewandt. 


Es sollen hier noch zwei mögliche Erweiterungen der Theorie erwähnt werden, 


die eine weitere Verbesserung der Profileigenschaften ergeben könnten. Einmal kann man 


in (15a) für die a, komplexe Werte zulassen, so daß die Geschwindigkeit auf der „freien 
Stromlinie* nieht mehr konstant bleibt. und dann die Konstanten so zu bestimmen versuchen, 


daß man einen günstigeren Druckanstieg auf der Saugseite erhält. 


Es wurden eine Reihe 


von Profilen der Klasse a, +0 bereehnet: jedoch verziehten wir auf eine Wiedergabe dieser 
Profile, da die Verbesserungen gegenüber den entsprechenden Profilen mit reellem a, nicht 
Weiter kann man versuchen, durch entsprechende Formgebung der Druck- 
seite, die bei den hier untersuchten Profilen zum Auftrieb recht wenig beiträgt, eine Ver- 
Über diesbezügliche Untersuchungen, die noch nicht 
abzeschlossen sind, soll in einer späteren Arbeit berichtet werden. 


wesentlich sind. 


besserung der Profile zu erreichen, 


1S0 
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KLEINE MITTEILUNGEN 


Kurze Bemerkungen zur äußeren Balli- 
stik. (Graphische Bestimmung des Abgangsfehler- 
vinkels bei beliebigen Erhöhungswinkeln eines Ge- 
chützes. Vorsiehtsmaßregeln beim Gebrauch der 
abellen von Fasella für kleine Werte von 


a 
2 


1. Schießt man aus einem Geschütz unter einem 
rerebenen Erhöhungswinkel A (Abb. D), und legt 
man einzelne Punkte der so erhaltenen Flugbahn 


= 





Abb. 1. 


in Ihrer 
Winkel & 


l,aare im Raum (Abszisse .r. Ordinate y,. 
der Visierlinie zu dem Punkte mit dem 
arete J 

Mn 
kann man die Aufsatzwinkel 4 e als Funktion 
der Abszisse .r in einem kartesischen Koordinaten- 
system darstellen (Abb. 2 Punkt 1. 2, 3. 4, ...) und 
die so erhaltenen Punkte mit Hilfe eines Kurven 


Horizont & photogrammetrisch fest, so 


(h-e)” 


Äurvenlage 








Abb. 2. 
holzes dureh eine Kurve verbinden. Der Winkel 
h—e stellt. genauer gesagt, die Projektion des 


Winkels zwischen Visierlinie und Seelenachse 
auf die Vertikalebene durch die Seelenachse in 
derjenigen Lage, in welcher sie sich vor dem Schuß 
befand. dar. Der Winkel « (Abb. 1) zwischen der 
’rojektion der Visierlinie auf die Vertikalebene 
dureh die Anfangstangente 7 der Flugbahn und 


dieser Anfangstangente ist bekanntlich um den 
unbekannten Abgangsfehlerwinkel ö größer. Zur 


Entfernung 20 gehört aber ein Winkel « =V. 


somit muß die Kurve Il in Abb. 2 die Ordinaten- 
achse £=0in der unbekannten Ordinate h — & Ö 
schneiden. 

Nun entnimmt man aber einer Bemerkung in 


Cranz,Ballistik 1. 5. Aufl.. 8.531. daß die Kurve Il 


für r=0 einen 
Neieuneswinkel 


vom Luftwiderstand unabhängiren 
eeren die \bszissenachse besitzt. 


Jessen tanzens ‚_ „ Ist. (Wenn h—e in Grad- 
A) Fi 
mab,.z in meter aufgetragen ist. ist der Neizungs 
winkel 
7 IS) v 
h £ | —e .c) 
2 dp I 


are |te 


ns 


Z meteı 


Aus dieser Bedingung ergibt sich die Möglich 
keit einer sehr genauen Bestimmung des Abezanges 
fehlerwinkels ö, wenn die Flugbahnpunkte mit der 
üblichen Genauigkeit aufgenommen und vom Wind- 
eintluß befreit sind'). 


Man lege nämlich (Abb. 2) ein biegsames Kurven 
holz durch die Punkte 1.2, 3.4 usw. (Kurvenlage ]), 
schiebe dann von unten her ein Dreieck, dessen 
Re y 180 i 
Seite PQ unter eeeen die 


. 


2 

neigt ist, parallel mit sich gegen das Kurvenholz 
hin und drücke das Kurvenholz elastisch nach 
oben, bis der Berührungspunkt zwischen Kurven- 
holz und Dreieck in die Ordinatenachse fällt (Kurven- 
lage II). Die Lage des Berührungspunktes auf 
der Ordinatenachse ergibt den Abgangsfehler 
winkel 6. 


r-Achse ge 


2, Die Beziehung wird eewöhnlich aus der 


= Un" 

Gleichung der Flugbahnparabel mit gleicher An- 
fangszeschwindigekeit =, im luftleeren Raum .her- 
eeleitet, mit welcher die ballistische Kurve in ihren 
Anfangspunkt drei unendlich nahe Punkte gemein- 
sam hat. -Wir wollen hier Beziehung aus 
einem Grunde, welcher im nachfolgenden klar 
werden wird, aus den Siaceischen Flurbahnfor- 
meln herleiten. Diese Formeln lauten bekanntlich 
(vel. z. B. Cranz 1,5. Aufl.. S. 187 


(diese 


cz [Alk)— A(v,) / | 
/ r to ( x P BER 
' 7 2 008°? 4 |D(u) D(v,) 
und 
2 
, Dim) D(n,) 2 
Pr 
Dabei ist 
udu du 
Diu - 5 J(w) 2q | 
Sstrım Juflu) 
In) u 
I) I. 
. /ım) 


in) ist der Verzörerung dureh den Luftwiderstand 
proportional und 4-=9.8047 die den Tafeln zugrunde 
eelerte mittlere Schwerebeschleunizunge. 


Aus (1) und (2) folgt zunächst 
y=rtigyg | 


J(v, Din 
Di Va)” 


x? Inu 4(v, 


2c08s°gq Din 


und mit 


Y 
- to: t2.(g a) 


I) Entstammen die Punkte 1, ?, 3, . verschiedenen 
Scehiissen, so müssen sie vor ihrer Auftragung auf gleiche 
"og und auf gleiches Luftgewicht reduziert werden. 
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wird hieraus 
to 4 [7 te q | 
I Inu Ai Jir,'Diu Div, | . 8). 
c0s°’gq Din I) Vo) ° 
Für kleine «a ist 
[Fi 
to 7 (! tt7gq - 
COS q 
und aus 
lu A(v J(%,) Di) Di(v,) 
Din D(v,)\° 
wird in der Grenze 
lım u A io J/(v,) Dia) Div,)| 
Hu "o Dim) D(v,) - 
Somit folgt aus (3 
r ‘({ 
lım 
\ () \ rl 
Hm 4(v,) J(v) Di Div, 
2 u=i Din) Div)? 
Da der rechtsstehende Grenzwert für u DM 
) - i \ 
zunächst die Form annimmt. differenzieren wir 
Zähler und Nenner nach «#. Mit den oben ange- 
führten Werten der Interrale A. J und D folgt 
Jin) Hu 
uU J(v,) 
a iu) fu) 
lim lım 
. ( lır & ’ 7 
) „H 0 2 D u D(rv,) 
(2) 
1 » JS J(%,) 
lim 
® u=ra 21DiW) Div,) 
\lso noeh immer 
€ H Jita) v0) 
lim Rind (4). 
u— ro 2|Diw Din) 0 
Abermalize Differentiation nach x ergibt 
1 
2g 
: r ON „H r ( 
lim lim lim J 
x 0 iz 2 u "o 2 H Zu vo U“ 
tim 
Somit !st 
da ] (] 
d.r 2 © 
3. Wir fanden 
Km AI Aiv,) J(v, Diu Den,) [7 
u "on Diu D 2) 2 vn 
I) A [BAR /(v, Diu Dir) 
Diu Dir, 


ist aber die Fasellasche Funktion },(”,. = ). 
2 


Die Fasellasche Funktion 


‚( & Alu) — A(v,) — J(v) [Did — D(w,)) 
[2 “ = f 
ea D(u) D(v,)| 
2\ a 
f; (v, e) ep’ 
hat dagegen für _, 0 den Wert 0, 
Vergleicht man das Anwachsen von f und 
von /, für kleine ,, so findet man, daß f, (v0, ‚| 
2. 


nahezu linear wächst. die zweiten Differenzen wer- 

den so klein, daß sie die letzte Stelle des linear 

interpolierten Funktionswertes bei genauerer Inter- 

polation unter Zuhilfenahme der höheren Differenzen 
. * * } L 

nicht mehr beeinflussen. Bei Av. ‚) begeht man 

B 


aber bei linearer Interpolation für ,< 100 unzu- 
er 


lässig große Fehler, weil die zweiten Differenzen 
beträchtliche Werte annehmen. Man erhält also 


Am Dr 
Ivo, = für ee 100 viel genauer, wenn man 


I\ . : 7_: 
-) in bisher gewohnter Weise berechnet. 


f; | "o. 


( 


und dieses mit _, multipliziert, als wenn man direkt 
r 


in die Tafel für / | er eingeht. Und gerade beim 
r 


Rechnen in Teilbögen kommt man öfters in die 


l,age. Werte für , <100 aufschlagen zu müssen. 
er 


Um über die Größe des Fehlers, welchen es hier 
zu vermeiden gilt, Aufschluß zu bekommen, be- 
zeichnen wir | mit 4:7, (9... 0) mit az}, (vo, 100) 

ce T 100 


mit b. Die zweiten Differenzen sind bei /,. wie 
man sich leicht überzeuren kann, so klein, dab 
man, ohne einen Fehler in der letzten Stelle zu be- 
zehen, linear interpolieren darf, also setzen darf: 


Kn 


', (vor) a- 


Ib da). 


Da nun 


ist. wird somit praktisch genau: 


; TE a er ’ 
\,—;)=1W04la+4(b a), 
c 


In der bisherigen Weise arbeitend würde man 
ddlareren direkt in die Tafel für v, ‚eingehen 
ei 


und dort finden /(2,.0) = 0: F(r,. 100) 100 5b und 


z 1 i 
setzen: Ivo. — 1004-b. Der Fehler. den man 
vr 


bereht, ist also 10 2b —- [100 Aka -—+-A(b - a))| 

10 (b -a)A(1 4). Der relative Fehler wird 
100 (db a)\7-(1 7) b-a)(i 7) 
1WAla+ib -a)) a+ilb-a 

u * .. ”. * a ”. T 

Er wird am größten für kleine . also für 


x () und wächst, wie man aus Zahlenwerten für «a 
und b entnehmen kann, auf über 3°/.. 


Ka 279 





ew. 
ch. 


nd 


e 
eT- 
ar 
ET- 
‚en 


an 


U- 


en 
Iso 


an 


im 
lie 


n. 


er 
IP - 


0) 
ie 
db 


e- 


Kan 


> 


... 
. .€® 
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d_ ; . f, ist aber das Siaceische T(«) — T(w). Dem- 

Bei 5 - 100 bleibt der Fehler schon unter 1a ury nach ist 


ann also vernachläßigt werden. 

Beispielsweise findet man in Fa sella bei linearer 
ıterpolation F (600, 35) 0,000 975. Unter Zu- 
‘Ifenahme der höheren Differenzen würde man 
500. 35) == 0.000 961 finden. 

Dareren ist /, (600, 35) = 0,0000 2746 bei linearer 
ıterpolation und f, +35 = 0,000 961. (Die höheren 
ifferenzen bei /, ändern das Resultat 0,000961 
icht mehr ab.) Den gleichen Erfolg erzielt man 
‚uch. wenn man nieht mit Formel (1), sondern mit 
'ormel (1a) rechnet, die ja auch 


Dr Mn 
I, | Vor u 
( 


Y ig > 608? 4 


-sehrieben werden kann. 
{, Wir haben oben Gl. (4) 


vn Ju) J(v,) J 
Hl "o > ID ( u) D (Vo) 0° 
hestimmt. Demnach ist 
lim J (u) J (vo) ki Jin) J(v,) 2g 
U=Vo Dim) Di ®g) Uetg LE do“ 
:=08 ce’ 


Nun ist /(a) — J (ev) gleich der Fasellafunktion 


j 


Mn Pr vo . F : 
lvo oe) und auch /, wächst für kleine r keines- 
| z 
I | "gs e 


” 


Be, ’ 
wees linear mit „. Bildet man dagegen 
r 


( 
so findet man. daß die zweiten Differenzen viel 
j E fl, (9,0) 2gy 
kleiner bleiben. während’ z „ genau defi- 
’g“ 
ce’ 


niert Ist. 
Man wird deshalb auch bei der Fasellaschen 
Tafel VI genauer arbeiten. wenn man über die 
: 2004. 
"unktionswerte 7, (2,.0)=0 die Werte B a 0m 
0 


für allemal anschreibt, zwischen diesen Werten und 


denjenigen für 100 interpoliert und das Inter- 


= j 
‚lationsresultat mit / multipliziert. 
r ce’ 7100 I 
7. B. ist in bisheriger Weise linear interpoliert 
/, (600. 35) — 0,00197.. 


Hat man dagegen über 7, (600.0)=0 den Wert 
200.9,8047 


600? 
‚wischen diesem und 7, (600,100) = 0.00564. so er- 


=, 00545 angeschrieben und interpoliert 


2 " r u . 
hält man für „==35 als Interpolationsresultat 
e 
(100552 und dieses mit 0,35 multipliziert ergibt 
00193. Der relative Fehler kann bei /, bei der 
bisherigen Methode bis zu 3,85 °/, anwachsen. 

>. Auch die Fasellasche Gleichung für die 
"luezeit 


ec’ 


t (ro. ;) 


COS q ( 


‚ besser 


schreibt man für kleine 


I, |» . Zu. 
Mi | 2) ec’ ] (v1) fj (24,) 
r DW) - Du, 
e’ 
2.2 "du 
mit T (a) ; 
ru) 
Auch hier wird 
Tu) — T (u) 0 
lim 
„H "oo Din) Diu,) 0 
oder differenziert 
1 
? (u) | | 1 
lim lim - 
[2 "oo u „H "og u u 
(u) 
somit 
l.r 
It» 
!ı COS Y 
x 
ACH} 51 r 
Da auch linearer zunimmt als Al. , 
z \ Ü 
M 
( 
" : x 100 
tut man gut, über 7, (7,0) 0 die Werte von 
Am 
anzuschreiben, und für , I00 zwischen diesen 


- 
Werten und 7, (®,. 100) zu interpolieren. Das Inter- 


polationsresultat mit , multipliziert ergibt 
ce 7 100 


=: - . 
den genaueren Wert für /, | "os |)» Bei dem Ver- 
\ ®” / 


100 


eleich der Werte für mit /,(7%. 100) wird man 


0 
oft finden, daß letztere Werte bei Fasella um 
eine Stelle zu wenig genau sind. Es empfiehlt 


: 1 

sich deshalb. gelegentlich noch = ; (v9. 200) und 
| 

5) 


«) 


/, (2,. 900) zu bilden. Ein Blick auf den Verlauf 


2 . 
‚=Vd bis „300 zeigt dann 
= . 


wie 7, (9. 100) zu berichtigen ist. 


dieser Werte von 


Im allgemeinen empfiehlt es sich, statt der stark 
abzerundeten /,-Werte auf die eizentliche Siacei- 
tabelle zurückzugreifen. bei welcher allerdings 
zweckmäßig auch die «-Werte um eine Stelle er- 
eänzt werden. 


6. Mit Rücksicht auf die Betrachtungen unter 
Ziffer 3 erscheint es empfehlenswert. in dem Aufsatz 
des Verfassers (Zeitschr. f. anzew. Math. u. Mech.. 
August 1931. Heft 4. S. 254) die ohnehin mit einem 
Druckfehler behaftete Formel (5). die richtig lauten 
mub 

CyE£ So Yu = lo. J > 


y=-zrigg = . 
2cos’p Ss? 98047 \ Cy 


zu ersetzen dureh die ohne weiteres aus ihr folgende 


2° rd 2 
- L f, w 0» ‚| ? 


y—=ıig9 eier — 
2cos?’p s? 9,8047 Cy 


also mit /, statt mit / zu arbeiten und auch das 
Rechenschema entsprechend abzuändern. 


0. v. Eberhard in Essen-Bredeney. 283 
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BUCHBESPRECHUNGEN 


(Die hier angezeigten Bücher sind dureh die VDI-Buchhandlung Berlin NW 7, Ingenieurhaus, zu beziehen.) 


Dr.-Ing, GUSTAV FLÜGEL, o. Prof. für Dampf- 


turbinen, Propeller und Strömungsphysik an der 


Teehn. Hochschule Danzig. Die Dampfturbi 
nen, ihre Bereehnung und Konstruktion. 
Mit einem Anhang über die Gasturbinen. 
\All 324 S, mit 297 Fire. im Text. Leipzig 1931. 
Joh. Ambros. Barth. Preis geb. 37M. 


Das äußerst interessante Buch legt «las Haupt 
vewieht auf die Beherrschung «ler 
Dampfturbinen. wie sie dem Berechnungsinge 
nieur umd auch dem auf dem Prüfstand oder bei 
UÜbergeaben tätieen Ingenieur not tut. Ohne die 
Theorie weiter zu treiben. als es für die Technik 
erforderlich ist. wurde sie dennoch in überlerener 
Weise überall zugrunde gelegt. wobei alle be 
kannten feineren Betrachtungen der Literatur be 


reechnerische 


rücksiechtiet und auch viel Neues gebracht wird. 
Insbesondere als neu hervorzuheben ist «lie für 
den Ingenieur so besoniers wiehtire,. ausführliche 
Behandlunze der Strömunges- und sonstiren Ver- 


luste,. weil sie dem Konstrukteur notwendige Hin- 
weise eibt, ferner die Theorie des Verhaltens der 


Dampfturbine bei geänderten  Betriebsverhält 
nissen. «die dem  Prüffeld und dem Betriebs 
ingenieur wiehtie ist, während «lie Näherungs 


rechnung bei Änderung der Bauart wieder «den Ge 
stalter und den projektierenden Ingenieur betrifft. 
\uch mehr allgemeine Aufgaben, wie die Dampf- 
expansion in Düsen, die Radreibung. der Gefälls- 
rückgewinn u, a, werden in neuartiger Weise ge 


klärt, ebenso erhalten bauliche Frazen. wie «die 
Berechnung flacher und zewölbter Böden. die Be 
reehnunz von Schweibensehwineungen. neue Be 


euehtung., 

Sehr eut und teilweise 
die zur Erhöhun®e «der Wärmewirtschaftlichkeit 
dienenden KEinriehtungeen behandelt. erst «ann 
tolet die Neubereehnung und der Bau «der Dampf 


auch neuartie wenden 


turbinen. ferner «die wiehtiesten Festirekeits- und 
Schwineungsreehnungen und emdlich ein kurzer 


\nhane über den Stand des Gasturbinenbaues. 

\lles in allem ein Werk. das (lem höher stehen 
den Ingenieur wirkliche Frewde und Nutzen be 
reitet und für «den Studierenden in seiner Knapp 
heit sicher hochwillkommen ist, 
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FRANZ OLLENDORFF, Berlin. Potential 
[elder der Elektroteehnik. VII 
mit 244 Abb. im Text. Berlin 1932. Julius Springer. 
Preis 32 M. 

In noeh stärkerem Maße. als 
heren. auf jeweils engere Gebiete beschränkten 
Büchern Verfassers «(ler Fall war. zeigt das 
vorliegende dessen grobe Fähiekeit. die verschie 
lenartiesten Fragestellungen «der Elektrotechnik 
für eine mathematische, speziell potentialtheore 
Formulierung zueänelich zu machen. Die 
dabei notwendigen Vereinfachungen der physika 
lisehen Verhältnisse sind zwar sicher nieht immer 


395. 


dies in «den frü 


(les 


tische 


leieht einleuehtend und werden manehmal (lie 
eleentliche Schwieriekeit bei «der Lektüre «les 
Buches biltlen, während «ie mathematischen An 


das hin 
(ruppe 


forderun®en selbst «(durehaus nieht über 
auseehen. was der «dafür interessierten 
von EKlektroteehnikern zugemutet werden kann. 
Jedtenfalls aber zeiet sieh überall in Buche 
eine eanz individuelle Einstellung des Verfassers 
‚u den behandelten Aufgaben. ein Streben, in der 
mathematischen Formulierung mehr zu sehen. als 


lem 


ein bloßes Hilfsmittel zur auantitativen Beherr 
sechune, was eerade im Hinbliek auf den techni- 


sehen Ausgangspunkt zu begrüßen ist. In diesem 


durchaus individuellen Charakter des Buch: 
scheint mir sein Hauptwert zu liegen, weniger | 
der Anordnung der vielfältigen Aufgaben. di: 
nieht nach der physikalischen Zusammengehörig 


keit. somdern mehr nach dem mathematische: 
Schwierigkeitsgrade erfolet. Beeinnend mit de 
Konstruktion von Potentialfeldern aus angenon 


menen Ladungsverteilungen und dann zu Ranıl 
wertaufgaben verschiedener Art übergehend: dies: 
letzteren werden unterteilt nach der Anwendung 
komplexer Funktionen. sowie nach der Koor 
dinatenwahl bei den räumlichen Problemen. Siehe:ı 
lieet hier weder ein schulmäßiges Lehrbuch dei 
P’otentialtheorie. noch ein elektrotecehnisches Lehr 


buch vor. Aber ein Buch. das den eleiehartie ein 
restellten Lesern sehr vieles und auch Instruk 
tives zu bieten vermag. und das auch Mathe 


(ebiet interessiert sinil. 
F. Noether. 2% 


die für «lieses 


eeben kann. 


matikern. 
Änrerung 


Prof. E. G. COKER. M. \., D, Se. F 
Prof. L. N. G. FILON, M. A. D. Se. F. R. S. 
\ Treatise on Photo-Elastieity. Cam 
bridge University Press. 1931. XVII + 720 S. mit 
vielen Abb. und mit 14 farbigen Tafeln. Preis 
30 sh netto. 


‚R.S.. und 


Insti 
hier 
Gegen 


Die Autoren, die an zwei verschiedenen 
tuten der Universität London wirken. haben 
die erste umfassende Darstellung des 
standes wereben. Das Buch gliedert sich in 
“ Teile. Der erste Teil ist eine kurze systema 
tische Darstellung «der physikalischen Optik vom 
Standpunkt der Maxwellschen elektromagnetischen 
l,iehttheorie. Die zweite bringt die mathematische 
lastizitätstheorie mit besonderer Berücksichti 


eunge «der ebenen Spannunges- und Formände 
runeszustände. Hier werden u, a. auch die 
nötigen geometrischen Vorarbeiten für die Be 
spreehung der photoelastischen Spannungsunter 
suchung szeleistet. Die eigentliche Darstellung 
der Methoden der Photoelastizität geschieht im 
Ill. Teil, und zwar im wesentlichen an Hand «der 
eeschichtliehen Entwieklung: hier sind auch (die 
interessanten Arbeiten von Maxwell u. a. übeı 


entsprechende Erscheinungen bei Flüssigkeiten 
eineehend berücksichtiet. Die weiteren 5 Teile 
enthalten die Anwendungen. Hier ist die Gliede- 
rung der umfangreichen Materie nur nach prak 
tischen Gesichtspunkten vorgenommen. Ein Teil 
ist der Balkenbiegung. einer von Kreisen und ge 
raden begrenzten Scheiben. ein weiterer Teil Kerh 


problemen sewidmet: «dann folet ein wichtige: 
Teil über Anwendung der Photoelastizität zu 


Klärune von Problemen der Materialprüfung (Aus 
wertung der Festigkeitsproben) und das Schluß 
kapitel über die Anwendungen auf zusammenge 
setzte maschinenteehnische Probleme und bau 
statische Aufgaben. 

Die Darstellung ist ausgezeichnet und besonders 


dlank «der zahlreichen zuten Abbildungen Jeich! 
lesbar. Bei der Auswahl des Stoffes fanden ın 


erster Linie die grundlegenden. klassischen Unter 


suchungen sowie «die eirenen Forscehungsergeh 
nisse der Verfasser und ihrer Schüler Berück 
siehtieung. Viele wiehtige Arbeiten blieben un 


herücksiehtiet. so z. B. die Arbeit von Mesmer, di: 
Untersuchungen von L. Föppl sowie die kinemato 
vraphische Methode von Tuz!. Seinem ganzen Auf 
bau und Charakter nach ist das Werk in ersteı 
Linie ein Lehrbuch der Photoelastizität und 
wird als solehes dem angehenden Fachmann au! 
diesem Gebiete auseezeichnete Dienste leisten 
Darüber hinaus dürfte es aber auch vielen Inge 
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ijeuren. Physikern und auch Mathematikern. die 
ich mit geometrischen Feldproblemen befassen, 
‚ertvolle Anregungen geben. Kine interessante 
‘euerung ist in der Bezeichnung der Formeln 
nd Abbildungen durehgeführt: aus der Nummer 
iner Formel oder Abbildung ist sofort zu er 


ehen. in welchem Kapitel und welchem Unter 
hsehnitte sie zu finden ist. 

Der verhältnismäßie nicht hohe Preis dürfte 
em Buche den Wege zum verdienten Erfolg 


‚esentlieh erleichtern. PaulNem6önyi 29 


Dr. phil.. Dr. ing. h. e. EMITL MÜLLER ;. o. ö. 
Professor an der Technischen Hochschule in Wien. 
Vorlesungen über darstellende Geo- 
netrie. Aus dem Nachlaß herausgege. von Dr. 
Iosef Leopold KRAMES, a. ö. Professor an der 
Deutschen Technischen Hochschule in Brünn. 
II. Bd.: Die Zyklographie. IX + 476 S. m. 208 Text 
ie, Preis 38 M. II. Bd.: Konstruktive Behand 
une der Regelflächen. VIII + 303 S. m. 155 Text 
fie, u. einer Tafel. Preis 5250 M. Verlag Franz 
Deuticke, Leipzig-Wien. 


Schon vor längerer Zeit sind aus «lem Nachlaß 


mil Müllers, von J. L. Krames herausgegeben. 
‚wei weitere Bände der Vorlesungen über «ar 


stellende Geometrie erschienen. >Sie sind zu be 
»rüßen als die fast einzie dastehenden Versuche 
einer  Weeiterbildune  darstellend - eeometrischer 
Methoden und Ergebnisse in unserer Zeit. Der 
erste «der beiden Bände behandelt die Zvklo- 
rraphie. das ist die Abbildung der Raumpunkte 
auf die Kreise der Ebene. in der erstmals von 
W. Fiedler angegebenen. später von Müller weit 
vervollkommneten Art. Der Hauptteil des Buches 
rjlt flächentheoretischen Untersuchungen. «lie aus 
len Beziehungen erwachsen. die zwischen krum- 
men Flächen und auf ihnen verlaufenden Kurven 
einerseits und den entsprechenden Kreismannig 
[altiekeiten in (der Ebene anderseits bestehen. 
ZAahlreiehe Ergebnisse, «lie gewiß. vom Standpunkt 
ler reinen Geometrie gesehen. nieht ohne Wert 
sind. werden auf diese Weise zewonnen. 

Der zweite Band (der dritte der ganzen Serie 
behandelt die konstruktiven Probleme der Regel 
tlächen. und zwar sowohl alleemeiner wie soleher 
vom 3. und 4. Grade. Hier sind viele Einzel 
heiten meu und von Interesse. Sehr erwünscht 
wäre in diesem Zusammenhange eine Untersuchung 
larüber. wie weit prinzipiell die Methoden «der 
(larstellenden. d. h. der konstruktiven Geometrie 
reichen: welehe Konstruktionsaufgaben sich mit 
len gewohnten Hilfsmitteln (das sind nieht mur 
/irkel und Lineal. da auch z. B. Schnitte einer 
(eraden mit beliebigen Kurven benutzt werden 
ausführen lassen. Solange derartige Fragen nicht 
beantwortet. ja nicht einmal gestellt werden. be 
hält die darstellende Geometrie immer noch «len 
Charakter einer Sammlung von Einzelaufgaben. 
die nieht «dureh einen genügend allgemeinen Ge 
lanken zusammengeehalten werden. 

Jeder an der Geometrie Interessierte sollte in 
liese beiden Bände. die übrigens vom Verleger 
ıufs beste ausgestattet sind. Einblick nehmen. 


Mises. 298 


Prof. Dr. phil.. Dr, jur. ALFRED MANES, Hono 
rarprof, a. d. Universität und an der Handels 
hochsehule Berlin, Versicherunieswesen. 
System der Versicherungswirtschaft. 5. völlige 
veränderte und erweiterte Auflaee in 3 Bänden. 
nit Einschluß der Sozialversicherunge. 2, Bi.: 
(»üterversicherung. B. G. Teubner. Leipziz unıl 
Berlin 1931. X +316 S. Preis geb. 20 M. 3. Bd.: 
"’ersonenversicherune. B. G. Teubner. Leipzig 
ind Berlin 1932. IX + 336 S. Preis eeb. 22 M. 


Wenn auch im Gegensatz zu dem in dieser Zeit 
sehrift (Bd. XI (1931) S. 401) besprochenen 1. Banile 
im 2. und 3. Band von mathematischen Dingen gar 
nicht die Rede ist. so erscheint eine Besprechung 
in dieser Zeitschrift im Interesse der mit Versiche 


rungesmathematik beschäftigten Leser doch sehr 
anzebracht. Ist doch für das Stwdium der Ver- 


sicherunesmathematik in erster Linie die >Sach- 
kenntnis auf dem Gebiet der Versieherungswirt 
schaft notwendie. Es gilt auch hier das von 
Tsehuprow geprägte Wort, daß der Mathe 
matiker auf dem Gebiet der Statistik die Sach- 
kenntnis haben müsse, um nieht nur mathema 
tischer Kunstsehlittschuhläufer zu sein. An diese 
Mahnung soll um so mehr immer wieder erinnert 
werden, als vor einieen Jahren eine italienische 
Versiceherunesmathematikerin in einem in der Zeit 
schrift für Versicherungswissenschaft erschienenen 
Aufsatz behauptet hat. «der Studierende könne 
nach der von ihr dargelerten Methode alle ver 
siceherungesteehnischen Aufgaben lösen, ohne von 
der Versicherung selbst etwas zu verstehen. Die 
(üterversicherune teilt Manes ein in Sachver 
sicherung und Vermögensversicherung, Jene zer 
fällt in die Hauptteile Transportversieherung. 
Feuerversicherune. Hagelversicherung, Viehver 
sicherung. Deliktversiecherung,. Wasserversicherung. 
sonstiee  Versieherungszweiee, unter denen vor 
allen «die von dem Mathematiker Höckner ins 
l‚eben gerufene allerdings jetzt verunglückte Sach 
lebensversicherung erwähnt sei. In der Feuerver- 
sieherunge gewinnen. wie erzänzend doch bemerkt 
werden soll. mathematische Überlerungen immer 
mehr Bedeutung. Bei der Transportversiecherung 
hätte das von Burrau in seinen Grundlagen der 
Versicherunesstatistik gerebene. sehr interessante 
Beispiel einer Prämienbereehnung. wie sie mit 
eutem Erfole in Norwegen während des unbe 
sehränkten U-Bootkrieres für eine Versicherung 
Torpedierunge eingeführt wurde, erwähnt 
werden können. 


veren 


Bei der Vermözensversieherung, unter der Haft 
ptlicht Kredit-. Betriebsverlust- und Rückver 
siceherunz behandelt werden. «dürfte insbesondere 
ie letzte künftie dem Mathematiker Probleme 
stellen. (ilt doch alleemein von ihr, was 
Bereer in seinen 1925 erschienenen Prinzipien 
ter Lebensversicherungstechnik unter 
lieet keinem Zweifel. daß die Lebensrück versiche 
rune, wie sie heute betrieben wird. auch nieht an 
nähernd der Forderung gerecht wird, die ihr ge 
stellte Aufgabe in theoretisch einwandfreier und 
praktisch brauchbarster und billiester Art zu er 
füllen. Die amerikanischen (Aktuare) Versiche 
rungesmathematiker beschäftigten sieh in neuerer 
Zeit stark mit allen auf die Güterversicherung be- 
zürliehen Fraeen. Sie haben sogar einen be 
sonderen Verein. auf (dessen wertvolle Abhanıl 
lungen Burran vor zehn Jahren in einer Abhanil 
lune .Beiträge zur Theorie und Praxis der Ver 
sicherung außer der Lebensversicherung“ (Zeit 
sehrift für die gesamte Versicherungswissenschaft 
1922) die europäischen Aktuare aufmerksam ge 
macht hat. Bei der von Manes angegebenen 
amerikanischen Literatur zum 2. Band hätten wohl 
auch diese Abhandlungen genannt werden können. 
Der 3. Band ist «der Personenversicherungz zewirl 
met mit den Hauptteilen Individualversicherung 
und Sozialversieherunge. ‚Jene zerfällt in Lebens- 
versicherung. Krankenversicherung, Invalidenver 
sieherune. Unfallversicherung. Diese behandelt 
die Sozialversicherunge in ihrer Gesamtheit. dann 
auch Krankenversicherung, Invalidenversicherung 
und Unfallversicherunge und schließlich noch die 
Arbeitslosenversicherung. wobei die bis Ende 1931 
bekannt zewordenen Tatsachen. Gesetze und Stati 
stiken berücksichtiet wurden, 


saet: „Es 
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Kin 19 Seiten umfassendes Sach- und Personen- 
register erleichtert die Benutzung dieses nun voll- 
endet vorliegenden dreibändigen sehr zuverlässi- 
een Werkes. das neben dem eroben von Manes 
herausgegebenen Versicherungslexikon allen Ma- 
thematikern, die irgendwie mit Versicherung zu 
tun haben, sehr empfohlen wird. Auf S. 255 des 
3. Bandes ist in Absatz 3 Zeile 2 ein Druckfehler 
zu verbessern: tausend statt zehntausend. 


W.Lorevy, Leipzig. 


Prof. HENRY SCHULTZ, The University of 
Chicaeo, Der Sinn der statistischen 


Nachfragekurven,  Veröffentlichungen der 
Frankfurter Gesellschaft für Konjunkturforschung. 
Heft 10. Verla®e Sehroeder. Bonn 1930. ”>. 
l’reis brosceh. 5.50 M. 


Der bekannte Satz. dab (der Preis eine Funktion 
tes Verhältnisses von Nachfrage und Angebot sei. 
ist schon deswegen unbestimmt, weil beide Größen 
ihrerseits vom Preis (der betreffenden Ware wie 
aller möglicherweise konkurrierenden abhängen. 
ls gilt diesen Satz zu präzisieren. 

Unter Anzebots- und Nachfragekurven versteht 
man mit dem Preis wachsende bzw. fallende Funk 
tionen. welche zleitend bzw. fest genannt wer:ilen. 
je nachdem sie zudem von der Zeit abhängen oder 
nicht. Für einen bestimmten Zeitpunkt gilt dann 
der Preis, an «dem sich die beiden Kurven 
schneiden. Unter Elastizität des Angebots bzw. 
der Nachfrage versteht man den Differential- 
quotienten «des Logarithmus «dieser Größen nach 
dem Logarithmus (les Preises. 

Diese theoretischen Funktionen können aus den 
Statistiken (der produzierten und abgesetzten 
Menzren und ihrer Preise nur dann abgeleitet wer- 
den, falls der Einfluß der störenden Faktoren. ins 
besondere «ler zeitlichen Veränderunz des allge 
meinen Preisniveaus und der Kaufkraft sowie (der 
Krsatzwaren eliminiert werden kann. 

Zur Ermittlung der Nachfragrekurve 
Jahresverbrauch der Vereinigten 
Zucker, also einem Massenpro4dukt ohne Ersatz 
artikel. dem durehsehnittlicehen Großhandelspreis 
für «die Jahre 18090 bis 1914 eerenübergestellt. So 
betrachtet. weisen Preis und Konsum keine Korre 
Istion auf. Es wird angeensmmen. daß (die 
störenden Faktoren langfristie sind und daher 
dureh Ausschaltung des Zeiteinflusses beseitigt 
werden können. Aus «dem Streuunesdiaeramm. 
welches die Preise dem Verbrauch zerenüberstellt. 
ereibt sieh. (lab fallende Preise im alleemeinen mit 
steirendem Konsum verbunden sind. Um den Einfluß 
der Zeit, mit der hier im allgemeinen die Preise 
steigen, auszuschalten. wird der Mittelwert der zu 
je zwei aufeinanderfolgenden Jahren gehörigen 
Differenzenquotienten bereehnet, Eine Gerade mit 
(dieser Neigung und dureh das arithmetische Mittel 
der Preise und «des Konsums. also eine Art Regres 
sionslinie, gibt die Abhängirkeit des Konsums von 
den Preisen unabhängeiz von der Zeit: Eine Er- 
höhunz des jährlichen Durehsehnittspreises um 
I Cent bringt eeteris paribus eine Verminderung 
des Verbrauchs um 217000 Tonnen mit sich. Die 
Abweiehungen der beobachteten Punkte des Streu 
unesdiaegramms von der Geraden werden dann als 
quadratische Funktion der Zeit dargestellt. 

Das gleiche Verfahren wird für die Beziehung 
zwischen dem Verbrauch pro Kopf und (dem durch 
eine Indexziffer dividierten Preis. also dem Real- 
preis, vorgenommen. Diese Zahlen weisen eine 
hohe nerative Korrelation auf. die allerdines zum 
Teil auf dem tendenziellen Steigen des Verbrauchs 
pro Kopf und Fallen der Realpreise beruht. Hier 
bei wird der Verbrauch als lineare Funktion «es 
Preises und der Zeit «dargestellt. wobei eine 
zyklische Komponente übriebleibt. Eine Erhöhung 


wird der 
Staaten an 


des Realpreises um 1 Cent führt ceteris paribus 
zu einer Abnahme des Jahreskonsums pro Kopf 
um 5.83 Pfund. Neben dieser Darstellung wird in 
beiden Fällen noch eine lineare Funktion zwischen 
dem Logarithmus des Preises, des Verbrauchs un 
der Zeit eingeführt. welche die Beobachtungen bei 
nahe ebensogut wiedergibt, jedoch bei Extra 
polation für «die Jahre 1915 bis 1925 auf schlech 
tere Werte führt. 

Irgzendwelche theoretischen Gründe für die eine 
oder andere Darstellung werden nicht vorgebracht. 


So krankt diese Untersuchung am allgemeinen 
Mangel aller Konjunkturforschung: wir wissen 


nicht. woher der Trend kommt, wie er aussieht 
und für wie lang er gilt. Daher bieten auch die 
Korrelationskoeffizienten zwischen «den Beobach 
tungen und ihren analytischen Darstellungen. 
deren hohen Werte für die beobachtete Zeitstrecke. 
weil selbstverständlich. bedeutuneslos sind. für 
die Extrapolation keinerlei Sicherheit. 

E. J. Gumbel (Heidelberg). 256 


OSKAR ANDERSON, Prof. an (der Handelshoch 
schule in Varna. Die Korrelationsreeh 
nune in der Konjunkturforschuneg, 
Kin Beitrag zur Analyse von Zeitreihen. Ver 
öffentliehungen der Frankfurter Gesellschaft füı 


Konjunkturforsehung. Heft 4. Verlag Kurt 
Schroeder, Bonn 1929. 141 S. Preis 8 RM. 


Der Autor gibt zunächst eine allzemein ver 
ständliche, an Tschuprows Unterscheidung von 
empirischen und apriorischen Größen angelehnte 
Kinführune in die Korrelationstheorie. 

Diese ist ursprünglich für zeitlich ungeordnete 
Beobachtungen aufgestellt. Aber bei allen wirt 
schaftlichen Reihen kommt zu der uns inter- 
essierenden zeeeenseitieen Verbundenheit die zeit 
liche Abhängiekeit der Daten hinzu. 

Um die Korrelationstheorie hier anwenden zu 
können, muß man ‘den Einfluß der Zeit aus 
schalten. An Stelle «ler üblichen Methode des 
Trend wählt der Autor eine formal abweichende 
Methode. die der  Differenzenbildung. Jedes 
Glied X, der Zeitreihe wird als zufällige Variable 
mit einem endlichen Erwartungswert m; unıd 
einem Fehler z; aufzefaßbt. Man bildet mun die 
endlichen Differenzen bis zur A-ten Ordnung so 
lange bis Ak m; Ak ır;. 

Das dritte Kapitel gibt Kritiken dafür, wann 
die Komponenten Ak m; verschwunden sind uni 
untersucht, wie aus den Werten Akr; die m; 
und .r; zu berechnen sind. Das Verfahren be- 
steht in folgendem: Von einer gegebenen Zeit 
reihe bilde man die aufeinander folgenden enıdl 
lichen Differenzen und berechne für jede dieser 


Reihen die Größen 0,7 == u (7) wobei 17° 


die zugehörige Streuung, Die aufeinander folgen- 
den Differenzen 07,?-+- 1—07,° vergleicht man mit 
ihren in einer Tabelle zerebenen Streuungen. 
Falls sie etwa den dreifachen Wert dieser Größen 
nieht überschreiten. ist es berechtigt anzunehmen. 
daß die Abweichungen der Reihenglieder von 
ihren Erwartungeswerten gegrenseitir unabhängige 
sind und daß letztere bei der erreichten Ord- 
nuneszahl 4 bereits genügend reduziert sind. 
Dieses o, kann als Annäherung an den mittleren 


Fehler des allen Gliedern zemeinsamen  Ver- 
teilunesgesetzes verwendet werden. Das Ver- 


fahren versagt. wenn die Streuungen 0,? mit der 
Ordnungeszahl # mehr als zulässig zunehmen. Die 
Berechnung der Erwartungswerte besteht, wenn 
das Verfahren überhaupt zulässig ist. in der üb- 
liehen mechanischen Ausgleiehune. Die Diffe- 
renzenmethode vermeidet also den Trend explizit 
anzuschreiben und läßt für seine analvtische 
Natur mehr Funktionen zu als üblich. Aber im 
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‚lizit existiert eine solehe Funktion «der Zeit bei 
ieser wie bei der üblichen Methode. Daher gibt 
lie Differenzenmethode keine wesentlich eröbere 
sicherheit. 

Das letzte Kapitel wendet diese Methode auf 
ie  Korrelationstheorie an. Der Koeftizient 
weier Zeitreihen ist ein kompliziertes Gebilde. 
stehend aus dem Koeffizienten zwischen «len 
athematischen Erwartungen und dem Koefti 
ienten zwischen den Fehlern (gleitender Koefti 
ijent). zwei ganz heterogenen Größen, zu denen 
ventuell eine dritte Komponente hinzu kommt. 
vobei die Größenverhältnisse «dieser Komponen 
en unbekannt sind. Das vorgeschlagene Ver 
'ahren besteht darin, neben «len beiden mittleren 
'ehlern auch die Produktsumme für die A-ten 
‚ifferenzen zu berechnen. um hieraus Annähe- 
rungen an den gleitenden Koeffizienten zu er 
‚alten. Aber schon bei einem künstlich kon 
struierten Beispiel, bei dem die apriorischen Maß- 
ahlen bekannt sind. ist «las Resultat, gemessen 
ım mittleren Fehler. problematisch. 

Es handelt sieh bei der vorliegenden Arbeit 
nieht um Anwendung der Korrelationsrecehnung 
auf die Konjunktur, sondern um eine theoretische 
Darstellung der Differenzenmethode und ihre An 
wendung auf Reihen und Korrelationen. Dem 
entsprechend behandeln «die meisten Beispiele die 
üblichen  wahrscheinliehkeitstheoretischen Sche 
mata und nur zwei Beispiele sind der Konjunktur 
theorie entnommen. Vermutlich liegt dies an 
ler Problemstellung selbst. Denn es ist fraglich. 
ob «diese logisch wie rechnerisch komplizierten 
\Methoden überhaupt einen großen ökonomischen 
\nwendungsbereich besitzen. Vorläufig fehlt 
jeder Beweis, daß die Kapitalisten so handeln. 
lab ihre in den systematischen Wirtschaftsreihen 
auftretenden Ergebnisse solehe Funktionen der 
Zeit sind, welche durch Differenzenbildung ver- 
schwinden. Daher ist die wirtschaftliche Bedeu- 
tung der dureh Ausschaltung soleher Funktionen 
übriebleibenden Schwankung bei einer und der 
Zusammenhänge bei mehreren Reihen proble 
matisch. E. J. Gumbel (Heidelberg). 290 


| 
| 


Ferner sind bei der Schriftleitung folgende 
biicher eingegangen (ausführliche Besprechung 
bleibt vorbehalten): 


MLADEN HEGEDUSIC, Das Wesen der 
Materiestrahlung. Neue Anschauungen 
über die physikalischen Gesetze umd Erschei 
nungen auf Grund (des dreifachen Prineipes der 
Materie. Kine neue Theorie der Gravitation. 
104 S. mit 16 Abb. Verlag Astra Club, Zagreb 
1931. Preis 4.20 M. 


Revue semestrielle des Publications mathe 
matiques vereenigd met Jahrbuch über die Fort 
schritte der Mathematik. uitzereben «door (le 
’reußische Akademie «ler Wissenschaften met 
bijzondere medewerkinge van het Wiskundig Ge 
nootschap te Amsterdam, Redactie Georg FEIGL, 
Hendrik de VRIES. Revue semestrielle 
les Publications mathematiques. 
Deel 37. Jahrgang 1932. Verlag Walter de Gruyter 
x Co. Berlin und Leipzig 1932. Preis Jahrgang 
1932 12 M. 


Dipl.-Ing. TH. BELJAKOW, Prof. an der Tech- 
nischen Hochschule für Baufach. Charkow. Das 
’rinzip der „virtuellen“ Biegelinie. 
alsGrundlagezurlLösung des Knick- 
problems. 
Technik, Charkow 1932. 71 SS. mit 26 Abb. im 
Text, 


Ukrainischer  Staatsverlag für 


H. KREY f. Prof. Dr.-Ing. e. h.. Oberregierungs 
und Baurat. Erddrucek. Erdwiderstand 
und Tragzfähiekeit des Baugrundes. 
(Gesichtspunkte für die Berechnung. Praktische 
Beispiele umd 118 Erddrucktabellen. 4. durch 


earbeitete und erweiterte Aufl. von J. EHREN 


BERG, Regierungs- und Baurat, Berlin. Berlin 
1932. Wilh. Ernst & Sohn. 345 S. m. 239 Abb. 
Preis geb. 26 M. 


H. W. TURNBULL, M. A. F.RS. and 
X, CU, AITKEN, D.Se.. An Introduetion to 
the Theory of Canonical Matrices. 
Blackie & Sons Ltd., London und Glasgow 1932. 
XI rt 192 8. Preis 17/6 sh. 


RICHARD DEDEKIND, Gesammelte ma 
thematische Werke, Herausgereben von 
tobert FRICKE +7 in Braunschweig, Emmv NOE 
THER in Göttingen und Öystein ORE in New 
Haven. 3. Bd. Braunschweize 1932. Verlag Vieweg 
& Sohn. 508 S. Preis geh. 41.40 M. geb. 43.65 M. 


WALTER DUBISLAV, a. 0. Professor an «der 
Teehn. Hochsehule Berlin. Die Philosophie 
der Mathematik in der Gegenwart. 
Heft 13 der Philosoph. Forschungsberiehte. Berlin 
1932. Junker & Dünnhaupt Vle. VIIE #88 S. Preis 
brosch. 3.80 M. 


Revue d’Acoustique. Schriftleitung: J. BRIIL 
LOUIN,. L. BRILLOUIN, R. LUCAS, A, MAR 
CELIN. Bd. I (1932) 5 Hefte. Preis 125 Fr, 
spätere Jahrgänge. 6 Hefte. 150 Fr. Les | 
Universitaires de France, Paris. 


.) 
’yr 


CSS, 


H. v. MANGOLDT’s Kinführunge in die 
Höhere Mathematik. Für Studierende uni 
zum Selbststudium. Vollst. neu bearbeitet und er 
weitert von Konrad KNOPP. ord. Prof. der 
Mathematik an (der Universität Tübingen. 2. Bi.: 
Differentialreehnung, unendliche Reihen. Elemente 
der Differentialeeometrie und der Funktionen 
theorie. 6. Aufl, Leipzig 1932, Verlag von S. Hirzel. 
XV +634 S. mit 108 Fig. im Text. 


Dr. L. BIEBERBACH, ord. Prof. an (der Univer 
sität jerlin. Analytische Geometrie. 
2, Aufl. (Teubners math. Leitfäden, Bd. 29.) Leip 
zie u. Berlin 1932, Verla&g von B. G. Teubner. IV 


.. 


+141 S. mit 44 Fig. im Text. Preis kart. 5.95 M. 

Dr. ALFRED BARNECK, Stwdienrat u. Priv. 
Doz. an der Teehn. Hochschule zu Berlin. Die 
Grundlagen unserer Zeitrechnung. 
2. Aufl. (Math.-Phys. Bibliothek. Bd. 29.) Leipzig 
u, Berlin 1932. Verlag B. G. Teubner. 49 S. mit 
% Fire. im Text. Preis kart. 1.08 M. 


Veröffentlichungen des Forschungsinstitutes der 
Rhön-Rossitten-Gesellschaft e. V. Nr. 5 Jahrbuch 
1930 und Abhandlungen, Herausgeg. von Direktor 
Dr. Walter GEORGII. Prof. f. Flugzmeteorologie 
an der Techn. Hochschule Darmstadt. München 
und Berlin 1932. Verlag von R. Oldenbourg. AIV 
HF 87 8. mit 164 Abb. Preis geh. 10 M. 


Dr. ERICH KAMKE, a. 0. Prof. an (ler Univer 
sität Tübingen. Einführung in die Wahr 


scheinlicehkeitstheorie. Leipzig 1932, 


Verlag& von S. Hirzel. VII + 182 S. m. 2 Abb. Preis 
geb. 11.50 M. 
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Ihre Ursache und ihre Behebung in exakter Dar 
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KURT REIDEMEISTER, 0. ö. Prof. a. d, Uni 
versität Königsberg i. Pr. Einführung in 
diekombinatorische Topologie, Bd. 86 
aus „Die Wissenschaft”, Sammlge. von Einzeldar 
stelluneen aus den Gebieten der Naturwissenschaft 
und der Technik. Hrsg. v. W. Westphal. X 
>)9 S, Braunschweize 1932, Verlage Frieilr. Vieweg 
& Sohn Akt.-Ges. Preis geb, 19 M. 

Dr. HANS FALKENHAGEN, Prof. an der Uni 
versität Köln. Elektrolyte. Leipzig 1932. 
Verlae S, Hirzel. XVI + 346 S. mit 104 Abb. Preis 
eb. 24,80 M. 

Dr.-Ing. ©. PFLEIDERER, !’rof, an der Techn. 
Hochsehule Braunschweig. Die Kreiselpum 
pen, 2. verb. Autl. Berlin 1932. Verlag Julius 
Springer. IX + 454 S. mit 335 Textabb. Preis geb. 


>2)50 M. 


Jahrbuch über die Fortschritte der Mathematik. 
Hrsg, v. d. Preuß. Akwlemie der Wissenschaften. 
Sehriftleiter Georze FEIGL. Bd. 561. Jahrg. 1950. 
Heft 1. 160 8. Berlin u. Leipzig 1932. Verlag 
Walter de Gruyter & Co. Preis 10.80 M. 

WILHELM STIEDA, Johann Albrecht 
Kuler in seinen Briefen 1766-1790. Ein 
Beitra@ zur Geschichte der Kais. Akademie der 
Wissenschaften in St. Petersburg. Leipzig 1932. 
Verlae@ S, Hirzel. 43 S. Preis eeh. 1.60 M. 


REINHOLD MÜLLER, Dr. phil.. Dr. rer, techn. 
h. e.. Prof. i. R. an der Techn. Hochschule Darm- 
stadt. Einführungin die Theoretisch: 
Kinematik. insbesondere für Stwlierende des 
Maschinenbaues. der Elektrotechnik und der Ma 
thematik. Berlin 1932. Verlag Julius Springer, 
IV+124 S. mit 137 Abb. Preis kart. 6.80 M. 


Dr. phil. W, FELGENTRAEGER, Überreg.-Rai 
a. D. Feine Waagen, Wägungen und 
Gewichte. 2. verm. u. umgearb. Aufl. von 


„Theorie. Konstruktion und Gebraueh «er feineren 


Hebelwaage“. Berlin 1932. Verlag Julius Springer, 
V1+3208 Ss. m. 117 Textabb. Preis zeb. 26 M. 


’rof. Dr. E. J. GUMBEL, Das Zufalls 
vesetz des Sterbens. 12. Ereänzungsheft 
des Dtsch. Statist. Zentralblatts. Leipzig 1932. 
Verlage B. G. Teubner. 65 8. Preis kart. 5M. 


Dr.-Ing. e. h., Dr. phil. L. PRANDTL, o. Prof. 
an der Universität Göttingen, und Dipl.-Ing. Dr. 
phil. A. BETZ, a. 0. Prof. an der Universität Göt- 
tinzen. Ergebnisse der Aerodvnami 
schen Versuchsanstalt zu Göttingen 
aneerliedert dem Kaiser Wilhelm-Institut für 
Strömungesforsehung). IV, Lieferung. Mit 234 Abb. 
und 127 Zahlentafeln im Text. München u. Berlin 
1932. Verlage von R. Oldenbourg, 153 >. Preis 


eh. 11.80 M. 


NACHRICHTEN 


Internationaler Mathematikerkongreß in Zürich. 

In der Zeit vom 4. bis 12. September d. J. 
fand ein Internationaler Mathematikerkongreb in 
Zürich unter Beteiligung von mehr als vierzig 


Nationen statt. Von deutscher Seite war insbe 
somdere auch die angewandte Mathematik stark 
vertreten. Aus der «eroben Zahl der deutschen 


Vorträee in den der anzewandten Mathematik ge 
wilmeten Sektionen seien die folgenden ange- 
führt: Riebesell (Hamburg). Was folet aus 
dem Misesschen Wahrscheinliehkeitsbegriff für die 
Versicherungsmathematik? v. Mises (Berlin). 
"raren der Wahrscheinlichkeitsreehnune: Pol 
laczek-Geiringer (Berlin). Bemerkungen 
zur Korrelationstheorie: G. Scehulz (Berlin). 
Über das Problem «ler Markoffschen Ketten: 
Sternberg (Breslau), Wahrscheinlichkeit un«d 
rfahrune: H. Sehmidt (Köthen), Zur Statik 
umd Dynamik elastischer Platten: Lotz (Göt 
tineen). Neure Probleme der Tragtlügeltheorie: 
Vietoris (Innsbruck). Zur graphischen Inte 
ration eewöhnlicher Differentialgleichungen: 
Cauer (Göttingen) Zur Theorie der Wechsel- 
stromsehaltungen: Korn (Berlin,  Mechani- 
sierung «der Wellenmechanik und der Quanten- 
theorie: W. Mueller (Prag). Laminare Ausbrei 
tungserscheinungen in Flüssiekeiten: Scehlich 
ting (Göttingen). Über «lie Stabilität der Couette- 
Strömung. 

Der nächste Kongreß soll im Jahre 1936 in Oslo 
stattfinden. 


Versicherungswissenschaftliche Sachverständigen- 
prüfung. 

Durch einen Erlaß «der preußischen Minister 
für Wissenschaft. Kunst und Volksbildung sowie 
für Handel und Gewerbe ist an «ler Universität. 
Technischen Hochschule und Handelshochsehule 
in Berlin eine  versieherungswissenschaftliche 
Sachverständigenprüfung eingeführt worden. Die 
Prüfung ist entweder eine wirtschaftlich-juri 
stische oder eine mathematische. Voraussetzung 
für die Zulassung ist die vollständige Absol 
vierung eines Hochschulstudiums einschließlich 
normaler Abschlußprüfung, also für Mathematiker 


der Dr. phil... Dr.-Ing... der Dipl.-Ing. oder der 
Stwdienreferendar. Die Prüfung, «lie in einen 
schriftlichen umd einen mündlichen Teil zerfällt. 
erstreckt sich in beiden Klassen auf allgemeine 
umd besondere Versicherungswirtschaft. ferner für 
die Bewerber der mathematischen Klasse auf Ver 
sicherungesmathematik und je nach Wahl des Be 
werbers auf eines der drei Fächer: Wahrschein 
liehkeitsreehnung, mathematische Statistik. Wirt- 
schafts- und Finanzmathematik. Die Prüfung 
wird vor einem Prüfungsausschub abgelegt, dem 
Dozenten «er drei genannten Hochschulen ange- 
hören. Der Vorsitzende. an den alle Meldungen 
zu riehten sind. ist zur Zeit Herr Dr, Büchsel. 
Verwaltungsdirektor bei «der Universität. stell 
vertretender Vorsitzender Prof. Dr. v. Mises. 


Persönliches. 


Am 25. Juli konnte Hr. 
laus Jolles seinen 75. Geburtstag, am Tage 
darauf sein SV0jähriges Doktor-Jubiläum feiern. 
Herr Jolles. der seit 1907 o. Professor der Dar- 
stellenden Geometrie an der Technischen Hoch- 
schule Charlottenburg ist. gehört zu den nam- 
haftesten Vertretern der synthetischen Geometrie. 

Am 7. August vollendete Hr. Hofrat Prof. Dr. 
PhilippForcehheimer sein 80. Lebensjahr. Er 
ist dureh seine Arbeiten auf dem Gebiete der tech- 
nischen Hydromechanik, namentlich durch sein 
l,ehrbuch der Hydraulik. «las kürzlich in 3. Auf- 
lage erschienen ist. weithin bekannt. 

In Stuttgart berine am 28. August Hr. Prof. 
Dr. R. Mehmke. der lange Jahre hindurch als 
Mathematiker und Geometer an der Technischen 
Hochschule in >tuttgart gewirkt hat. ein ge- 
sehätzter Mitarbeiter unserer Zeitschrift, seinen 
75. Geburtstag. Herr Mehmke war auch Heraus- 
eeber der bis 1915 erschienenen „Zeitschrift für 
Mathematik und Physik“. die in mancher Hinsicht 
Vorläufer unserer Zeitschrift war. 

Hr. Dr. Robert Sauer. bisher Privatdozent 
an der Technischen Hochschule München, ist zum 
n. b. a. o. Professor an der Technischen Hoch- 
schule Aachen ernannt worden. 


Gehr. Prof. Dr. Stanis- 
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Die Spannungsverteilung in Schwei- 
ungen. In dem Aufsatz von W. Hovgaard: 
)jeSpannungsverteilungin Schwei- 
ınzen“t) wird auf amerikanische Untersuchun 
n und Versuche aufmerksam zzemacht, die ver 
ehrte Klarheit in den Spannungszustand von 
ektrisch zeschweißten Kehlnähten bringen sollen. 
diesem Zusammenhang sei auf die Versuche des 
hweizerischen Vereins von Dampfkessel-Be 
tzern. Zürich, vorgenommen vom Unterzeiehneten 
0 und 1931, hingewiesen. Die Überlerungen 
nd Ergebnisse sind kurz die folgenden, 
Werden die Spannungen eines Flachstabes 8 
it Trennfuge zwischen den Stabhälften «durch 
ıfzeschweißte Laschen (Flacheisen) L von der 
nen Seite der Trennfuge auf die andere über 
aren. und sind die Laschen nur an den Flanken. 
echt an der Stirn mit «den >Stabhälften S ver- 
ımden, Abb. 1 und 6. so nimmt im Stab 8 (lie 
Spannung gegen die Fuge hin ab. in den Flach 
sen L zu, der Verlauf ist durch AD bzw. BC, 
\bb, 2, gekennzeichnet. Die Längung des Stabes 
‘t von derjenigen der Laschen verschieden, im 
Schnitt 7 ist der Unterschied infolge des Hook - 
chen Gesetzes 
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Abb. 1 bis 3. 


\bh. 1. Stab mit dureh Flankennähte aufgeschweißten 
asehen, das Bleeh ist dureh eine Trennfuge gehälftet. 


Abb. 3. 


Schematischer Spannungsverlauf 

bei Bleeh und Laschen. 

\hbh. 3. Innere relative Verschiebungen Ja, (max. Werte 
Jay und JSa,) bzw. äußere Sag, schematisch. 

Im Schnittpunkt E kann der Stab wegen 05 = 07 
»ingespannt und der Spannungszustand von hier 
us verfolgt werden. Die Ebene «durch E nennen 
vir die „Bezugsebene“, weil alle Verschiebungen 
ınd auch die Schubspannungen auf sie Bezug 
iehmen. Teilen wir Laschen und Stab von der 
ezugsebene aus in gleiche Teile /. so ändern 
ich bei der Belastung die ursprünglich kon- 
tanten Teillängen 2 unter Annahme proportio- 
aler Verhältnisse 


nach rechts 


tt(1—+a) t(1+2ua) ...t(1-+ na) 


nach links 
tt(1l—a) £(l 2a) 


!, Zamm Bd. XI (1931). S. 341 bis 348. 





t(l—nu) ) 


S nach rechts 


S nach links 


umeekehrt beim Stab 


tt(il bb) ti 





! ti] r 3) ti] T 25) er t(] T NM D) 
worin @ und 5 abhängige sind von Teillänge, Quer 
schnitt, Last und Elastizitätsmodul, 


Die erößten Verschiebungen Jay. «lie wir innere 
relative Verschiebungen nennen wollen, ergeben 
sich an den Enden der Laschen und bei der Fuge. 
Der Verschiebungesmechanismus kann mit «lem 
jenigeen ineinandergelerter Bürsten mit wenig 
elastischen Borsten und weichen, sehr dehnbaren 


Böden (Kautschuk) verglichen werden. Diese 
Verschiebungen treten als Wirkung des Hook 


schen Gesetzes in die Erscheinung. 


Zu einer zweiten Art der Verschiebung ge 
langen wir dureh die Berücksichtigung der Wir 
kungen der Biegung bei «der Übertragung der 
Spannungen von einem parallel geschichteten 
Körper (Blech) an den andern. Bei kurzem Hebel 
arm (Abstand der Mittelebenen) treten die Schub 
spannungen «des gebogenen Körpers (der Nähte 
hinsichtlich der Bruchgefahr hervor, die Biegungs- 
spannungen zurück. Zu einer Vorstellung g 
langen wir dureh den Vergleich, die Bürstenböden 
seien starr, die Borsten elastisch. Die rel. Ver 
schiebung des einen Blechs hinsichtlich des 
andern, sie sei „äußere rel. Verschiebung“ ze 
nannt,. Jas von Abb. 3. ist dann konstant, unıd 
„war auf der ganzen Länge der Flankennähte. 
Die „gesamte“ rel. Verschiebung ist 


e 


la a. + A 5 KR) 

Diese Überlerungen sind «durch den Versuch 
bestätiet worden. Die Verschiebungen von 
Flacheisen (Laschen) hinsiehtlieh ihrer Unter 
laxen sind durch Dehnungsmessungen (Abb. 4) 
festzestellt worden (an der EKidgenöss. Material 


prüfungsanstalt Zürich). Die eine Schneide des 
Dehnungsmessers wird auf das Flacheisen bei M 
(Abb, 7). die amdere ebenfalls auf das Flacheisen 
umd später auf (die Unterlage aufgesetzt, zu 
diesem Zweck ist diese mit Zapfen versehen. 
Den in Abb. 4 gezeigten Ergebnissen lag ein 
Stab mit Trennfuge zugrunde. Die Dehnungen 
las der Unterlage und Aa, der Lasche werden 
von einer Grundlinie aus aufgetragen, der Unter 
schied Sa bedeutet die örtliche Verschiebung des 
Stabs hinsichtlich der Lasche, er kommt in der 
schraffierten Fläche zum Ausdruck: die ge 
strichelt schraffierte Fläche, Abb. 7, ist «lie näm- 
liche, dargestellt über einer gemeinsamen Grund- 
linie. Abb. 8 zeigt den Spannungsverlauf os der 
Stabhälften in der Ebene E und o, der Laschen 
in D, dureh die Sehnittpunkte ist die Lage der 
bezugsebenen P, und P, bestimmt. Die ge- 
strichelten Flächen. Abb. 5. sind in der äußern 
Form (der Verschiebungsfieur Abb. 3 sehr ähn- 
lich. Abweichungen lassen «darauf schließen, daß 
die Spannungsänderungen in Stab und Lasche 
nicht proportional sind, wie bei Abb. 2 und 3 an 
eenommen, man vergleiche den durch Messung 
festgestellten Spannungsverlauf, Abb. 8. Auch 
die Faktoren «a und £ der Gl. (2) und (3) sind nicht 
konstant. Die kleinsten Ordinaten der Abb. 5 
entsprechen den äußern rel. Verschiebungen 4 @s 
der Abb. 3. die innern rel, Verschiebungen Aaı 
dem Unterschied 4a — Aas. 
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werden. (die rel, 
Verschiebuneen der hinsichtlich der 
aufgeschweißten Laschen den Schubspan 
nuneen der Flankennähte in ihrer Längsrichtung 
verhältniseleich. 


eetroffen 
Stabhälften 
seien 


Die Annahme darf 


BE te Br TE: 


\nnahme beruht insofern auf einer An 
näherung, als im nämliehen Querschnitt durch 
Flankennaht Konstanz der Schubspannungen 


Diese 
sin 


e111e 


voraussetzt, obwohl dies offenbar nieht zutrifft. 
Durch «lie Annahme eelingt es jedoch. sieh ein 
Bild über den Spannungsverlauf in der Längs 


riehtune einer Flankennaht zu machen, worauf es 


zunächst ankommt. 
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Abh. 4 bis 8. 

Abb. 4. Feststellung der gesamten rel. Verschiebungen 
bzw. der Schubspannungen der Flankennähte bei Stab 12 
(mit Fuge) 

\bb.D. Gesamte rel. Verschiebungen, Ja, herausgezeichnet. 


Probestahb. 


Abb. 6 und 7. 


\bb.s, Spannungen von Stab und Lasche. 
Die max. Schubspannungen treten bei Abb. 5 
(zestrichelt schraftierte Flächen) an den Enden 


der Laschen und bei der Fuge «er Stabhälften 
in Erscheinung. Die in Gl. (5) getroffene An 
nahme führt zur Proportion 

an 28 Be (6) 


(): F ermittelt werden. 
Schweißnähte zu über- 


tr kann aus der Gleiehung 
wobei @ die Last. die die 


bedeutet 
ein Mehı 


Hafttläche 
werden 


und # ihre 
Schubspannungen 


haben 
max, 


tragen 
Die 


faches «ler mittleren 
links Tax = 24T 
rechts 7 ET ET. 


IIAN 


Der Faktor von 7 wird je nach Länge und Be 
sehaffenheit der Flankennähte wechseln. Sin 
die Flacheisen (Laschen) auch durch Stirn 
nähte mit dem Untergrund (Blech) verbunden 
so ändert sich der Spannungszustand etwas, wii 
haben zefunden, daß (die größte Schubspannung 


dann bei der Fuge auftritt umd beträgt 
a 
Emax = 1,17€ Dis 23,27. 


Gänzlieh verschieden vom Spannungszustanı 
des Stabes mit Fuge ist der, wenn die Fug: 
auseeschweißt ist. die Erörterung an diesen 


Stelle würde jedoch zu weit führen t). 

Die Ausführungen von Hovzaard und die 
amerikanischen Versuche bestätigen innerhalb ge 
wisser Grenzen das Vorlierende. 


KH. Höhn, Zürich. 


Der ebene Spannungszustand. In seinem 


Aufsatz in Heft 3. 1932. dieser Zeitschrift ent 
wiekelt Herr Prof. Hager ein Näherungsver 
fahren, das weeienet ist. als Ersatz für meine 
in Heft 3. 1930. gegebene strenge Lösung des 
Biegungsproblems (des Balkens auf 2 Stützen zu 
dienen. Es ist jedoch nieht richtig, wenn der 
Verf. auf S. 141 oben behauptet. daß bei meiner 


lösung «die Schubspannungen am oberen Rande 


nicht, wie es sein muß. verschwinden. Man ver 
eleiche dieserhalb Abschnitt 5. Gl. (3b). meines 
Aufsatzes, die für „= + zu null wird, wie 


+) 


auch in der Erörterung der Randbedingungen im 


eleichen Abschnitt ausdrücklich gezeigt. 
Im übrigen kann natürlich die Ausarbeitung 


von Näherungsverfahren nur begrüßt werden. Für 
etwa quadratische Scheiben erweist sich die 
Übereinstimmung als gut, während für schlankere 
Balken die Konvergenz zur Navierschen Lösung 
in Wirklichkeit rascher erfolgt als nach der Nähe 
rune, Für stark gedrungene Tragwände dagegen, 
bei denen nur der untere Teil mitträgt. dürfte die 
Näherunge nur mit großer Vorsicht zu brauchen 
sein. Craemer,. Frankfurt a. M. 295 


Die Behauptung des Herrn Dr. Uraemer bezüg 
lieh der Schubspannungen an den Rändern nach 


seiner Formel ist riehtie. Wenn meine Verein 
fachung der Spannungsfunktionen für quadra- 


tische Tragwände Lösungen gibt, die mit denen 
der Reihenentwieklung gut übereinstimmen, hat sie 
ihren Zweck vollkommen erreicht, denn nach den 


polarisationsoptischen Versuchen von Dr. Bay?) 
kommt bei hohen Tragwänden nur der untere 


quadratische Teil für die Biegungsspannungen in 


Betracht, Harzer. 295a 
ı) Vel. Höhn, „Über den Spannungszustand und die 


Festigkeit von Kehlnähten‘, Braunschweig, Vieweg u. Sohn. 

2) Über den Spannungszustand in hohen Trägern, Stutt 
art 1931. Dr. Petry, Bericht über Scheiben und Schaben. 
Intern. Kongreß für Brücken- und Hochbau 1932. 
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